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12 - Vom Ziele der Wissenschaft - Wie? Das letzte Ziel der Wissenschaft sei, dem Menschen moglichst 
viel Lust und moglichst wenig Unlust zu schaffen? Wie, wenn nun Lust und Unlust so mit einem Stricke 
zusammengekniipft waxen, dass, wer moglichst viel von der einen haben will, auch moglichst viel von 
der andern haben muss,-dass, wer das "Himmelhoch-Jauchzen" lernen will, sich auch fur das "zum- 
Todebetriibt" bereit halten muss? Und so steht es vielleicht! Die Stoiker glaubten wenigstens, dass es 
so stehe, und waren consequent, als sie nach moglichst wenig Lust begehrten, um moglichst wenig Unlust 
vom Leben zu haben (wenn man den Spruch im Munde fiihrte "Der Tugendhafte ist der Gliicklichste" , so 
hatte man in ihm sowohl ein Aushangeschild der Schule fur die grosse Masse, als auch eine casuistische 
Feinheit fur die Feinen). Auch heute noch habt ihr die Wahl: entweder moglichst wenig Unlust, kurz 
Schmerzlosigkeit-und im Grunde diirften Socialisten und Politiker aller Parteien ihren Leuten ehrlicher 
Weise nicht mehr verheissen-oder moglichst viel Unlust als Preis fur das Wachstum einer Fiille von feinen 
und bisher selten gekosteten Liisten und Freuden! Entschliesst ihr euch fur das Erstere, wollt ihr also die 
Schmerzhaftigkeit der Menschen herabdriicken und vermindern, nun, so miisst ihr auch ihre Fahigkeit zur 
Freude herabdriicken und vermindern. In der Tat kann man mit der Wissenschaft das eine wie das andere 
Ziel fordern! Vielleicht ist sie jetzt noch bekannter wegen ihrer Kraft, den Menschen um seine Freuden 
zu bringen, und ihn kalter, statuenhafter, stoischer zu machen. Aber sie konnte auch noch als die grosse 
Schmerzbringerin entdeckt werden!-Und dann wiirde vielleicht zugleich ihre Gegenkraft entdeckt sein, ihr 
ungeheures Vermogen, neue Sternenwelten der Freude aufleuchten zu lassen! 



F. Nietzsche, Die frohliche Wissenschaft ("la gaya scienza"), 1882 



Introduction 



The subject of this thesis is the study of dissipative dynamics and their properties in 
particle physics. 

Physical systems can be considered closed when they are isolated from their envi- 
ronment; otherwise, due to interactions between them, there occur exchanges of energy, 
entropy and matter from system to environment and vice versa. Physical systems are 
termed open in such a case. 

Typically, the physical systems we will be concerned with have few degrees of freedom 
in comparison to the very (infinitely) many of the environment and matter exchanges are 
neglected. Closed quantum systems evolve reversibly in time according to the unitary dy- 
namics generated by the Schrodinger equation; in particular, there cannot be decoherence 
since Hilbert space vectors evolve in time into Hilbert space vectors. 

On the contrary, the time-evolution of open quantum systems is quite complicated due 
to the coupling with the environment. It can be formally derived from the unitary dynam- 
ics of the closed system formed by the open system plus the environment, by eliminating 
(or summing over) the degrees of freedom of the latter. 

The result is a so-called reduced dynamics for the open system alone that does not 
conserve its energy and entropy; the dynamics is thus dissipative and, moreover, does not 
transform pure states into themselves, rather it destroys quantum coherence. 

The reduced dynamics generally contains memory effects, that is it depends not only 
on the initial conditions and the observation time, but on the whole history of the system 
in between. However, if the coupling with the environment is weak, it is possible to 
approximate this dynamics with semigroups of linear transformations in which the memory 
effects have disappeared, namely by Markovian semigroups which constitute the class of 
time-evolutions studied in this thesis 0[2J|3] • 

Open quantum systems are met in many branches of physics; in particular, the ap- 
proach based on semigroups of linear transformations has been applied in several contexts, 
ranging from phenomena of quantum relaxation in magnetic resonance [I], to quantum 
optics jnilElIZllHl) to foundational aspects of quantum mechanics (HIE], to the study of 
non-standard dynamics in elementary particle physics |44U49l[TT7llll8j and of interference 
phenomena with photons or neutrons |10 | I11 | IT2] . 

Since reduced dynamics are decohering, the states of open systems must be described 
by density matrices; their eigenvalues represent probabilities, so that they have to be 
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positive. In order to have a consistent statistical interpretation of the formalism, the 
dynamical maps describing the time-evolution must preserve such positivity; technically 
speaking, they need to be positivity preserving or positive for short. However, when there 
is entanglement between the open system and its surroundings, positivity is no longer 
sufficient to guarantee the consistency of the physical description. A stronger request need 
to be enforced: the so-called complete positivity. Formally, complete positivity assures the 
positivity of states of composite systems in which there is a statistical but not a dynamical 
coupling between the system under study and a generic n-level system. 

Due to its rather abstract definition, complete positivity is sometime considered as a 
mere mathematical property rather than a physically motivated request. Indeed, the liter- 
ature about open quantum systems abounds in reduced dynamics that are not completely 
positive (and some examples will be discussed below) |125|ll26( l27j. 

The physical systems and the discussions presented in this thesis provide some ar- 
guments in favor of complete positivity by showing explicitly the unavoidable problems 
that arise by considering dynamics that are not completely positive. Moreover, since the 
models that will be presented represent concrete systems amenable to actual experiments, 
they offer the possibility of direct tests of this property. 

In Chapter after introducing the basics of positivity and complete positivity for 
generic time-evolutions, we develop the formalism leading to the reduced dynamics of 
open quantum systems; in particular the weak coupling limit and the resulting general 
structure of semigroups. 

In Chapter [21 we restrict our attention to open quantum systems described by two- 
dimensional Hilbert spaces, for the concrete physical systems we shall consider can be 
described in this framework. We introduce the coefficients that characterize general irre- 
versible reduced dynamics fulfilling the Markov property; in order to guarantee either the 
positivity or the complete positivity request, they must satisfy some constraints, expressed 
by sets of inequalities. 

In the next three chapters this formalism is applied to three specific contexts: the study 
of irreversible dynamics of neutral B mesons in a gravitational environment (Chapter E|, 
of a beam of neutrinos crossing a layer of matter (Chapter HJ, and finally of neutrons in 
interferometric experiments in presence of stochastic magnetic fields (Chapter [SJ). 

The analysis of a possible dissipative time-evolution for the neutral B mesons is sug- 
gested by general considerations about the behavior of quantum gravity (or, from a more 
fundamental level, string theory) at the Planck scale. The approach is phenomenological, 
that is the microscopic details of the interaction between system and gravitational bath 
are unspecified; in this sense, the adopted dissipative time-evolution for the B mesons is 
completely general. We have explicitly evaluated the expressions of the decay probabil- 
ities of neutral mesons for some selected decay channels and, by means of them, some 
asymmetries, that is quantities usually employed to evaluate the violation of the discrete 
symmetries (parity, time reversal and charge conjugation). These asymmetries can be di- 
rectly used to fit experimental data. In this respect, entangled B mesons produced in the 
decay of spin-1 resonances are particularly useful: in such a case, the dissipative effects 
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induced by the non-standard dynamics produce characteristic signatures, that, at least in 
principle, can be experimentally probed. Indeed, the experiments involving entangled neu- 
tral B mesons at the so-called B-factories have reached a high level of sensitivity, very close 
to the one required for detecting the non-standard, dissipative phenomena. This is a fur- 
ther motivation for the investigations presented here. Moreover, since complete positivity 
is tightly related to the quantum correlations expressed by entanglement, the B-factories 
offer the possibility of performing experimental tests of this fundamental property. 

The analysis of neutrino interactions and oscillations in matter described in Chapter 
is not purely phenomenological in character, rather it is based on a microscopic model. The 
formalism is nevertheless very general and, provided the Markovian approximation is valid, 
it can be applied to study dissipative phenomena in neutrino physics induced by various 
matter effects, e.g. the one induced by the thermal fluctuations inside the Sun. We have 
computed the modified form of the transition probability between two different flavors; 
for thermal fluctuations in the Sun we find no meaningful changements to the standard 
expression, in agreement with previous results coming from different approaches. 

Finally, in the last chapter, we devote our attention to neutron interferometers and 
discuss how these apparata can be used to study dissipative dynamics, in a controlled envi- 
ronment. The basic idea is to insert the interferometer in an external stochastic magnetic 
field. By suitably adjusting the parameters characterizing the fluctuating external field, 
one can then reproduce the physical situations adopted to justify widely used dynamics 
(e.g. of the Redfield-Bloch type). It turns out that such time-evolutions are not completely 
positive and sometimes not even positive; neutron interferometry would then allow direct 
experimental tests of the inconsistencies that such dynamics necessarily produce. 

Neutron interferometers have been proposed as one of the possible set-ups to study non- 
contextuality in quantum mechanics, by testing suitable Bell's inequalities. In presence of 
dissipative effects leading to loss of quantum coherence, one expects a weakening of the 
effectiveness of these Bell's inequality tests. As a byproduct of the previous investigations, 
we have studied in detail how the expressions of these inequalities change under the action 
of dissipative dynamics. 

In this thesis, elementary particle phenomenology and quantum interferometry are 
shown to provide, under specific conditions, very good laboratories to investigate the 
properties of reduced dynamics of open systems, in particular the notion of complete 
positivity. The motivations about relevance and necessity of complete positivity result 
strengthened by the arguments developed in what follows. 
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Chapter 1 

Open Quantum Systems 



In this chapter we introduce the concept of quantum open systems and their reduced 
time-evolutions. After some basic definitions, we focus our attention on the properties 
of positivity and complete positivity of reduced time-evolutions, stressing their relevance 
and the problems that arise if they are not positive or not completely positive. Complete 
positivity shall be the leit-motiv of this thesis, since it is usually believed that any dy- 
namical map describing the time-evolution of an open quantum system should enjoy this 
property in order to avoid any fallacy in the interpretation of the formalism. However, 
in the literature about open quantum systems many examples of not completely positive 
reduced dynamics are encountered; we will present some of them to show the problems 
they suffer. 

In Section fLSI we explicitly show how, by means of the projector formalism, it is possible 
to extract the dynamics of a subsystem of a larger system whose time-evolution is known. 
The reduced dynamics so obtained is in general rather complex because it contains memory 
effects negligible only when some particular conditions are fulfilled, typically involving the 
weakness of the coupling between system and environment. In view of this, in Section fl. 41 
we introduce the so-called Markovian dynamics, that is the class of time-evolutions in 
which memory terms are made to disappear; they are described by linear maps satisfying 
the semigroup composition law. This is the class of dynamics we shall consider in the 
following chapters. 

A different approach to the time-evolution of an open quantum system is presented in 
Section fT~71 based on a Langevin-type equation. 

In the final part of this chapter some considerations about the long-time behavior of 
states evolving under a Markovian dynamics are included. 

1.1 States, observables and dynamical maps 

In this section we introduce the basic elements to describe a quantum system: the states, 
the observables and their time-evolutions. 



7 



8 



1. OPEN QUANTUM SYSTEMS 



Given a Hilbert space $), we recall the following definitions: 

Definition 1.1 A linear operator A acting on fj is said to be bounded iff || -Alloc < 00 > 
where 

IIAHoo = sup \\Ail>\\, tpeSj. 
IHI<i 

Definition 1.2 A linear operator A acting on fj is said to be trace class iff A £ and 
||A||i = Tr(A^A) 1 / 2 < +oo. 

We introduce two complex linear spaces of operators acting on |14j : 

• B(Sj): the complex linear space of bounded operators, which is complete in the 
operator norm || • H^; 

• T($j): the complex linear space of trace-class operators, which is complete in the 
trace norm || • ||i. 

Remark 1.1 If fj is finite-dimensional, say dim(^) = n, then B(Sj) = T(Sj) = M n (C), 
the algebra ofnxii matrices. If dim(fj) = oo, T(Sj) C B{9)); in particular £>(fj) is an 
algebra whereas T(Sj) is not. 

The physical states of quantum systems with finitely many degrees of freedom can 
always be represented by positive (then Hermitian) trace-class operators with trace one 1 , 
called density matrices or statistical operators: 

peT(Sj), p>0, (p f =p), Tr(p) = l. (1.1) 

They constitute a convex set in T(Sj), denoted by V(Sj), whose extremal points are 
projectors on one-dimensional subspaces; in spectral representation, they read 

P = ^ Aj|et)(ej|; (e»|ej) = Aj ^ 0, } j X i = l. (1.2) 

i i 

Definition 1.3 A state p is said to be pure if the sum in (|1.2jl can be reduced to only one 
term (after a suitable choice of the basis vectors | e^) ) , otherwise it is said to be mixed and it 
represents a statistical convex superposition of pure states. In such a case the eigenvalues 
are the statistical weights of the corresponding eigenprojectors in the mixture. 

Remark 1.2 Density matrices are not uniquely expressible as convex linear combinations 
of other states. For instance any mixture ^ i p,iO~i with generic density matrices g% and 
Pi > 0, ^2iPi = 1) can always be spectralized and written also in the form ()1.2|) 

Pure states are fully characterized as follows. 

1 Or, as we shall shortly see, less than one if the state represents an unstable system 
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Remark 1.3 A state p G T(Sj) is pure iff p 2 = p. 

The system observables are represented by Hermitian operators A = A' belonging to 
B(S)). The expectation value of physical observables A with respect to a state p, in this 
formalism, is given by the trace operation: 

(A) p = Tr(pA) = ]T \iiei\A\ei). (1.3) 

i 

In order to describe the time-evolution of quantum systems, two approaches are pos- 
sible. In the Schrodinger picture a state evolves as p(0) — > p(t) = A t [p(0)] where At is a 
linear operator, At : V($j) — > V(Sj) and i is time; the one-parameter family of linear maps 
{At,t G 9\} is called time-evolution of the system; D\ is the range over which i can vary 
and it depends on the characteristics of the considered dynamics. 

In the Heisenberg picture, we introduce the one-parameter family of linear maps 
{A|,i G producing the time-evolution of the observables: A(0) — ► A(t) = A*[A(0)]. 
A*[ is called the dual map to Af, defined by the duality relation: 

Tr\p(t)A]=Tr\pA(t)]. (1.4) 
We distinguish the following classes of dynamical maps. 

Definition 1.4 A dynamical map is called reversible iff A t admits an inverse 
Otherwise it is called irreversible. 

Definition 1.5 A dynamical map is called trace-preserving iff Tr[p(t)] = Tr[p(0)], Vi G 9t. 

For the physical consistency of the formalism, any dynamical map must be Hermiticity- 
preserving, that is A<[p]^ = A t [p], Vp G V(J~L), Vi G 9t. Instead, reversibility and trace 
preservation may be satisfyed or not, depending on the physical systems we are dealing 
with. 

In particular, trace-preserving maps describe stable systems in which the total prob- 
ability is conserved; on the other hand, the time-evolution of unstable systems (that is 
decaying systems) is modeled by maps that do not preserve the trace of states: the total 
probability has indeed to decrease in time. 

The standard quantum dynamics is given by the time-evolution of stable closed sys- 
tems, that in the Schrodinger picture reads (putting H = 1): 

p(t) = A t [ P m=e-* Ht p(0y Ht (1.5) 

or, in differential form, 



p(t) = -i[H,p(t)}, 



(1.6) 
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where H = is the Hamiltonian describing the system. In the following we denote by 
-^#["] = — "] the Liouville operator which generates of the time-evolution maps (or 
dynamical maps): 

A t = e LHt Vt G R. (1.7) 
The same evolution in the Heisenberg picture reads 

A(t) = A* t [A{0)} = e iHt A{0)e~ iHt , (1.8) 

which is generated by 

A(t) = i[H,A(t)], (1.9) 

where A is any observable of the system. 

In such a case, the set of linear operators {A t ,t G R} (and the set of their dual maps 
{Aj,i G R}) forms a group of transformations. The group structure implies that for each 
At there is an inverse map AT" = A_t, that is the dynamics is reversible; moreover the 
composition law A tl+f2 = A tl o A i2 , Vii,i2 G R is satisfied, as can be seen from 1)1.7]) . The 
invariance of the trace under cyclic permutations produces the trace-preserving property 
for the dynamical maps. Moreover, it is Hermiticity-preserving and it maintains the 
coherence of the states over which it acts: indeed, from 1)1.5)1 . p 2 [t) = p(t). This is 
another expression of the fact that Hilbert space vectors are mapped into themselves: 
|V) |^> =e- %m \^). 

Another example of dynamical map is the one describing a closed unstable system. The 
description of decaying systems is performed introducing a non-Hermitian Hamiltonian H 
in the Weisskopf-Wigner form (2H1, that is composed by a Hermitian part and an anti- 
Hermitian one: 

H = M- % -r-, m = m\ r = r f (1.10) 

with positive M and r characterizing masses and natural widths of decaying states. In 
this case, in the Schrodinger picture, the states evolution reads 



p(t)=A t [p(0)]=e- im p(0)e iHH (1.11) 



or, in differential form, 



pit) = -iHp(t) + ip{t)rf = -i[M, p(t)] - ^{r, p(t)} (1.12) 



leading to the property 

d 



it Tt[p(t)] = -Tr[p(t)r]<0. (1.13) 
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In the Heisenberg picture the time-evolution of the observables is described by equa- 
tions analogous to (jl.lljl and (|1.12j) . 

In this case the range 9\ over which the dynamics makes sense is not the real axis 
R but rather the positive semi-axis R + . Indeed for t < equation ()1.13j) should imply 
Tr[p(i)] > 1 since Tr[p(0)] = 1, and thus the statistical interpretation of the eigenvalues 
of p(t) would be compromised. Since this time-evolution is defined only for positive time, 
it is irreversible; moreover it is not trace-preserving, as can be seen from (|1.13|) . and it 
maintains the coherence of the states, i.e. Hilbert space vectors transform into Hilbert 
space vectors, though not normalized, since H ^ W. 

In the study of open quantum systems we will meet irreversible time-evolutions either 
trace-preserving or not, and decohering, that is causing the transition from pure to mixed 
states. 

We have to discuss the fundamental property of the preservation of positivity of the 
states subject to a dynamical map. For its relevance and its prominent interest we reserve 
to this issue the next section. 

1.2 Positivity and Complete Positivity 

Let us start with the definitions of positivity and complete positivity: 

Definition 1.6 A linear transformation T : B(Sj) — > B(Sj) is said to be positive (or 
positivity-preserving) iff A ^ =>• T(A) ^ 0, \/A G B(Sj). 

Definition 1.7 A linear transformation T : B(Sj) — ► B(Sj) is said to be n-positive iff 
r (g> I n : B{$) (8) fin) B{$) (g> S) n ) is positive, with n € N and I n the identity on the 
n-dimensional Hilbert space Sj n . 

Definition 1.8 A linear transformation T : B(fi) — > B(fi) is said to be completely positive 
(CP) iff it is n-positive Vn G N. 

Remark 1.4 Complete positivity is a stronger property than positivity which amounts 
to 1-positivity. 

Remark 1.5 Historically, the notion of complete positivity was introduced for linear 
transformations acting on bounded operators, that is, in a dynamical context, not for 
the time-evolution of states, At, but for its dual Ajf, the time-evolution of operators. In 
this thesis, the notion of complete positivity will be referred to dynamical maps on states in 
the sense that a map At : P(Sj) — > V(Sj) is CP if and only if its dual map Ajf : B($j) —> B($j) 
is CP. 

In the context of dynamical maps describing generic time-evolutions of states, the 
property of positivity is required in order to preserve at any time the positivity of the 
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eigenvalues of the statistical operator representing the system. As previously stated, these 
eigenvalues represent the weights of a statistical superposition; then they must be positive 
in order to guarantee the correct probabilistic interpretation of the formalism. For any 
non-positive time-evolution Af there necessarily exists a positive initial state p(0) mapped 
by At in a non-positive one, p(t) = At[p(0)], leading to physical inconsistency of the 
formalism. 

It is less intuitive to understand the relevance and necessity of complete positivity for 
time-evolutions of states. Let Af be a generic dynamical map describing the time-evolution 
of a state ps', following the definition of CP, we have to consider the map At <8> I n , where 
I n is the identity in a n-dimensional Hilbert space and n is any natural number. This 
map represents the factorized time-evolution of a composite system T = S + S n , where S 
is the system described by ps and evolving under At while S n is a generic n-level system 
characterized by a trivial Hamiltonian H n = 0. Since the dynamics of the composite 
system has the factorized form At (g) I n there is no dynamical coupling between S and S n . 
Now, the positivity of any state pr, describing the composite system T, must be preserved 
under the action of At<g)I n , but this is not guaranteed by the positivity of A t : a completely 
positive At is needed. 

To clarify this point, we define a fundamental property of the state describing the 
composite system T, tightly connected to the notion of complete positivity. 

Definition 1.9 Given a state px describing a composite system T = S\ + S2, it is called 
separable if it can be written as pt = Yli j c ijPs! ® z 9 ^ w ^h Cjj > 0, Vi, j; otherwise it is 
called entangled. 

It turns out that the request of complete positivity is necessary when the evolution of 
entangled states pr is taken into account. Indeed, any factorized statistical operator px 
of the bipartite system T = S + S n , evolving under the time-evolution At <g) I n , where At 
is positive-definite, is mapped into a positive operator. In fact, 



Then it is the presence of quantum entanglement that asks for complete positivity, 
that is only entangled states of S + S n can lose their positivity under the dynamical map 
At (g) I n when At is not completely positive. We shall shortly come back to this point after 
discussing the structural characterization of completely positive maps. 

The following theorem is due to Kraus JH] who used a result by Stinespring |17j . 
obtained in the general framework of completely positive maps on C*-algebras. 




(1.14) 



with Cjj > 0, so that, since At ^ 0, 



p T (t) = (A t ® l„)[pr(0)] = Ws(0)] ® /4 n (0) > 0. 



(1.15) 
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Theorem 1.1 (Kraus, Stinespring) A linear transformation T : B(S)) — > S(^) is CP iff 
T(B) = Ylk&K A\BAk, yB G S(.£j), where if is a finite or countably infinite set; the 

operators G satisfy X^fce^: ^l^ fc ^ ^* ^he se * {^fc>^ ^ ^} i s no * uniquely 

determined by V. 

Corollary 1.1 A dynamical map A t : V(fi) -> is CP iff A t [p] = J2keK A kP A {, 

Vp G with J2keK A Uk < I 

Corollary 1.2 For a trace-preserving completely positive map the operators in The- 
orem must satisfy J2keK A \ A k = I- 

When S)s is finite-dimensional, say dim (^5) = d, one doesn't need to check all the 
maps At ® In, but can use the following result. 

Theorem 1.2 (Choi |18ll9j ^ Given a finite d-dimensional Hilbert space fj^, a linear trans- 
formation T : B(S)d) — ► B(Sjd) is CP iff it is d-positive. 

Remark 1.6 Note that the standard evolutions for systems, either stable or unstable 
(equations (|1.5j) and (jLlljl ). are in the form of Corollary thus they are completely 
positive. 



Remark 1.7 A simple example of a positive but not completely positive map is given 
by the transposition in a two-dimensional Hilbert space, r : M2(C) — > M2(C). Positivity 
holds for r[p] ^ 0, Vp G V(^H?)\ to see that complete positivity does not hold, consider the 

where lib' 



singlet state p - 
4x4 positive matrix: 



(| ei ) <8> \&-}) — \e%] ® |ei))/v2. p is represented by the 










/ 



(1.16) 



After explicit computation, we find that p T = r (g> I2 is given by 



/ 



1\ 



V 




10 
10 
-1 / 



that is no longer positive, since one of its eigenvalues is 
completely positive. 



(1.17) 



T/2. Consequently, r is not 



Using Theorem 11.21 we stress the tight connection between complete positivity and 
quantum entanglement. 
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Proposition 1.1 If At is not completely positive, there exists an entangled state pr in 
S + Sd whose positivity is not preserved by Af Id- 
Proof. If At is not completely positive, the same holds for its dual A*. Then, because of 
Theorem II. 2| there must be a positive bounded operator A S Md(C) (g) Md(C) such that 
(A| ® has a negative eigenvalue. Let pr be the corresponding eigenprojection, then 

Tr [p T (A* t ® I d )[A]] = Tr [(A t ® I d )[p T ]A] < 0. (1.18) 
Thus, (At Id)[pr] must have negative eigenvalues since A ^ 0. 

To summarize, positivity of a dynamical map At guarantees that the eigenvalues of 
any state ps evolving under At remain positive at any time. Complete positivity of the 
dynamical map At guarantees that the eigenvalues of any entangled state pr of S + S n 
remain positive at any time, where S n is an arbitrary n-level system, possibly very far 
from S, and the coupling between S and S n is statistical, not dynamical. 

Because of the abstractness of the definition of CP, the necessity of positivity of maps 
of the form At <g> I n , and thus of complete positivity of Af, is often refused as a technical 
artifact rather than a physical constraint. 

In fact, the motivations for considering the system of interest S jointly with an ar- 
bitrary n-level system S n , not evolving and without interaction with S, are not obvious. 
All the more so, since this coupling should motivate the special form that, according to 
Corollarv ll.il the time-evolution would inherit. 

In Section 11.61 we will show that complete positivity of At is a necessary and sufficient 
condition to preserve the positivity of factorized maps in the form At <8> Af. This gives 
the opportunity to study the notion of CP in the treatment of the dissipative evolution of 
entangled pairs of neutral mesons, given in Chapter [2j More in general, time-evolutions 
of the form At <8> At are definitely physically plausible and less abstract than At ® I n . 

Nevertheless, in Chapter [3 we shall deal with a physical setup reproducing a factorized 
dynamics in the form At (g> I n and involving entanglement between the degrees of freedom 
of two subsystems: then, it will be possible to show explicitly, even with At 0l n , that lack 
of complete positivity of At produces necessarily physically concrete inconsistencies. 

1.3 Reduced dynamics: the Generalized Master Equation 

In the previous sections we have described the general features characterizing the time- 
evolutions of quantum systems p G V(7i). Now, we proceed to describe the time-evolution 
equations for open systems. 

We denote the system of interest by S, the environment by B and assume that the 
Hamiltonian Ht of the compound system T = S + B is given. Such a Hamiltonian consists 
of contributions H$ and Hb from system and environment alone, plus an interaction term 
Hi : 



Ht = Hs + Hb + Hj. 



(1.19) 
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Hi is usually assumed to be of the linear form 

H I = J2vi®Vi 3 , (1.20) 

3 

with Vg, V B Hermitian operators of S, respectively B. 

We want to describe states and observables belonging to S only (because of lack of 
control over all the degrees of freedom of B). The formalism we follow in order to obtain 
the dynamics of 5 when the dynamics of T is known, consists in the projector technique 
developed independently by Nakajima |21| and Zwanzig [22] (see |23| for a review). There 
are also different approaches (2411 .25 . 26 , 2 f|l28|l29|l3()j to describe the time-evolution of open 
systems; in Section ITTTl we shall present a stochastic approach (125j . 

The Hilbert spaces characterizing S, B and T are such that Sjt = fis ® We are 
interested in observables A G B(Sjs), whose expectation values are: 

(A) s = Tt{A PT ) = Tt s Tv B [(A®I B )p T } = Tr 3 (Aps), (1.21) 

where the statistical operator of S is obtained from the one of T by the partial trace over 
the degrees of freedom of B: 



p s = TT B (p T ). (1.22) 

We are interested in reduced dynamics of the system S alone, that is in linear maps At 
describing the time-evolution ps(0) — ► A t [ps(0)]. In full generality, the state of S at time 
t, denoted by ps{t), is obtained from the unitary evolution of T = S + B, Ut = e~ lHrt , as 
follows: 



p S (t) = Tt B [U tPT (p)U}]. (1.23) 

However, ps(t) is not always expressible as A t [ps(0)], because there is no general 
prescription for extracting ps(0) from the initial global state />r(0). Indeed, there can be 
two different initial global states pj>(0) ^ Pt{Q) that lead to the same reduced state of S, 
p s (0) = p h s (0), and neverthless at time t p s (t) / Ps(t)- In this case, the dynamics A t is 
not properly defined because a state of 5 at time t = would evolve in two distinct states 
at time t. It is thus necessary to define an assignment map $ : V(S)s) V{S)t) giving) 
for each initial state ps(fi), a state of the total system T such that 

Ps (t) = A t [ps(0)) = Tr B (U t $\ps(P)pt), (1-24) 

and the diagram in Figure [L3l commute. Following this procedure, it is possible to define 
the reduced dynamics A^. The conditions usually imposed on the map $ are: 1. 

preserves mixtures: $ Yli ^iPs = ^2% ^i^lPsl'i ^- 3> is consistent: Tib(^[ps]) = PS\ 3. 

is positive for any positive ps- It has been shown that, for a two-dimensional Hilbert 
space Sjs, the only assignment maps satisfying these three conditions for all initial density 
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Ps (0) e V(H S ) 



At 



Ps (t) E V(H S ) 



$[p S (0)] E P(W T ) 



C/^[ P5 (0)]C// E 7>(W T ) 



Figure 1.1: The procedure required in order to define consistently the reduced dynamics A t acting on 
the states of the system 5* using the projector technique. 

matrices are the product maps $[ps] = ps <8) Pb, where pb is a fixed reference state in 
the bath [21] . This factorized state is usually justified in the context of a weak coupling 
between S and B; however there are physical situations in which such a factorized initial 
state cannot be used [25] . 

Remark 1.8 Since Af in (|1.24j) is obtained by the composition of three maps, its proper- 
ties depends on those of the trace over the environment degrees of freedom, of the standard 
evolution and of the assignment map <I>. If the choice $[ps] = Ps ® Pb is assumed, the 
reduced dynamics At is completely positive since it is the composition of maps that are 
CP. However, in some cases the initial state has not the form of a product state, namely, 
when there exist initial correlations between S and B [241125] . In [21] an argument against 
the complete positivity of a reduced dynamics has been discussed. It has been argued that, 
with correlated initial conditions, dynamical maps At can be properly defined only on a 
subset of initial system states and the linear extension of these maps on all possible initial 
states may not be physically realizable and may not even be positive. In this case the 
condition 3 is not fulfilled; an example of this procedure is given by the assignment map 
proposed in [26] . In literature many examples of not completely positive dynamical maps 
are encountered |2"7ll28U29ll3T?ll3*T] : however, in [32] it has been shown that it is possible in 
any case to preserve the complete positivity of At, giving up either the condition 1 or 2; 
it has also been stressed that, beyond the weak coupling regime, there exists no unique 
definition of the reduced dynamics At for the system S. 

Let us now present the Nakajima-Zwanzig procedure. We assume the dynamics of the 
total system T to be given; furthermore, to be as general as possible, we do not consider the 
latter to be necessarily closed. Namely, we consider, for the states pt of T, time-evolution 
equations of the form 



where Lj- = L$ + Lb + Li is the generator of the time-evolution of pt', only when T 
is a closed system, Lt is proportional to the commutator with the total Hamiltonian 
Ht = Hs + Hb + Hi, see equation (|1.6|l . 



p T (t) = L T [p T (t)] 



(1.25) 
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We define the two projectors P, Q : V(S)t) V($)t) that map states of T into them- 
selves as follows: 

P[ PT ] = Tt b (pt)®Pb = PS®Pb; Q = I-P. (1-26) 

They depend on a normalized reference state ps of the bath, Tr b(pb) = 1> that one 
usually choose to be invariant under the dynamics of B : Lb[pb\ = . The following 
manipulations are independent from the choice of this state. Inserting in equation 1)1.25(1 
p T (t) = P[pr(t)] + Q[p T (t)}, we obtain 

P[p T (t)} = PL T P[p T (t)]+PL T Q[p T (t)}, 

(1.27) 

Q[f>At)} = QL T P[p T (t)]+QL T Q[ PT (t)}. 
Formally, solving the second of these equations, one gets 

Q[p T (t)] = e QLTt Q[ PT (0)] + f* tie QL ^QL T P\p T {t - t')\. (1.28) 

Jo 

Then upon inserting (|1.28|) in the first equation in Q1.27JI . we find a closed equation for 
the relevant part P[pT(t)]; after a partial trace over B we are left with the equation of 
motion of S alone, a so-called generalized master equation (GME): 

p s (t) = L e s ff [p s (t)} + f dt'K{t')\p s (t - t')] +l(t) (1.29) 
Jo 

with 

tf/f^Ls + TTBiLilpB]); 

lC(t)[Ps] = TT B [L ie ^ LTt QL T {ps ® PB ]); (1.30) 
l(t) = Tr B (L ie QLTt Q[p T (0)])- 

These expressions have been obtained using that PLs^bQ = 0. 

The GME is the equation we were looking for; it describes in general the behavior of 
an open system S because no assumptions have been made about the nature of the system 
and of the bath, and the interaction is completely arbitrary; we have only assumed the 
existence of a reference bath-equilibrium state ps in V(S)b)- Notice that in absence of 
interaction between S and B (i.e. if Lj = 0) this equation reduces to the time-evolution 
for S isolated from B, that is p s (t) = Ls[ps(t)]; in fact L e J^ = L$, IC(t) = and Z(t) = 0. 

Let us now describe the various terms appearing in the GME. L e J* generates the 
effective time-evolution in p$, redefined by Tr b{Li[pb\) which represents a mean field 
term. KL{t) is the kernel of the convolution term, that is the memory term depending on 
the whole history of ps- The inhomogeneity 2(t) is related to initial conditions: if they 
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are factorized, />r(0) = ps(0) ® p B , then Qpx(0) = and consequently Z(i) = 0. In the 

following we will assume this condition to be fulfilled. 

The reduced dynamics Aj is well defined and it is completely positive. 

If T is a closed system it is possible to simplify (|1.29|) and IJl.MOj) . In such a case the 
operators L are proportional to the commutators with the Hamiltonians and the following 
relations hold: [P,L S ] = 0; PL B = L B P = 0. Then the GME becomes: 



p s (t)=L e / f [p s (t)]+ / dtX(t')[p s (t-t% 



LflPs] 

mips] 



: -i[H s + Tr B (H I p B ),ps], 
TrB(Lie QLTQt QL J \p s ®p B } 



(1.31) 



Let us focus our attention on cases in which the interaction between S and B is 
weak. We introduce a small dimensionless parameter e in the interaction Hamiltonian; 
consequently Ht = H$ + H B + eHj. In this way we can write a perturbative expansion 
in e for the kernel K(t) in equation 1)1.31(1 |40j : 



with ^ t n ^ 



JC(t) 
. ^ ti ^ t, and 



iC (t) + X> n 



n=l 



dt n K n (t\ti ...t n ) 



(1.32) 



/C (t)[p 5 ] = Tx B (L}(t)QL}[p s ® PB ]), 

^(t|ti . . . i n )[ps] = Tr B (LKt)Q4(^i)Q • • • QL}(t n )QL}[p s ® p B ]) 

where 

L}(tj)[-] = -*[flj(*i), •] = -z[e i(Hs+//B) '^ /e - i( ^ s+Hs) ^, •] (1.34) 

is the Liouville generator in the interaction picture. KL n {t\t\ . . . t n ) contains the multi-time 
correlation function of the bath 



/',.,;...,.. it- h,... t n ) = Tt B (e iHBt V^e- iHBt e iHBt W^e~ iH ^ . . . V^Pb), (1-35) 
where, as alredy mentioned in (11.20(1 . we have assumed the interaction Hamiltonian to be 

Usually, one proceeds with a Born approximation of the GME, namely with neglecting 
in the kernel expansion (|1.32|) all the terms except /Co(i); in this case only two-time 
correlation functions appear. 
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1.4 The Markovian approximation. The Weak Coupling 
Limit 

As we have seen, the GME obtained in the previous section, eq. Q1.31JI (or, more in general, 
eq. (jl.29[) ) contains a memory term that depends on the kernel KL{t). If KL{t) decays rather 
fastly, it is possible to extract the density matrix from the convolution integral, then 
obtaining an equation without memory term [Tl l33ll3"Ij : in such a case it is possible to 
introduce a time-independent generator of the reduced dynamics. 

Definition 1.10 The dynamical map {Aj,t G 9^} is said to be Markovian iff there exists 
a linear bounded operator L on the states space V(7i) which is time- independent, such 
that At = e Lt Vi £ IK; L is the generator of the dynamical map. 

We want to extract a Markovian time-evolution from the GME presented in the pre- 
vious section: in this way we obtain a more manageable time-evolution for the system 
considered. 

Remark 1.9 A Markovian time-evolution {At,t E £R} satisfies the semigroup composition 
law A^ 1+ t 2 = A^j o At 2 , Vti,t2 £ £H) as can be seen from Definition 11.101 

As alredy stated, in order to perform a Markovian approximation of the dynam- 
ics (jl.29|) or (|1.31|) we need the decay time of the kernel KL{t) to be much shorter than 
the typical evolution time of the system S. We find this situation when the interaction 
Hamiltonian can be parametrized by a suitable coefficient e which is then let to vanish: 
Hj — ► eHj. The total Hamiltonian is thus scaled as Ht — > Hf^: 

H £ T = H S + H B + eH u (1.36) 
and the corresponding dynamical map is denoted by A|: 

Af [p s (0)] = Tr B (Uf Ps (0) ® PB U^), (1.37) 

where Uf = e - *^ 4 . 

As we shall shortly see, for small e the reduced dynamics A| becomes Markovian |35j 
over times scaling as e~ 2 : this procedure is called Weak Coupling Limit. 

Since the memory term is of order e 2 (being quadratic in the interaction Hamiltonian), 
the evolution has to be studied over a time-scale of order e~ 2 . It is then appropriate to 
rescale the time as t £ = e~ 2 t; in other words, the evolution is observed over a time long 
enough so that the memory effects can be neglected. On this time scale the convolution 
term in the GME heuristically behaves as: 

i e / r+oo \ 

1C{t')[ps{t £ - t')]dt> « l J IC(t')dt'j [ Ps (t £ )} = K[p s (t s )} (1.38) 
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and thus the GME becomes 



p s (t) = L e s ff [ Ps (t)} + e 2 f K{t')\p s (t - t')]dt' « (L e s ff + e 2 K)[p s {t)l (1.39) 

J o 

whose solution is 

p s (t) = A^ s (0)] = e^+^^W]. (1.40) 

Thus L = ifg^ + e 2 K is the time-independent generator, and the time-evolution Af 
is Markovian. This is the basic idea underlying any Markovian approximation of the 
generalized master equation; there are however many ways to perform the Markov limit, 
leading to Markovian reduced dynamics with different properties |33] , 

Remark 1.10 The prescription (|1.38|) generally breaks the complete positivity and even 
positivity of the reduced dynamics [HS] 

The convergence condition of the weak coupling limit depends on the time decay of 
the correlations of the bath. If we suppose Hi = Y,j V 3 S (g> V 3 B this condit ion reads 



dt\h jk (t)\(l + t) s < oo, V5 > 0, (1.41) 



o 



where 

h jk (t) = Tt b {e iHBt V 3 B e- iHBt V^ PB ) . (1.42) 

The condition Q1.41JI can be fulfilled only if the reservoir B is infinite [331134] . 

The Markovian approximations of the GME usually does not preserve the positivity 
of the states over which they act |35j . The following theorem defines Markovian approx- 
imations that do not suffer from lack of positivity. Noticeably, they turn out to be not 
only positive but also completely positive. 

Theorem 1.3 (Davies [36J Suppose the condition (|1.41j) is valid and the spectrum of H$ 
is discrete, then 



lim sup ||Af[p s (0)]-e^ +£ ^) f [p5(0)]||i=0, (1.43) 



where || • ||i is the norm introduced in Definition 1 1 . 2 1 and 



K\p s ]=\mi^-f dte- L ^Ke L ^[ps] (1.44) 



1.5. QUANTUM DYNAMICAL SEMIGROUPS 



21 



with 

/•oo 

K[p s }= / dt' e- Ls1/ Tr B (L I e < - Ls+L ^ t ' Li\p B ])\ps}. (1.45) 
J o 

The Markovian evolution defined by (jl,43j) is completely positive. 

In ljl.43[) we have to consider the averaged evolutor rather than K since the free 
motion and the dissipation have different time scales. 

Another kind of Markovian approximation is the so called Singular Coupling Limit (or 
Singular Reservoir Limit), introduced in m this case we do not rescale the Hamiltonian 
Ht but the two-time correlation functions of the environment: hjk(t) — > h E - k {t), with 
h e - k — ► CjkS(t) as e — ► 0. It has been shown in jHHl that such a case is equivalent to a weak 
coupling limit with a rescaled Hamiltonian in the form 

H £ T = e 2 H s + H B + eH I) (1.46) 

where the free Hamiltonian H$ is scaled by a factor e 2 and then the uncoupled and 
dissipative evolution are of the same order. 

Remark 1.11 The singular coupling limit is particularly relevant in the study of Marko- 
vian quantum open systems. Indeed, it is possible to show that any family of completely 
positive maps satisfying the semigroup composition law, describing the dynamics of a 
n- level system, can be derived by a singular coupling |39| . 

Remark 1.12 Another relevant property of the singular coupling is that, in this limit, 
the correlation functions of order higher than 2 vanish, that is the generalized master 
equation reduces to its Born approximation, and it is Markovian |4Uj . 

Historically the Markovian approximations have been studied considering two typical 
times, ts and t b , characterizing the system and the bath respectively; the density matrix 
ps can be extracted from the integral kernel JC(t) under the hypothesis t$/t b — > oo. There 
are two possibilities to realize this condition: either t$ — ► oo and t b remains constant, or 
t b — > while Ts remains constant [31]. The two alternatives correspond to weak coupling 
and singular coupling respectively. 

1.5 Quantum Dynamical Semigroups 

In the previous section we have introduced the set of dynamical maps we are interested 
in, namely, the set of Markovian time-evolutions describing the irreversible evolution of an 
open system under the influence of an external environment. It turns out that the structure 
of the generators of these dynamical maps is fully characterized by two theorems that we 
shall present after introducing the notion of quantum dynamical semigroups. 
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Definition 1.11 A set of time-evolutions {At,t 0} forms a quantum dynamical semi- 
group iff the following properties are fulfilled: 

l.i) Aj is completely positive Vt ^ 0; 

2.i) A t is trace-preserving: Tr(A t [p]) = Tr(p) Vi ^ 0, Vp G 

3.i) A t satisfies the forward in time composition law: A tl +t 2 [p] = ° At 2 [p] ^ti,ti ^ 0, 



4.i) the At are strongly continuous: limt_*o || A^[/o] — p\\\ = Vp G V(Sj). 

In the Heisenberg picture, defined by (|1.4|) . the corresponding properties of the dual map 
A* are: 

l.ii) A^ is completely positive Vt ^ 0; 

2.ii) A^ I = I, where I is the identity on B(Sj); 

3.ii) A* satisfies the forward in time composition law: A| 1+t [A] = A| x oA* 2 [A] Vii, ti ^ 0, 
G B(Sj); 

4.ii) A^ ultraweakly continuous: lim t _ Tr[p(A*L4] - A)] = Vp G P(£), V>1 G 
5.ii) A* is normal. 

Remark 1.13 The properties 3.i) and 4.i) imply that At = e Lt , where L is the gener- 
ator of the semigroup (and p(t) = L[p(t)]); then the maps {At,t ^ 0} are Markovian. 
Analogously, A* = e L *' (i.e. A(t) = L*[A(t)]) and L* is the generator of the Markovian 
semigroup {A*,t ^ 0}. 

Remark 1.14 The property 5.ii) guarantees that A* is the dual map of At |33j . 

The form of the generator of a quantum dynamical semigroup, if it is bounded, as a 
linear operator on the space of states V(TC), is completely determined as follows. 

Theorem 1.4 (Lindblad [U) Given a generic separable Hilbert space $), a bounded linear 
transformation L : V(Sj) — > V(f)) is the generator of a quantum dynamical semigroup on 



where H = and the Vj G B(H) satisfy £V VJVj G B(S)). The dual generator L*, in the 
Heisenberg picture, then reads 



vp g vm 



r{sj) iff 




(1.47) 




(1.48) 



1.5. QUANTUM DYNAMICAL SEMIGROUPS 



23 



Independently, Gorini et al. arrived at (|1.47|) in the case of finite dimensional system S. 

Theorem 1.5 (Gorini, Kossakowski, Sudarshan |42j ) Given a finite n-dimensional Hilbert 
space, S) n , a bounded operator L : V($Jn) - > V($3n) is the generator of a quantum dynam- 
ical semigroup on V(Sj n ) iff 

L[p] = -i[H, p] + ^ E ^ F ft + l F k>P F l]) - (1-49) 

k,l=l 

with H = W , {Fk, k = 0, . . . n 2 — 1} a basis in M n (C) (the space of n x n matrices with 
complex entries), Fq = l n /^/n, Tr(i^) = VA; ^ 0, Tr(i^'i^) = 5jj and C = [c^i] a positive 
matrix (whose entries depend on the choice of the basis {F^}). The Heisenberg generator 
L* : B(Sj n ) — > B{9) n ) is written as 

n 2 -l 

= + c fe j (f/L4, F fc ] + [F/, AjFfc) , [c M ] > 0. (1.50) 

fc,i=i 

Corollary 1.3 Given a n-dimensional Hilbert space, fj n , with n finite, a bounded opera- 
tor L : V($) n ) — ► V(9)n) is the generator of a trace and Hermiticity-preserving semigroup 
of maps iff it has the form in (jl.49[) . with H = H' the effective Hamiltonian and C = [cki] 
a Hermitian matrix. 

Remark 1.15 The form of the generator of a generic semigroup does not change if we 
relax the request of complete positivity, but the positivity of C is lost. 

The entries of the matrix C are not all independent because of the Hermiticity (or 
positivity) condition. A relation between these entries we shall use in the following is 
expressed by Lemma ll. II 

Lemma 1.1 The entries of the matrix C in the generator L of a trace and Hermiticity- 
preserving semigroup of maps acting on V(Sj), eq. (jl,49j) . satisfy the constraints 

|cw| ^ -(cjfcft + cu), k,l = 1, . . .n 2 - 1 (1.51) 

Proof: Expressing the operators in (|1,47|) as Vj = Y^=q v k^ k we nave c kl = J2j v k v i*- 
Then, applying Schwartz inequality, 

|Cfcj| ^ ^CkkQl < \{c k k + cu). 
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Remark 1.16 By rearranging equations (|1.47l) and (|1.49|) . the generator L can also be 
written as 

L[p] = -i[H, p]+Y^ {VjpV} - ±{p, v}Vj}\ (1.52a) 

3 

n 2 -l / 1 \ 

L[p] = -i[H, p]+Y, cm ( F kP F} - -{p, F}F k } ) ; (1.52b) 
k,l=i ^ ' 

all these expressions for L are generically called Lindblad forms. 

Remark 1.17 The Hamiltonian H is not the free Hamiltonian H$ of the system of in- 
terest; indeed, in general it contains dissipative terms produced by the interaction with 
the environment: contributions of first order in this interaction, the mean field term 
Ttb(Lj[pb]) in (|1.3Ujl . and of second order, as we shall see explicitely in Chapter |IJ 



Let us now analyze separately the contributions that appear in (|1.52a|) . 

The commutator with the Hamiltonian H generates a time-evolution trace and coher- 
ence preserving. The anticommutator between p and the positive operator R = ^ ■ VjVj 
may be interpreted as a friction term; indeed, let us consider the time-evolution produced 
only by the anticommutator and the commutator: 



p{t) = -i[H, p] - -{R, p} = -iKp{t) + ip(t)K\ 

£ (1.5o) 

p(t) = e- iKt p(0)e iKh , 

where K = H — \R. This evolution has the Weisskopf-Wigner form and it is typical of an 
unstable system, as previously remarked; in fact it causes a loss of probability, as we have 
seen in Section H~T1 

j t Tr[p(t)] = Tr[p(t)] = - Tr[p(t)R] ^ Vt > 0. 

Notice that, despite loss of probability, the time-evolution (|1.53|) preserves the coher- 
ence of states. In fact, in order to get decoherence, the term ^ ■ VjpVj in (|1.52aj) . called 
noise, is needed. Starting with a pure state, p = \ip)(ip\, 

±Tr[p(t) 2 ]\ t=0 = 2Tr[/j(0)p(0)] = 2^Tr[V>(0)V/ - V}V jP (0)] 

j (1.54) 

= -2Y,mv]v^) - (m\^^\Vj\^)) < o 

3 
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Then, in general p(t) 2 ^ p(t), so that pure states become mixtures. Moreover, the 
noise term balances the probability loss due to pure friction and the overall evolution is 
probability preserving: 



1.6 Factorized dynamics of quantum dynamical semigroups 

In this section we consider two systems of the same kind, Si and S2, in weak interaction 
with a same environment. We further suppose that the composite system Si + S2 evolves in 
time according to a reduced dynamics of factorized form At® At, where At is the dynamics 
of either systems S± or S2 and {A.t,t ^ 0} is assumed to form a Markovian semigroup so 
that the same is true for the set of maps {A t ® A t , t ^ 0}. We focus our attention on the 
positivity properties of these maps; in particular, in Theorem 11.61 we shall prove that, in 
the case of S± and S2 n-dimensional systems, the positivity of At ® At necessarily requires 
the complete positivity of At |43j . 

Given a semigroup {A t ,t ^ 0} generated by (|1.49|) . the justification why At should 
be completely positive and thus the matrix C = [cki] positive, is based on the fact that, 
otherwise, Af ® I n would fail to preserve the positivity of entangled states in V(Sj n ® $) n ) 
(see Proposition II. ljl . 

However, while the first factor in S) n ® Sj n refers to a concrete open quantum system 
evolving in time according to At, because of the possible interaction with a certain envi- 
ronment, the second factor represents a mere possibility of entanglement with anything 
described by a n-dimensional system and generically out of physical control. 

Instead, we argue that complete positivity is necessary to avoid physical inconsistencies 
in compound systems consisting of two n-dimensional systems that interact with a same 
environment, but not among themselves, neither directly, nor indirectly, that is through 
the environment itself. In such a case, the two systems are expected to evolve according to 
semigroups of linear maps {Af ® At,t ^ 0}, where At is the single open system dynamics 
obtained when only one of them is present in the environment. 

A necessary request for the physical consistency of such dynamics is that the At ® Af's 
preserve the positivity of all separable and entangled states of the compound system, which 
now describe physically concrete and controllable settings. 

Then for Markovian dynamics of open quantum systems the argument in favour of the 
necessity of complete positivity results strengthened with respect to the argument based 
on At ® I n (see the definition of complete positivity in Section I1.2JI . 

Theorem 1.6 If {Af, t ^ 0} is a semigroup of linear maps on V{9) n ) generated as in p.49|> . 
then the semigroup {At®A t ,t ^ 0} of linear maps over "P(fj n ®fj n ) is positivity-preserving 
if and only if {At , i ^ 0} consists of completely positive maps |43j . 



d_ 

dt 



Tr[p(t)] =Y,^[v 3 p{t)V] - \{V}Vj,p(t)}] =0, Vt > 0. 



(1.55) 
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The proof of Theorem II .61 is divided into several steps. We need just show the only if 
part; indeed, if At is completely positive, At <8> I n and I n <8> At are both positive and such 
is the composite map A t <8> A t = (A t (8) I n ) o (I n <8> At). Actually, it is completely positive 
(see e.g. (13]). 

Before starting with the proof of this theorem, some remarks are needed. 

Remark 1.18 If the Af's preserve the positivity of states of Sj n , At (8> At preserves the 
positivity of separable states of Sj n (8 fin'- this follows by the definition of separability 
(Def.EHJ). 

Remark 1.19 For generic positive linear maps At on the states of fj n , it does not follow 
that, if At (g) At is positivity-preserving, then At is completely positive. A counter example 
is the transposition r over M2(C), presented in E,emark ll.7l t(8>t is positivity-preserving, 
but r is not completely positive. We notice, however, that r cannot be among the At of 
a continuous semigroup over the states of M^{G) since it is not connected to the identity 
map. 

Remark 1.20 Actually, there are experimental situations that are describable by semi- 
groups {At (g> Af,t ^ 0}. For instance, as we shall see in Chapter neutral mesons may 
be imagined to suffer from dissipative effects due to a noisy background determined by 
Planck's scale physics. As decay products of spin 1 resonances, these mesons are produced 
in maximally entangled states and, while independently flying apart back to back, they ar- 
guably evolve according to semigroups {A t ®A t ,t ^ 0}. In such a context, whether At® At 
is positivity-preserving is crucial for concrete physical consistency [HH151Hd1I47U481I1^] . 

The proof of Theorem 11.61 proceeds as follows. 

Lemma 1.2 If {A t ,t ^ 0} is a semigroup generated by (|1.49j) . the set {At (g) A t ,t ^ 0} 
consists of positivity-preserving maps only if 

C M = (<P\(L®I n + I n ®L)U){^\M ^0 (1.56) 

for all orthogonal vector states \cf>), \ip) in where L is the generator on the right hand 
side of (jl.49f) and I„ is the identity operation on M n (C). 
Proof: The request of positivity preservation demands that 

G^(t) = {<f>\ (At® At) U) m \<f>) >0, (1.57) 

for all and G fi n <giS) n . Choosing (<f>\ip) = 0, if d</^(t)/di|t=o < 0, then Q^(t) ^ 
is violated in a neighborhood of t = 0. Thus (|1.56|) follows. 

Lemma 1.3 In the hypothesis of Lemma 11.21 let {\j),j = l,n} be an orthonormal basis 
of f) n , and <I>, \& the n x n matrices <1> = [<pij], ^ = [tpij] consisting of the coefficients of 
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the expansion of \<fi) and 1^) with respect to the basis {\j) (g) \k);j,k = l,n} of 9) n ® 9) r 
Then 



n 2 -l 



C M = c ab [Tr(V&F a ) Tr($*t F t) + Tr(($^) T F a ) Tr((^ t $) T F fe t )J , (1.58) 
0,6=1 

where C = [c a b] is the matrix of coefficients and F a ,Fi, the traceless matrices appearing 
in (|1.49|) . while X T denotes transposition of X with respect to the chosen basis. 
Proof: Let \(f>) = YTj,k=i Vjk\ j) ® \fy, I V 7 ) = Y%,k=i ^jk\j) ® then, one calculates 

C^ =^2^2^2(V*j<Pkli>pj'>Prl + Vjimi>jpi>ir)( i \ L [\p)( r \}\ k ) 

13 kl W (1.59) 

^^[rt i (wt )ir+( $t$ ) . ( $t $W(i | £[b)(r |p ) . 

The commutator and the anticommutator in the generator L drop from this equation 
since, given any K G M n (C), |p)(r|)|&) = Ki p 5 r k- In ()1.59j) . we can further sum over 
either r = k or i = p; in either = 0, we find Tr^ 1 " = (Tr^" 1 ")* = and 

the result follows. 

Lemma 1.4 In the hypothesis of Lemma |l.2( the matrix C = [c a b] m (|l-58[) must be 
positive definite. 

Proof: With any w = {w a , a = 1, n 2 — 1} G C" 2 " 1 , we consider W = \ *YTa=\ w *a^a, which 
is a traceless n x n matrix. If matrices \l/ and $ exist such that = W and $ = W T , 
then, from Lemmas II .21 and 11.31 and from the orthogonality of the matrices F a it follows 

n 2 -l 

Cftf = ^ c ab w* a w b ^ , (1.60) 

a, 6=1 

whence the positivity of C = [c ab ] and the proof of Theorem 11.61 Any matrix W and 
its transposed with respect to the given basis, W T , have the same elementary divisors; 
therefore, they are similar to the same canonical Jordan form and thus similar to each 
other [HOI. Let $ such that Q^W® = W T , that is we take as vector \<f>) G fj n <S> fj n 
the one whose components faj are the elements of the similarity matrix transforming 
the given W into its transposed W T . It then follows that ^ = and moreover 

\J/tcj> = Q^WQ = W T , which is what we need. 



1.7 The stochastic approach 

The description of an open quantum system can be performed employing an approach 
different from the projector technique presented in Section (jl.3j) . that is using a stochastic 
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equation |125j . The point of view is similar to that developed by Langevin to describe 
the Brownian motion: we introduce a time dependent term in the system Hamiltonian: 
Hg(t) = Hq + Hj(t), where Hj(t) characterizes the interaction between S and reservoir. 
The statistical operator that describes S is obtained by an averaging procedure over the 
noise: ps(t) = (Eg(t)), and £s(i) is the stochastic density matrix of S: 

± s (t) = L s (t)[Z s (t)] = (L + Lj(t))[E 5 (t)] (1-61) 

with Liouville operators Lo[Ss(t)] = ~ ^[Ho,^s(t)] and Li(t)[T,g(t)] = — i[Hi(t), Y*g(t)]. 
In the interaction picture ljl,61[) reads 

±W) = L\m s (t)}, (1.62) 

where 

S l s (t) = e itHo ^s{t)e- itHo , L\{t) = e- Lot L 7 (t)e iot . (1.63) 
The solution of (fl~62l is 

S l 5 (i) = TefiLWWYkip), (1.64) 
where T means time-ordering. Since p s (t) = (Ti s (t)), it follows that 

oo 
k=0 

with 

M k (t) = f ds x ds 2 ... H' 1 ds k {L\( Sl )4(s 2 ) . . . L}(s k )). (1.66) 
io jo io 

We assume that the initial state £s(0) is decoupled from the noisy environment so 
that p s (0) = (£s(0)) = S 5 (0) (and then p J 5 (0) = S J S (0)). Eq. (tO^j) is a series expansion 
in the interaction, containing the correlation functions of the fields defining Hj(t). If we 
suppose these stochastic fields to be Gaussian, only the even terms in (jl.65|) survive. 

Setting M(t) = ^Z^Mkit) = Xwc-^2fc(i)> we rewrite the differential form of equa- 
tion (|1.65j) as follows: 



p s (t) = M(t)M-\t)[f? s {t)\ = (M 2 (i) + M 4 (i) - M 2 (t)M 2 (t) + . . . J [p^(t)] (1.67) 
using the formal expressions 



M(t) = I + M 2 (t) + M 4 (t) + ..., M~ 1 (t) = I - M 2 (t) + [M 2 (t)] 2 - M 4 (t) + . . . . (1.68) 
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Assuming the system S to be a n-level system, it is possible to express the inter- 
action Hamiltonian as Hj(t) = Ylk=o Vk{t)F k , where F k is the basis of M n (C) in- 
troduced in (|1.49l) and V k (t) are Gaussian stochastic fields with correlation functions 
hki(t — s) = (Vk(t)Vi(s)) . Then, up to the second order in the interaction, eq. (jl.67() 
assumes the form 

™ 2 -i ft 

PW) = " E / dsh M (s)[F k (t), [F t (t - s),pUt% (1-69) 
k,l=i Jo 

with F k (t) = e Lot [F k ] = e itH ° F k e~ itHo . 

Turning back to the Schrodinger picture, ps(t) = e~ ttH ° p l s (t)e ltH ° and we find: 

n 2 -l 

p s (t) = -i[H ,p s (t)} + d kj {t)[F k , [Fj,p s (t)]], (1.70) 

k,j=i 

with 

n 2 -l ,t 

dkj{t) = ~ T2 I h kl {s)ui j {-s)ds (1.71) 
i=i Jo 

and uij(-s) defined by F t (-s) = Y^Tq 1 u ij(~s)Fj. 

The Markovian approximation of the reduced dynamics (|1.70j) is performed with the 
hypothesis that the correlation functions h k i{t) decay so fast that we can substitute +oo for 
the parameter t in (|1,71|) : in this way the coefficients d k j(t) become time independent. This 
assumption, as we have already seen, is the same as considering the system evolution for 
times long enough, so that the transient terms (the memory contributions, non Markovian) 
disappear. 

The reduced dynamics so obtained, although Markovian, in general is not completely 
positive nor positive (in fact the generator of the dissipative part in general has a different 
form than that expressed in Theorems 11.41 and II. 5|) : the stochastic approach suffers from 
this failure. Equations of motion in the form 1)1.70(1 . in the two-dimensional case, after the 
Markovian limit, are called Redfield-Bloch equations )27U125lll26j : in Chapter we shall 
explicitly show instances where lack of positivity occurs. 

It has been argued that it is possible to consistently use these time-evolutions, even if 
they are not positive, provided they act only on a suitable subset of "P(-fj n ), constituted by 
statistical operators whose positivity is preserved. Such a restriction is usually explained 
in terms of a slippage of initial conditions |3QII31j : the initial transient, neglected in 
the Markov limit, is supposed to select from V($j n ) exactly the set of states which remain 
positive under the dynamical map considered. In such a perspective, the long-time Markov 
approximation of the dynamics, possibly not positive, should be coupled to a redefinition 
of initial conditions, generated by the action, on "P(i^ n ), of the short-time non-Markovian 
dynamics. In such a way the positivity of ps(t) should be saved. However, it can be shown 
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that the prescription needed to guarantee the consistency of the action of a Markovian 
dynamical map At in general does not hold when two systems are immersed in the same 
environment and evolve with the factorized dynamics At <8) At, starting from an entangled 
state: a new prescription has to be enforced [HJ- This argument can be repeated when 
three or more systems, living in the same environment, are taken into account: there does 
not exist a general prescription, valid for any number of systems interacting with the same 
environment. In this sense, complete positivity appears at least as a cheaper tool to avoid 
physical inconsistencies. 

1.8 Entropy and long time behavior in open quantum sys- 
tems 

Typically, an open system exchanges energy and entropy (also matter in general) with 
its environment. Because of the vastly different number of their degrees of freedom, one 
expects the noise coming from the environment to always increase the randomness of the 
system and thus its entropy. 

For later use, in this section we shortly review some results about the long time behavior 
of a state evolving according to quantum dynamical semigroups, and of its entropy. 

Definition 1.12 Given a statistical operator p G T > (Sj), its von Neumann entropy S(p) is 
defined by S(p) = — Tr(plogp). 

For a review about entropy and its properties, see pHI; here it is sufficient to remark that 
the von Neumann entropy is always non-negative, S(p) ^ 0, Vp G ^(Sj), and S(p) = if 
and only p is a pure state. Furthermore, if S) is finite dimensional, say dim(fj) = n, the 
entropy is greatest for the maximally mixed state p n = I n /n: S(p n ) = logn. 

The time-evolution of entropy is related to the long time behavior of states evolving 
under quantum dynamical semigroups, to the existence of invariant states, and to the 
approach to equilibrium of the considered open systems. 

Definition 1.13 A state po E V(7i) is said invariant (or stationary ) under the action 
of a quantum dynamical semigroup {At, t ^ 0} iff A t po = Po, Vi ^ or, equivalently, 
L[po] = 0, where L is the generator of the semigroup, expressed in (|1.47|) . 

Definition 1.14 The dynamical semigroup {At,t ^ 0} is called relaxing iff there exists 
an invariant state po £ V(7i) such that \/p G V(H), limt_>oo A t p = po- 

A complete classification of all types of stationary states associated to relaxing semi- 
groups is given in [221; see a ^ so |53U54j . For later reference we are only interested in the 
following considerations. 

As we have seen in Section H. 51 the dual time-evolution A^ always preserves the identity 
I, and then its generator, denoted by L* , maps I into the null operator: L*\J] = 0. In 
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general, the same does not hold for the evolution map At and L[I] 7^ 0. However, as we 
shall shortly see, there is a tight connection between the entropy time-evolution of the 
system and the invariance of I, technically called unitality. 

Definition 1.15 A quantum dynamical semigroup {At,t ^ 0} is called unital iff the 
identity operator I 6 V(7i), is invariant under its action. 

Remark 1.21 If Aj is unital, the maximally mixed state p n is invariant. 

Remark 1.22 If Af is unital, the operators A/, defined in Corollary 11.11 must satisfy 

Remark 1.23 If A^ is unital, the operators Vk defined in Theorem 11.41 must satisfy 

Now, it turns out that, for an n- level system, the unitality of the dynamics is a necessary 
condition for the non-decreasing behavior of its entropy. Indeed, if unitality is not fulfilled, 
the maximally mixed state p n is not invariant under the reduced dynamics considered: 
p n = L[p n ] 7^ 0; therefore p n {t) = A t [/? n ] 7^ p n and then, necessarily, the entropy decreases: 

S[ Pn (t)} < S{ Pn ). 

In the next chapter we shall prove that, for n = 2, the unitality of quantum dynamical 
semigroups is also a sufficient condition in order to have a non-decreasing entropy. 
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Chapter 2 

Bidimensional Systems 



In view of later applications, in this chapter we restrict our attention to quantum open 
systems described by two-dimensional Hilbert spaces. The parameters describing gen- 
eral Markovian dynamics are introduced; positivity and complete positivity of the time- 
evolution are expressed by sets of inequalities among them. 

Although the equations of motion can be integrated exactly, in many cases it is prefer- 
able to express their solutions as perturbative expansions in the dissipative parameters, 
provided they are small with respect to the parameters describing the dynamics of the 
system alone, in absence of the environment. This computation, to leading order in the 
dissipative parameters, is performed in Section |2~H 

2.1 Evolution of states 

The physical systems we shall consider in the next chapters will be described by means 
of two-dimensional Hilbert spaces. Therefore, their states will be represented by density 
matrices corresponding to positive elements of M2(C) with unit trace. We will use two 
different basis for the algebra M2(C), namely: 

• the set {Pi,i = 1, . . . 4}, where 

--Go). *=(S!). Ms;). «=(;s). <-) 

with P ] l 2 = P h2 and Pl A = P 4 ,3; 

• the set {<7j, i = 0, . . . 3}, that is the set of Pauli matrices plus the identity: 

CT ° =i2= i)' cti= (i o)' a2= o 1 )' ct3= -0' 

(2.2) 
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satisfying crj = crj, i = 0, . . . 3; Tr(<7j) = 0, i = 4, . . . 3; Tr(crjcrj) = 25^, i, j = 0, . . . 3. 
Accordingly, states p will be written, with respect to these two basis, as 

4 

P = Yl P iPi 
i=l 
3 

i=0 

We use upper and lower indices to distinguish the coefficients in the two different 
expansions. We have p. L £ C and p4 = pj, while p l 6 R. 

We are interested in Markovian dynamics for the states p; they are fully characterized 
by Theorem 11.51 and Corollary 11.31 The basis {Fi = a{/y/2,i = 0, ... 3} fulfills the prop- 
erties required in Theorem IT31 F = h/y/2, Tr(i^) = Vi + 0, Tv(F^ Fj) = <%; then, the 
more general form for the generator of a Markovian dynamics is 

L[p] = (L H + L D )[p], (2.4) 

where Ljj[-] and Ld[-] are the Hamiltonian, respectively dissipative part, of the generator. 
Lh , for a stable system, is proportional to the commutator with a Hamiltonian H = H< 
that can be chosen traceless without loss of generality, then H = Yli=i hi&i, hi £ R; 
can be written in the Kossakowski form: 

Lh[p] = ~i[H,p], L D [p] = i °kl (vkPVi - ^{p,<ria k }j (2-5) 

with C = [cm] a Hermitian matrix. 

For unstable systems H is non-Hermitian (i.e. hi G C) of Weisskopf-Wigner type (see 
the discussion in Section ll.ljl and the time-evolution does not preserve the trace of the 
states p; the generator has thus the form ()2.4j) with a modified Hamiltonian part: 

L H [p] = -iHp + ipH^ (2.6) 

and will be used when dealing with neutral mesons systems; at present we restrict our 
attention to trace preserving dynamics with Lh in form (|2.5|) . 

It is convenient to introduce a vector notation and write the density matrix p as a 
4-component vector \p) = (p° , p 1 , p 2 , p 3 ) whose entries are the coefficients of the expansion 
of p in the basis (|2.2|) : in fact, since equation (|2.4j) is linear in p, it can be rewritten in 
vectorial form as a kind of Schrddinger-like equation: 



f Pi Pi\ 

V At P2 J 

(2.3) 

/ p° + p 3 p 1 - ip 2 \ 

V p 1 + ip 2 p° - p 3 )' 



-\p(t)) = C\p(t)) = (H + V)\p(t)}, 



(2.7) 
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where 7i and T> are 4x4 real matrices, representing the Hamiltonian, respectively the 
dissipative part of ([2.4)1 . These matrices have the form 



n 



( 











\ 








^3 


-h 2 







-h 3 





hi 




I o 


h 2 


-h x 








V 



( \ 

u a b c 



v 
\ w 



a (3 
P 7 / 



(2.8) 



In ([2.8)1 we have introduced nine real coefficients to express the dissipative part T>. 
They depend on the entries of the matrix C in ([2.5)1 : explicitly [T]l53]: 



/ o 

2 Im c 32 
2 Im ci 3 
\ 2Imc 2 i 





-C22 - C33 

Re C12 
Reci 3 





Re C12 

-Cll - C33 

Rec 23 





Re ci 3 
Re c 23 
-en - C22 / 



(2.9) 



conversely the matrix C can be written in terms of these dissipative parameters: 



C = 




—2b + iw —2c — iv 
a + 7 — a —2/3 + iu 
—2(5 — iu a + a — 7 



(2.10) 



Remark 2.1 The Hamiltonian contribution W is skew-symmetric since it represent a 
commutator. 



Remark 2.2 The entries of the first row of both 7i and D vanish since we have required 
the dynamics to be trace preserving. Indeed, from ()2.8j) it follows that p°(t) = 0, as it 
should be since Tr[p(t)] = 2p°(t) and then dTr[p(t)]/dt = 0. 



Using the basis ([2.1)1 p gets represented by the 4-vector 
evolves in time according to 



d_ 

dt 



\p(t)) = c\p(t)) = (n + f>)\p(t)), 



= (pi,P2,P3,Pa), which 



(2.11) 



where Ti. and T> are respectively the Hamiltonian and the dissipative contribution in this 
representation 1 . From the relations in ([2.3)1 it follows that the two vectors representing 
the state p are related by 



\p(t)) = V\p(t)), V 



( 1 

10 

1 -i 

\ 1 i 



1 \ 

-1 



/ 



(2.12) 



J We use hats to characterize states as well as operators in this representation 



36 



2. BIDIMENSIONAL SYSTEMS 



Thus, 



n = vnv- 1 



and 



with 





h 2 + ihi 
\ h 2 — %h\ 



V = VVV 



-i 



— h 2 + ih\ 

h 2 — ih\ 

-h 2 — ih\ —2ihs 

-h 2 + ih\ 



—7 — w 
7 + w 
-C* - U* 

-c-u 



7 — w 
w — 7 
C* - U* 
C-U 



-h 2 - ih\ \ 
h 2 + ihi 


2ih 3 J 



-C -C* \ 

c c* 

-A B 

B* -A ) 



(2.13) 



(2.14) 



A = a + a, B = a — a + 2ib, C = c + if3, U = u + iv. 



(2.15) 



Before discussing solutions of the equations of motion we discuss the positivity of the 
dynamics. 



2.2 Positivity and Complete Positivity 

The requests of positivity or complete positivity for the time-evolutions described in the 
previous section are expressed by different constraints on the dissipative parameters ap- 
pearing in D. 

Remark 2.3 The positivity of A G M 2 (C) is equivalent to Tr(A) ^ and Bet(A) ^ 0, 
since we are in two dimensions. Since states p have Tr(p) = 1, so that p° = 1/2, they 
are positive if and only if Det(p) = 1/4 — X^=i(/ "') 2 ^ 0. A state p is pure (i.e. a one- 
dimensional projection) if and only if Det(p) = 0. A mixed state has Y^=i(p^) 2 < 1/4. 

Theorem 2.1 If the trace and Hermiticity preserving time-evolution (|2.7j) with 7i and T> 
in (|2.8f) is positive, then the real entries of T> satisfy the inequalities: 

a ^ ( aa^b 2 

a ^ ; < a-f^c 2 ; aaj + 2bc[3 - a(3 2 - ac 2 - jb 2 > 0. (2.16) 
7^0 I 07 ^ (5 2 

Proof: p(t) ^ if and only if Det[p(i)] ^ (Remark 12. 3() . Choosing the initial state 
\p a {0)) = (1/2, p\p 2 ,p 3 ) to be pure, that is £? =1 (p>') 2 = V 4 , th en dBet[p a (t)}/dt\ t=0 ^ 
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must hold if p a (t) has to remain a state. Explicitly, 

1 3 
±Det[p a (t)]\ t= o = -2j2fj(t)P j (t)\t=o 



j,k=l j=l 

4 ^2 Ji\ 



Analogously, starting with the pure state \pb(0)) = (1/2) ~p\ ~ P i ~ P ) we find 

3 3 

DetfnJ^lL-n = -2 ' 

dt 



i Det[p b (t)]| i=0 = -2 + ^U^'(0)p fc (0) + £>;oP J '(0) ^ 0; (2.18) 

j,fe=l j=i 



adding equations (|2.17|) and 1)2.18(1 and considering that H is a skew symmetric matrix, 
we are left with 

3 

V jk p k (0)^(0) ^0. (2.19) 



Introducing the sub-matrix D as 

/ \ 



a & c 

V=\ba/3\, (2.20) 

c /3 7 



the condition (|2.19|) becomes T> ^ 0; which is equivalent to (|2.16() by considering the 
minors of T>. 

Theorem 2.2 The trace and Hermiticity preserving time-evolution (|2.7() is completely 
positive iff the entries of T> satisfy the inequalities: 

2^ = + 7-0^0 ( RS^b 2 + w 2 /A 
25 = + 7- 0^0 ; I RT ^ c 2 + v 2 /4 ; 
2T = a + a-7^0 ( ST^/3 2 + u 2 /4 

RST - 2bc(3 + ±6m; + i/3vw + |cmu - S (c 2 + ±u 2 ) - R {(5 2 + ±u 2 ) - T (b 2 + ±u; 2 ) ^ 0. 

(2.21) 

Proof: The thesis follows from the fact that the time-evolution is completely positive if 
and only if the matrix C in (|2.10() is positive fTheorem ll.5|) . 

Remark 2.4 Note that the constraints expressed in these two theorems are not all inde- 
pendent. Also, the constraints (|2.2H) imply those in (|2.16|) . Indeed, CP is stronger than 
positivity. 
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2.3 Positivity and CP for entropy increasing evolutions 

There are physical situations in which the reduced time-evolution of states does necessarily 
increase the von Neumann entropy (see Definition 11.12(1 and states becomes less ordered 
as the time passes. For instance, in Chapters El and 03 we shall describe irreversible time- 
evolutions that are entropy increasing. 

As anticipated in Section ll.8l any entropy non decreasing time-evolution can be fully 
characterized as follows. 

Theorem 2.3 The trace and Hermiticity preserving time-evolution (|2.7j) with TC and T> 
in (|2.8|) has a non decreasing von Neumann entropy iSu = v = w = and T> ^ 
(i.e. the time evolution (|2.7|) is unital), where u, v and w are the parameters defined 
in equations (|2.8|) and (|2.9|) and T> is the matrix in eq. (|2.2U|) . Moreover, the matrix C 
in (|2.10j) is symmetric |39U55j . 

Proof: The eigenvalues, respectively eigenvectors, of a generic initial state represented by 
the vector |p(0)> = (1/2, p 1 , p 2 , p 3 ) with \p\ 2 = Y^=i(p j ) 2 < V 4 , are 

p± = i(l±2H); \p±) = ^(\p\,±p\±p 2 ,±p 3 ). (2.22) 

Consider now 

J t S[p(t)]\t=o = -j t Tv[p(t)logp(t)}\ t=0 = -Tr[p(i)|t =0 logp(0)] (2.23) 

where we have made use of the trace preserving property of the time-evolution. Evaluating 
the trace over the states \p+) and \p~) we find 

j t S[p{t)]\ t=0 = -logp + Tr[p(t)| t=0 |p + )( /0+ |]-logp_Tr[p(i)| t= o|p-}(p-|] = 



1 

R 



3 3 

I^^(o)+^^y(o)^'(o) 



i=l i,j=l 



, P+ 

log — 

P- 



(2.24) 



having taken into account the skew-symmetry of 7i. 
Suppose u = v = w = (that is T>iq = Vi) and V ^ 0. Then 

j t S[p(t)]\ t=0 = p E ^W(0)log^ (2.25) 

and, as /0+ ^ p_ , the von Neumann entropy is not decreasing. 

Vice versa, suppose ^S[p(t)]\t = Q ^ and choose p 1 = e, p 2 = p 3 = 0, with |e| 1/2. 
Then 

d „ 2 (\ 



di S[p(t)]\ t=0 = -[-ue + as 2 ) log ^. (2.26) 
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Since \u\ ^ a (see Lemma 1 1.1 (I . the right hand side of (|2.26|) could be made negative by 
choosing sign(e) = — sign(tt) and |e| ^ \u\/2a, unless u = 0. The same argument implies 
v = and w = choosing p 2 = e, p 1 = p 3 = and respectively p 3 = e, p 1 = p 2 = 0. Then 
eq. ([2.25ft holds, whose positivity imply T>ij ^ 0. 

Theorem 2.4 The trace and Hermiticity preserving, entropy non decreasing, time evo- 
lution (|2.7j) is positive iff the real parameters of T> satisfy the inequalities (|2.16j) . 
Proof: The if part is a direct consequence of Theorem l2.11 Vice versa, considering a generic 
initial state |p(0)> = (1/2, p 1 , p 2 , p 3 ) with Det[p(0)] ^ 0, we have 

d 3 

- Bet[p(t)]\ t=0 = Ai/(0)^'(0) > (2.27) 

since T> ^ 0. Consequently Det[p(i)] ^ 0, Vi; thus the evolution is positive (Remark 12. 3j) . 
Indeed, suppose there is a time t for which Det[/?(t)] < 0. By continuity, there must be 
a time t' , with < t' < t such that Det[p(t')] = and dBet[p(t')]/dt < 0. But this can 
never happen, if T> ^ 0. 

Theorem 2.5 fl^ l4ini4l^ l4T ^ l4T |l4^l4lfll55] The trace and Hermiticity preserving, entropy 
non decreasing, time-evolution (|2.7|) is completely positive iff the real parameters of T> 
satisfy the inequalities 

2^ = + 7- 0^0 ( RS^b 2 

25 = + 7- 0^0 ; < RT^c 2 ; RST - 2bc(5 - Sc 2 - R(3 2 - Tb 2 ^ 0. (2.28) 
2T = a + a- 7^0 [ ST ^ [3 2 

Proof: The thesis follows from Theorem 12.21 putting u = v = w = 0. 

2.4 Perturbative expansion 

We want to integrate the equation of motion for p, i.e. we want to solve either (|2.7|) 
or (|2.11j) . In vectorial notation the solutions will be denoted by 

\p(t))=Q t \p(0)}, g t ^e ct = e^ n+ ^ t (2.29) 

and 

\p(t))=Gt\p(0)), Qt^e tt = e ( ^ + ^ t (2.30) 

respectively. Both Qt and Qt are 4x4 matrices and Qt = VQtV -1 , with V expressed 
in (l2~T2l . 

Although it is always possible, at least in principle, to compute the form of Qt and 
Qt without approximations, often these explicit expressions are rather involved and not 
particularly useful in order to capture the essence of the dissipative processes. 
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Usually, however, the parameters appearing in the dissipative part of either ()2.7|) 
or 1)2.1 1J) can be considered small with respect to those characterizing the Hamiltonian 
contribution. It is thus possible to express both Q t and Qt as series expansions in these 
parameters , 



Qt = Qt +Ql+Q 2 t+-.. (2.31) 



with 



g° = e nt , Q}= [ t ds 1 e n ^Ve Hs \ 
Jo 



Q 2 t = f d Sl r ds 2 e H{t - Sl) Ve Hisi - S2) Ve ns2 
Jo Jo 



(2.32) 



Equivalent expressions hold for Q%. 

Remark 2.5 Note that Q\ = VQ ] t V~ l for any j, that is the similarity relation holds at 
any order in the perturbative expansion. 

In this section we write the time-evolutor to leading order in the dissipative coeffi- 
cients. Indeed, in the coming chapters, we will encounter physical situations in which the 
dissipative parameters are small with respect to the parameters describing the free time- 
evolution (i.e. without environment) and the formalism here developed will be applied. 
We perform the computation in the representation ()2.11)) . 

We suppose H as general as possible, so that we allow for a non Hermitian part. Hence 
its eigenvalues ei and e 2 of H can be complex and its eigenstates non orthogonal: 

H|€ 1 > = ei|ei> ■ 

€i = rrii - -7i m,,^ eR, <ei|e 2 ) / 0; (2.33) 

H\e 2 ) = e 2 |e 2 ) 

rrii and 7$ are the energies, respectively the widths, of the eigenstates. 

Since the underlying Hilbert space is two-dimensional, in any given basis the two 
eigenvectors can be written as |ei) = (pi,qi) and \e 2 ) = (p 2 , —q 2 ), with the normalization 
conditions 

M 2 + M 2 = l, \p 2 \ 2 + \q 2 \ 2 = 1. (2.34) 

Remark 2.6 If H is Hermitian then (ei|e2) = and 7, = 0, that is ej G R. The condition 
of orthogonality of the two eigenstates reads 



PiP 2 ~ QiQ 2 = 0. 



(2.35) 
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It is convenient to work in the basis that diagonalizes the matrix Ti,, hence T> will be 
transformed in a more complex matrix 2 ; the whole computation will result simpler |49j . 
Accordingly we choose as basis for our bidimensional Hilbert space 9)2 the pair of eigen- 
states of the Hamiltonian, given in 1)2. 33 J) . H is diagonalized by means of an invertible 
transformation U (unitary if H is Hermitian); its diagonal form will be denoted by an 
upper index d, as any operator represented in the basis ()2.33)) : 

H d = U^HU; H d =( € ' V U = ( Pl P2 ) . (2.36) 
It is useful to introduce the two ratios 

ri = — , r 2 = — ; (2.37) 
Qi Q2 

this allows to write U = RQ, with 

*=(? -i)' ^(o »)• (238 > 

so that the similarity transformation connecting H and H d depends solely on the param- 
eters (J2~3T1) : 

H d = BT l HR. (2.39) 

Remark 2.7 Note that R is not unitary even if H is Hermitian. In this case U is unitary 
and then r\r\ = 1. 

The evolution of a state p, expressed in 1)2.4)1 with Hamiltonian and dissipative contri- 
butions as in 1)2.5)1 or 1)2.6)1 . in the new basis takes the form 

p d (t) = (L d H + L d D )[p d (t)], (2.40) 

where 

p d (t) = R^p^R-^, L d H>D [] = R- l L H , D [R ■ flt]fl-it. (2.41) 
Using the vector representation of p d we have 

j t \p d {t)) =C d \p(t)) = (H d + V d )\p d (t)), (2.42) 

where \p d (t)} is the vector whose components are the coefficients of the expansion of p d (t) 
in the basis 1)2.1)1 . It is obtained from \p(t)) by the linear transformation 

\p d (t))=Tl\p(t)), (2-43) 



2 Notice that the choice of the representation of the state and the choice of the basis of the Hilbert space 
are completely independent 
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whose explicit expressions can be obtained from (|2.41|) ; 

2 



K 



1 



/ 1 
1 

1 

V i 



l r 2| 

N 2 

-r 2 r\ 



* 


r 2 \ 


— r* 


-n 


* 


r2 




-n / 



(2.44) 
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-71 











\ 







-72 


















-r_ 







V 











-r + 


/ 



(2.45) 



v d = nvn- 1 



In + r2\ 

The 4x4 matrices in (|2.42|) can now be explicitly obtained: 

n d = nnw 

with T± = T ± iAm, T = (71 + 72V2 and Am = m-2 — mi; and 

/ A + Ai £ + <r A + <5i A* + 5* \ 

H + e A + A 2 $ + 0i $* + (t>t 

4> 2 - «52 - A* + w G + 6 

\ fa - $ flj - A Q* + 0* fi + a;* / 

where small and capital Greek letters refer respectively to contributions related to the 
entropy behavior (i.e. functions of u, v and w in l|2.8|) ) and remaining parameters (i.e. 
those that are functions of a, b, c, a, (5 and 7 in (|2.8jl ). Their expressions are listed in 
Appendix both for stable and unstable systems. 

Remark 2.8 For later convenience we observe that Ai, A2, cr, £ G R and 6\, 82, (j>i, $2, 
uj, 9 £ C; on the other hand A, S, S, !] e R and A, <£, G G C. 

In the basis of the eigenvectors of H and in the chosen vectorial representation, the 
evolution operator is denoted by Qf. Up to first order contributions in the dissipative 
parameters, Qf 



(2.46) 



gdO + gdl. from (j^J) we find 



i ( /0 

h 





-72* 








-r_t 



\ 



e 

\ e 

The next order contribution is obtained, accordingly with (|2.32|) . as 



-r+t 



(2.47) 



St 
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The operator Qf so obtained contains terms which are linear in t; they disappear if 
we reconstruct the exponential dependences out of them. For example, summing the 
contributions from Qf° and Gt, the entry in first column and row of Qf for an unstable 
system reads 

0f) u = [1 + (A + A^iK^ « e (A+Ai)* e -7it. 

by looking at this expression we observe that, at the first order in the dissipative contri- 
butions, what we physically mean by width of the state |ei) is not the 71 appearing in the 
Hamiltonian H, but rather 71 — A — Aj . Accordingly, we redefine 71 so that 

(Qf) n = e -(7i-A-A 1 )t ^ e - 7 it_ 

In the same way also 72 and Ar± = Ar ± 2iAm get redefined. The final, explicit 
expressions of the components of Qf, both for stable or unstable systems, are collected in 
Appendix iBl 

Remark 2.9 In general, the n-th term in the expansion (|2.31|) contains polynomials in 
t of order n; these terms represent truncated exponential behaviors and it is possible to 
consistently re-absorb them at any order in the perturbative expansion by reconstructing 
the exponentials they represent. A truncated perturbative expansions does not hold for a 
too long time since the neglected terms become more and more relevant as the time passes, 
then the reabsorbtion of the contributions polynomial in t could seem useless rather then 
physically justified. However, this procedure allows to redefine the parameters rrii and 7, 
by dissipative terms; in such a way they coincide with physical quantities (see Chapter El 
for further details). 

Remark 2.10 The time-evolutors in the physical basis are obtained as Gt = !Z~ l GflZ 
and Qt = V~ 1 GtV m the two different representations adopted. However their explicit 
expressions are unnecessary for the later developments; indeed in the formalism of density 
matrices the mean value of any observable, represented by a Hermitian operator A, is 
obtained by a trace operation over the evolving state p(t); this operation is invariant 
under any change of basis in the underlying Hilbert space: 

(A) p = Tr[Ap(t)} = Tr[A d p d (t)}. (2.48) 

It turns out to be more convenient to represent the operator A in the basis of the eigenstates 
of H and to compute (A) p via (|2.48j) . rather than re-expressing p d (t) in the physical basis. 

Remark 2.11 In the case of a stable system, the reabsorption of the terms linear in time 
is more difficult because it involves also the off diagonal entries of Gf ■ The final expressions 
are reported in Appendix IBl 

These results will be used in the coming chapters for the description of the phenomenol- 
ogy of stable and unstable particles immersed in various specific environments. 
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Chapter 3 

Neutral B Mesons 



In this chapter we concentrate on the study of the neutral B mesons system treated as an 
open system. 

This system is presently under both theoretical and experimental investigation since it 
enables the study of some specific properties of the Standard Model, e.g. the topic of the 
violation of the CP symmetry [S71|5H] (charge conjugation plus space reversal), but also 
of the CPT symmetry (CP plus time reversal) and the flavors oscillations. In particular 
single B mesons or entangled pairs of B mesons physics are intensively studied at the 
colliders or at the so-called B mesons factories |87j . 

Usually, the neutral B mesons system is considered as a closed system. Instead, in 
the following, we assume this system to be immersed in a thermal bath of gravitational 
origin. We shall adopt a phenomenological approach, that does not make any hypothesis 
about the microscopic nature of the interaction between B mesons and environment; 
we suppose solely the coupling to be weak and the entropy non-decreasing. It turns 
out that the resulting irreversible dynamics for the B mesons can be parameterized by 
six phenomenological coefficients that affect in a very specific way the relevant physical 
quantities (decays rates and asymmetries) jl^l • hi principle, an experimental study of 
these quantities enable us to evaluate the magnitude of these dissipative parameters. 

As usual, we shall describe our system by a bidimensional Hilbert space, using the 
formalism developed in Chapter |21 and in particular the time evolution operator calculated 
in Section ITU 

This dissipative approach has been also applied to describe the neutral K mesons 
system [So1l67U68ll69U44lli5ll4l)ll47U48| . For another approach to decoherence in the neutral 
B mesons system, characterized by only one dissipative parameter, see |7()j . 

3.1 The physics of B mesons 

In the standard phenomenological approach the neutral B meson system can be described 
by means of a state or, more in general, a statistical operator, in a two-dimensional Hilbert 
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space f)2- The effective Hamiltonian H of this system is a non Hermitian operator of the 
Wigner-Weisskopf type [37 ] l58j. 

Hi Hi 



H = M--r; H=^ R4 £ ) (3.1) 
Several basis in $)% are of interest for the discussion that follows: 

• The set {|B°), |B°)} of the physical states (called also flavor states); 

• The set {|B^), |B#)} of the eigenstates of the Hamiltonian of the system: 

H\B L ) = e L \B L ), e L = m L -^ L 
H\B H ) = e H \B H ), e H = m H -^ H 



(3.2) 



where m and 7 denote respectively the masses and the widths of the two eigenstates 
and the subscripts L, H means "light" , "heavy" respectively, since it turns out that 
rriL <C ran [571135] . These states are not ortogonal: (Bl\Bh) 7^ 0; 

The set {|Bi), IB2)} of the eigenstates of CP, the operator performing the CP 
transformation (charge conjugation plus parity): 



CPlBi) = |Bi) 

, (3-3) 
CP|B 2 ) = -|B 2 ) V ; 

The eigenvalues of CP are 1 and —1 since (CP) 2 = do. The CP operator exchange 
the two physical states: 



CPIB ) = e^|B°) 

_n (3-4) 
CP|B°) = e-^|B°) 

and (f is a phase depending on the convention adopted in defining the physical states. 

Remark 3.1 In the treatement of kaon phenomenology the distinction between the two 
eigenstates of H is based on the large difference between the two widths and the terms 
"short" and "long" are used. 

The relations between the basis (|3.2[) and (|3.3j) are: 

\B L )=p L \B°)+q L \B^), \p L \ 2 + \q L \ 2 = l 
\B H ) = PH \B°) - ^|B U ), \p H \ 2 + \q H \ 2 = 1 
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and 



B 2 ) 



Bi> 




(3.6) 



The distinction between these three basis is one of the reasons of the physical relevance 
of the B mesons system. 

Since the physical states, produced in the interactions, are linear superpositions of 
the eigenstates of H, an initially pure state, e.g. |B°), would develope a non- vanishing 
component on the other state, e.g. |B°), as time passes. We have then transitions between 
the states |B°) and |B ), called flavor oscillations, as the probability of transition between 
the two flavors behaves like a periodic function. 

On the other side the distinction between eigenstates of CP and H is a consequence 
of CP violation: it is not possible to diagonalise at the same time both these operators. 
In the standard Model, CP violation is due to the presence, in the Cabibbo Kobayashi 
Maskawa matrix 1 , of a complex phase that cannot be reabsorbed by a redefinition of 
the fields describing the quarks [57ll58j : this violation is expected to be greater in the 
B mesons sector with respect to the K mesons case. CP invariance is equivalent to 
\ph\ = \qn\ = \pl\ = 

The adopted phenomenological approach can further account for violations of other 
symmetries that are preserved by the Standard Model, like CPT invariance and the phe- 
nomenological rule AB = AQ. In this relation, B is the beauty quantum number and Q 
the charge of the decaying quarks. 

Further, it is possible to describe time-evolutions different from the standard ones, tak- 
ing into account a possible interaction with an external environment. As discussed below, 
the modifications of the relevant quantities produced by the irreversible dynamics turn 
out to have a characteristic signature: the open system dynamics cannot be confused with 
the others possible symmetries violations. To be as general as possible in our discussion, 
we shall consider all of these different possibilities together. 
After the definition of the two ratios 




PH 



(3.7) 



the Hamiltonian (|3.1j) can be rewritten as 





1 This is the matrix giving the mixing in the hadronic sector of the Standard Model 
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For future reference we find useful to define three complex parameters, a, £ and 6 |57l 

EH]: 

H z H ~ 1 Q r L - r H Hi - H 2 

o = r H r L = — , £ = , 6 = ■ = ; (3.9 

Hi \a\ + 1 r L + r H e L - e H 

as we shall see in the coming sections, by means of them it is possible to parametrize in a 
convenient way violations of the discrete symmetries. 



3.2 The gravitational bath 

The physical motivation of an open quantum system treatement of the neutral B mesons 
system (or a generic neutral mesons system) comes from quantum gravity. A satisfactory 
and complete theory of quantum gravity has not been found yet. Nevertheless, general 
physical arguments involving quantum fluctuations of the gravitational field lead to the 
conclusion that space-time could lose its continuum character at the Planck's scale. Sev- 
eral model have been proposed to account for this striking conclusion, for example the 
existence, at the Planck's scale, of virtual black holes |59 |l()()|lfiH l()2j: a l so the dynamics 
of string and branes could cause the appearance of a foamy nature of the space-time at 
microscopic scales |63 ( l64 ( l65] . In all these cases the space-time could act as an effective 
environment and all the physical systems would necessarily interact with it. Since the 
degrees of freedom characterizing such an environment are not accessible, any physical 
system should exhibit an open system dynamics, leading to loss of quantum coherence 
(decoherence) [HniESElESISIlEIlIZSlESlEllESlEIl 

The microscopic details of the interaction between physical systems and gravitational 
bath are unnecessary to parametrize this dissipative evolution: we adopt a phenomeno- 
logical approach, that is independent on the underlying fundamental dynamics. Indeed, 
the form of the reduced dynamics can be uniquely fixed by the requirements of being a 
completely positive, Markovian, entropy increasing evolution. 

The Markovian hypothesis is justified by the fact that the action of the gravitational 
thermal bath on the neutral mesons system is very weak. The magnitude of the dissipative 
parameters may be estimeted to be of order E^/Eb, where E$ is the characteristic energy 
of the system and Eb is the fundamental scale for the bath. In our case E$ = wb ~ 5 GeV, 
the neutral B mesons mass, and Eb = Mp ~ 10 19 GeV, the Planck mass (the scale at 
which the space-time is expected to exhibit a kind of granularity). The magnitude of the 
dissipative parameters has thus an upper bound of about 10 -18 GeV. The Markovian 
dynamics is justified, since this scale is very small respect to the scale characterizing the 
usual Hamiltonian dynamics of the B mesons system. 



Remark 3.2 The gravitational thermal bath should affect any physical system. However, 
its effects on the majority of these systems would be completely negligible, so the systems 
characterized by energies in the sub-nuclear range or higher are the only that could feel 
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it. For example, for atomic phenomena, characterized by energies of about 100 KeV, the 
dissipative parameters would be as small as 10 -27 GeV. 

The hypothesis of entropy non decreasing evolution is based on the choice of time- 
evolutions leading to states less-ordered in time. The correlations contained in the coher- 
ence of states are lost as time passes. 

The accuracy of present experiments on neutral B mesons is very high, so that one 
expects the dissipative effects induced by gravitational phenomena to be in the reach of 
present or planned experiments. This is certainly one of the motivations that justify the 
coming considerations. 



3.3 The dissipative evolution 

The equation of motion for the state p representing the B mesons is assumed to have the 
form of a quantum dynamical semigroup, with a non standard contribution Lq in the 
generator, due to the gravitational bath: 

p(t) = -iHp(t) + ip(t)rf + L D [p(t)\. (3.10) 

with H as in (|3.1|) and expressed in (|2.6|) . 

Working in the basis of physical states {|B ), |B )} and using the vector notation 
introduced in the previous chapter, we have 



d_ 
dt 



\p(t)) = (H + V')\p(t)) 



(3.11) 



with 



H 



( 2Im#i 






2 Im H 2 

-iH 3 

m 



-iH A 
i{H* 2 - Hi) 




-iH 3 




(3.12) 



i(H{ - H 2 ) j 

the matrix containing the standard Hamiltonian terms, the masses being redefined to take 
into account possible Hamiltonian dissipative contributions, while the matrix containing 
the non Hamiltonian contribution, causing decoherence of states and entropy increase, is 



V' 



( 


—a 


a 


-a 


-C* 


\ 




a 


—a 


a 


a* 




-C* 


C* 


-A' 


B' 




V 


-C 


a 


B'* 


-A' 


) 



(3.13) 



where 



A' = 7 + a, B' = a. - 7 - 2i(3, C = c- ib. (3.14) 
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The matrix in (|3.13l) differs from that appearing in Q2.14[) because here we assume the 
dissipative contribution to have the form D in (|2.8j) . with u = v = w = 0; V is obtained 
from T> in Q2.14JI by a change of basis, from the eigenstates of CP to the basis of physical 
states. 

The dissipative parameters a, b, c, a, (5 and 7 are small with respect to the parameters 
characterizing the Hamiltonian contribution; then it is possible to apply the formalism 
developed in Section E3J The only difference is the use of T>' instead of T>. The parameters 
a, A' , B' and C take the place of 7, A, B and C respectively. 
In the basis {|B^), |B#)}, that diagonalizes H, the statistical operator becomes 



p d = R-ipR- 1 !, R 



1 -1 



while ()3.11|) takes the form 



d 



j-y{t)) = {n a + v' d )\p d (t)) 



(3.15) 



(3.16) 



where 7i d is a diagonal Hamiltonian part, see (|2.45j) . and T>' d is a more involved dissipative 
part: 



/ A' 



V 



id 



\ 



A' 

-A'* 
-A' 



A' 

9J 

&* 



A'* 
$'* 

& 

n' 



\ 



(3.17) 



the coefficients A', S', S', A', 0' and are expressed as the corresponding parameters 
in Appendix EI with the substitutions 7 — > a, A — > A', B — > B' and C — > C . They are 
linear superpositions of a, b, c, a, j3 and 7 with coefficients that depend on and rn, 
or, equivalently, on 9 and a. In the next section we will show that 9 and a are related to 
violations of the discrete symmetries, that are expected to be small; it is then justified to 
neglect the terms where these coefficients multiply the dissipative parameters. We finally 
find 



A' 
A' 

n' 
e' 



i{Re[(a- 7 + 2i/?y] - (a + 7)}, 
-A', 

-^{2 Re[(c - ib) v 7 ^] + i Im[(a - 7 + 2i/3)a] } , 

-A', 
-A', 

+ + 2a + Re[(a-7 + 2i/3)c7]}, 
--|2a- (a + 7) -Re[(a-7 + 2i^)fj] +4ilm[(c-^) v / ^]}. 



(3.18) 
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Since we have request that entropy be non decreasing, the remaining coefficients in 
Appendix lAl vanish, (see Section f2.3[) 

By means of the coefficients (|3.18|) we can compute the time evolutor Of acting on our 
states; see Appendix IbI 

!/&>(*)) = ^IPbo(O)), l4o(i)) = ^I4o(0)) (3.19) 
with the initial conditions 

p B o(0) = |B )(B |, Pb o(0) = |B°)(B |. (3.20) 

The operators p^ (t) and p'Lo (t) so obtained will be used to compute the decay rates 
of the neutral B mesons. 

As alredy observed in E,emark 12.101 it is convenient to perform all the computations 
in the basis of the eigenstates of H. 



3.4 Invariance properties 

Any symmetry transformation on our system can be modeled by a unitary or antiunitary 
transformation acting on the basis vectors |B°) and |B ). In our formalism, we have to 
consider the corresponding transformation on the vector \p), realized by a 4 x 4 matrix IA: 

\p)^U\p). (3.21) 

The request of invariance of the formalism under the transformation expressed by U has 
the form 

UiH + V')!}- 1 = (H + T>'), UCH + V'p- 1 = {H + V') ] (3.22) 
for unitary, respectively antiunitary, transformations. 

For the discrete CPT, T and CP transformations the form of the matrix U is 



(o :)■ Mo Mo <-> 

where the cr, are the Pauli matrices and the 2x2 identity and 




(3.24) 



contains the phase appearing in the CP transformation of the physical states, eq. (|3.4|) , 
In the standard case T>' = 0; the conditions of invariance under these discrete transforma- 
tions reads 



(3.25) 
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for the CPT symmetry, 

\r L r H \ = 1 => |er| = 1, £ = (3.26) 
for the T symmetry; the CP invariance is obtained imposing both CPT and T invariance, 

tl = r H = VZ, a = e~ 2i v. (3.27) 

Remark 3.3 The parameters a, £ and characterize possible violations of the discrete 
symmetries during the time evolution of our system; they are called violations in mixing 
to distinguish them from violations of the same symmetries that can occour during the 
decay of the B mesons and which are due to the decay mechanism (see Section [3. 5|) , 

Also the dissipative part T>' can be resposable for violations of the discrete symmetries; 
choosing Ti = the conditions of invariance are 



c = b = 
C 1 = e 2ilp C'*, B'* 



c 



D' 



(3.28a) 
(3.28b) 
(3.28c) 



holding for CPT, T and CP transformations respectively. 

We consider now the invariance properties of our formalism under a redefinition of the 
phases of the two physical states. If the states |B ) and |B ) undergo an independent 
change of phase, 



IB 



|B°) 



correspondingly the 4- vector 



-+e^|B°), |B u )^e\ 
undergoes the transformation described by 

\ 



(3.29) 



Us 



( 1 
1 _ 

\ e -n<p-w j 
that causes a phase shift of the off-diagonal entries of the density matrix p: 



(3.30) 



P3 



Pi 



'PA- 



(3.31) 



This transformation is not an invariance for TC; indeed, the off-diagonal elements #3 
and H4 of the effective Hamiltonian H do change. In particular rjy, rt are found to 
transform as 



tl 



(3.32) 
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consequentely, 

a e 2i ^-®o (3.33) 

whereas £ and 9 are invariant. The phase redifinition is not an invariance also for T>': the 
dissipative coefficients are found to transform as 

C> e-^^-^C", B' -» e 2i ^B' (3.34) 

whereas a and A' are invariant. Then equation 1)3.11)1 is not invariant under although 
the physically meaningful quantities (i.e. mean values of observables) are. Some problems 
can arise when this equation is solved by perturbative methods: the physically meaningful 
quantities could exhibit an explicit dependence on the phases <j) and <f>. 

However, combining the results 1)3.32 J) and 1)3.34)1 it is possible to show that the entries 
of the matrix T)' d are rephasing invariant. The same is true for the matrix 1rl d since 
it is diagonal. These results ensure the invariance of our formalism under any phase 
redefinition, also for a truncated perturbative expansion of the time evolutor. This is 
another advantage in the use of the eigenstates of the effective Hamiltonian as basis of the 
two-dimensional Hilbert space. 



3.5 Observables 

In this section we want to introduce the observables usually used to describe the physics of 
the neutral B mesons system. In the formalism of density matrices, these are represented 
by Hermitian operators Of related to the decay of the B mesons into a generic final state 
/. In the |B°), |B°) basis, these operators are represented by 2 x 2 matrices: 

Ox 3 
o 4 2 



°f= n n ( 3 - 35 ) 



whose entries, by construction, are related to the two independent decay amplitudes 
A(B° — ► /) and _4.(B° — > /): 



Ox = \A(B° - f)\ 2 , 3 = [A(B° - f)]*A(W - /), 
2 = \A(B° - f)\ 2 , 4 = [A(B° -» f))*A(B° -+ /). 



(3.36) 



with 



A(B°^/)^(/|T|B ), A(W^f) = (f\T\W), (3.37) 
T being the scattering operator. 
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Accordingly to Q1.3[) . the mean value of Of is obtained by a trace operation over the 
state describing the system at time t, regardless of the basis of the Hilbert space considered: 

(O f (t)) = Tr[0 /P (i)] = Tr[O d fP d (t)}, (3.38) 

where O f is t he matrix representation of Of in the |B^), |B#) basis: Of = R^OfR and 
R as in (|3,15j) . This mean value gives the probability rate for the decay of the B meson 
into the final state /. 

In the analysis of the decays of neutral B mesons, it is customary to introduce the 
following rephasing invariant phenomenological parameters [571158] : 



A 



f _ q L A(B° -> f) f= l_ 
L ~ PL A(B° — >■ /) ' ^- A / 

f _ q H A(B°^f) f J_ 
By means of them the matrix Oi can be written as 



(3.39) 



O d f = \A(B°^f)\ r < 



l^(B°^/)| : 



|r L | 2 |l + A{| 2 rHrl{ i-X f H )(l + X f L r 

r L r* H (l + \{)(l-\ f H )* |r H | 2 |l-A^| 2 

-(i+ M {)(i-M^r ii-4i 2 



(3.40) 



A class of decays widely studied is the set of the so-called semileptonic decays, namely 
the decays where the final state / is constituted by a charged hadron h, a charged lepton 
/ and its associated neutrino v. B°, B -» hrl + v, h+rv. In this case, it is convenient to 
parametrize the decay amplitudes as follows: 

*4(B° — > h~l + v) = MhO- — Vh), (3.41a) 

A(B° ^h + rv)=M* h (l + y* h ), (3.41b) 

^(B° -» h + l~V) = ZhA(T? -» /i + / - ^), (3.41c) 

„4(B° -» /TZ+i/) = x ft „4(B° /TZ+z/), (3.41d) 

where is a common factor. Possible violations of CPT invariance and of the AB = AQ 
rule are parametrized by the coefficients Xh, Uh and Zh- In particular, non vanishing Xh and 
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Zh indicate violations of AB = AQ, because this rule would forbid the decays B° — > h + l~V 
and B — ► h~l + v, while a non-vanishing yh signals violation of CPT invariance, since this 
symmetry implies equal probability for the processes B° — ► h~l + v and B — > h + l~V. 

We assume the violations of the discrete symmetries (in mixing or in decay) and of the 
AB = AQ rule to be small. In this way we can neglect all terms in wich the coefficients 
6, £, Xh, Vhi z h are multiplied with the dissipative parameters a, b, c, a, (3 and 7. With 
this assumption, for the semileptonic decays, the parameters introduced in (|3.39f) assume 
a simplified form: 



Al" = \ h H = V^x h = Xh, nf = ii h u = \fcz h = n h . (3.42) 

Another class of decays we shall examine are characterized by a final state / that is 
an eigenstate of CP, with a definite CP parity Q: 



CP\f) = Cf\f), C/ = ±l- (3-43) 
At lowest order one finds L = Q [BT] . 



3.6 Decay rates for a single B meson 

In our formalism, the probability rate for the decay of an initial B° or B° meson in a final 
state / at time t is given by 

V f (B°;t) = Tr[OjpUt)}, V f (B°;t) = Tr[Ojp^(t)}. (3.44) 
Together with these probability rates, it is useful to define integrated rates as well: 

Vf(B ) = -J Vf(B°;r)dT, V f (B°) = - J V f {B°;r)dr. (3.45) 

As an example, we will explicitly list the probabilty rates for the semileptonic decays. It 
is convenient to define the parameters 

Am AT . . 

, S=— (3.46) 

where Am = mn — mL, T = (7^ + 7l)/2 and AT = ji — ^h] it is also convenient to 
rescale the time as r = tT (that is, we measure the time in units of B lifetime). Note that 
in the case of B^ mesons 5 is expected to be small, of the order of 10" 2 ; this allows to 
simplify some expressions in the coming sections. 
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and B° 



Explicitely, one finds for the four relevant semileptonic decays B° 
h~l + v, h + l~v; for example, 



hrl+v, h+l- 



r h -(B°;r) 



\Mh\ 2 e r |sinh5r 



+ cosh 5t 



sinwr 



+ COS LOT 



\ru\ 2 ~ Vl\ 2 _ 
\th + r i\ 2 
\ru\ 2 + \tl\ 2 
\th + r h\ 2 
Im.(rnr* L 



Re \h 



AT 



+ 



2 5 

Rey h + - — -Re A' 



L \rH + r L \ z 
Re(r L r* H ) 
\th + r L \ 2 



T5 2 +co 2 
Re8' 



+ 



2A?n 



+ 



e -(H'+A')f _ Re 



Vh 



T6 2 + uj 2 



Re A' 



(3.47) 



and the others are reported in Appendix O 

In these expressions we can identify the contributions coming from the standard evo- 
lution, from the hypothetical violations of CPT and AB = AQ and from the dissipation. 
The probabilities Vh are directly accessible and can be studied in experiments performed 
at colliders: a fit of the time dependence of these quantities should permit, in line of prin- 
ciple, the evaluation of dissipative parameters A', A', £1' and ReO'. In the next section 
we shall define some quantities, called asymmetries, more suitable to extract informations 
about the dissipative parameters. 

Analogously we can compute the probability rates for the decay of B° and B° into final 
states with a definite CP parity, for example D + D~, ir + ir~ , J/\I/ K. Their expressions 
are more involved and can be found in Appendix [O 



3.7 Asymmetries for a single B meson 

In order to reduce systematic errors, two experimental studies of the neutral B mesons 
system are usually performed using suitable combinations of the decay probabilities intro- 
duced in the previous section. They are called asymmetries and give a direct information 
about interesting physical quantities. 

Given the process B° -» B°, its CPT conjugate is B -» B ; given instead the process 
B° -> B its T conjugate is B — > B . Since, in absence of violation of the AB = AQ 
rule, h + and h~ signal the decay of B and B° respectively, CPT and T violations are 
studied by the asymmetries 

_ 7V(B°;t)-7V(BV) 
Acpt(t) = — — — -, (3.48) 

7V(B U ;t)+7V(B°;t) 
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respectively 



A t (t) 



(B°;r) 



7V(B°;r 



n-(B°;r)+n + (B°;r) 

In the standard case (i.e. without dissipative evolution) they vanish if CPT and T are 
preserved. In our case, using the expressions in Appendix we find 



(3.49) 



A C pt{t) 



(l r i| 2 — l r i?| 2 ) sinh^r + 2Im(rf,r|^) sinuT 
{\tl\ 2 + l r H| 2 ) cosher + 2~Re(rLr* H ) coswr 

[Im(A^ — jj,h) sinwr + Re(Afr 



+ 



cosh 5t 
cosh 5t 



COS UJT 

cos UJT 4 



^ ft )sinh*r]+ ( 3 - 50 ) 



cosh St + cos ujt T 5 2 + uj 2 



Re A' 



and 



A t (t) 



r L r H 



1 



1 

+ 



2 Re y h + 



cosh 5t — cos ujt 



sinh <5r 



+ sin ujt 



lm(X h - fi h ) - - 



Re(/i h 
6 



ts 2 + uj 2 



ImA' 



+ 



(3.51) 



T5 2 +lo' 2 



ImA' 



We observe that these quantities are affected by the presence of dissipation, since 
they contains both the real and imaginary part of A'. However, in their expressions the 
dissipative parameters appear multiplied by the coefficient <5, that, as alredy noticed, is 
expected to be small; the dissipative effects are suppressed with respect to the standard 
contributions. As a consequence, these asymmetries are not in general suitable for probing 
the presence of dissipative effects in the neutral B mesons system. 

To isolate the dissipative contribution, more complicated combinations of the proba- 
bility rates are needed. For example, the asimmetry 



7V(B°;t 



7V(B°;t) 



7V(BV)-7V(B°;t 



TV (BO; r) + V h - (B°; r) + V h+ (B°; r) + V h+ (B°; r) 



(3.52) 



used in the standard case to fit the mass parameter Am, in our case presents a charac- 
teristic dependence on the dissipative terms; in the limit of a vanishing 6, it takes a very 
simple form: 



AAm( T ) = (R e $) 2 (1 ~~ cos UJT ) + e n r cosur — sinwr 



1 

Am 



Re6' + lm{X h + fi h ) 



(3.53) 
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Assuming the validity of the AB = AQ rule (that forbids the decays B° — > h + l v and 
B — > h~l + v), and neglecting the quadratic dependence on the CPT-violating parameter 6, 
the different time behaviors in A/^ m (r) should allow a determination of the non-standard 
parameters £1' and the real part of 0'. The actual accuracy of such a determination 
highly depends on the sensitivity of the measure of -A Am ("?")• From the results of available 
simulations one can expect an accuracy in the determination of the various terms in 

mir) °f sbout a few percent |79|I80| . This sensitivity is alredy enough to give interesting 
bounds on fi' and Re 0'. 

The assumption 5 = has some important consequences. Indeed, in the limit AT = 
the non-Hermitian part r of the Weisskopf-Wigner Hamiltonian H becomes proportional 
to the identity. In this case the eigenstates of H, |B#) and |B^) becomes orthogonal, 
so rLr* H = 1, producing |c| = 1 (or £ = 0) and, accordingly to (|3.9[) . Im.9 = 0. Then 
T violation in mixing is forbidden and CPT violation in mixing is parametrized by Re#. 
This fact has important consequences for experimental tests of CPT and T invariance in 

mixing |miEllE3IBllE51ll>l- 

Another interesting asymmetry is defined using the decay of both B° and B in an 
eigenstate / of CP: 



Mr) 



V f (B G 



T f (B°;T)+T f (& > ;T) 



(3.54) 



in the standard case, this observable is particularly suitable for studing violations of CP 
invariance. In presence of dissipation, we get 



Mr) 



2ReA f „ „ 2( f T A , 

=- — ^Re# + ^ImA'+ 

1 + \Xf\ z loT 



sin ujt 



+ COS LOT 



21mA 



l + |A f l 2 



/ +^ReA' 



\X f \ 2 



l + |A f | 2 



2_Cf 
2 Re A/ 



+ 



(3.55) 



Re6- ^ImA' 



where the condition 5 ~ has been used; this enables us to write A^ and A£ in terms of 
a unique amplitude ratio A f : 



A/ 



: A(1? 



f) 



A(&> - /) 



(3.56) 



and 



A{ = A/(1-Re0), A^ = A/(l + Re9). (3.57) 

The presence of dissipation in (|3.55|) manifests itself in the modification of the coeffi- 
cients of the oscillating terms and in the presence of a r-independent piece. Dedicated B 
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experiments at colliders should be able to identify the various r-dependence and therfore 
provide stringent bounds on the parameter A'. 

Asymmetries analogous to the ones just discussed can be constructed using integrated 
rates, see eq. (j3.46|) . rather than the probability rates. The quantities we obtain in this 
way are time independent and are denoted with a prime index to distinguish them from 
the previously defined ones. Two examples are 



A' Ar 



1 + 



l + w 2 (Re#) 2 -u 



Am 



Ree / + Im(A /i + ^ h ) 



+ 



r( 



I l^ y [(^-l)0' + 2^A']} 



(3.58) 



and 



2 Re A 



l + \X f 
1 



? Re0 + %ImA'+ 



l+UJ 2 

1 — lAv 



21mA 



+ 



1 + |A/| : 

2 Re X f 



+ 



2C/ 



Re A' 



+ 



(3.59) 



1 + IA/P 1 + |A, 



Re0 + ^ImA' 



The result in equation Q3.58JI is particularly interesting: assuming the validity of the 
AB = AQ rule and CPT-invariance, a measure of A' Am (l + w 2 ) not compatible with 
1 whould signal the presence of non standard effects. Such a test can be performed at 
colliders with high accuracy, in particular using the next generation dedicated B mesons 
experiments. 

The integrated asymmetry (|3.59|) is not very useful in probing the presence of dissipa- 
tive effects in neutral B mesons decays: since there is no time-dependence, the measure 
of this asymmetry would give little information on the magnitude of A'. 



3.8 Decay rates for entangled B meson 

In the analysis of neutral B mesons physics, the study of pairs of mesons produced at the 
so-called B-factories, that is e + e~ colliders tuned at the energy of the T(4S) resonance, is 
of particular relevance. The pairs of B mesons are produced by the decay of this resonance, 

T(45) -► B°B°, (3.60) 

in an entangled state. Since T(45) has spin 1 whereas the two B mesons are spinless, 
they are produced, in the T(45) rest frame, in an antisymmetric spatial state. In this 
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reference frame B° and B are flying apart with opposite momenta p and — p; the state 
encoding this configuration is the singlet state: 

\iP) = ^|(|B , -P) \B°,P) - |B°, -p) ® |B°,p», (3.61) 

or, in the density matrix formalism, the operator p = in this case a 4 x 4 matrix. 

Its irreversible time evolution is assumed to have a factorized form: 

p(ti ) t 2 ) = (A tl ®At a )b(0)], (3.62) 

where we consider two possibly different proper times t\ and ti for the two mesons since 
in general they won't decay at the same time. 

Remark 3.4 In the standard treatement of the neutral B mesons system the time- 
evolution is assumed to have the factorized form (|3.62ft since there is not interaction 
between the two particles in the pair; in the open system approach an interaction can 
be induced by the common thermal bath, but it disappears at the lowest order in the 
dissipative effects; whence the hypothesis of factorization. 

As in the previous sections, we prefer to express the statistical operator in the basis 
diagonalizing the Hamiltonian contribution, then 

P d (h,t 2 ) = (ir 1 <g> R~ 1 ) P (t 1 ,t 2 )(R- 11 <8> ir lf ), (3.63) 

which in terms of the basis {Pi, i = 1, ... 4} introduced in Section \2. 11 reads 

P d (h,t 2 ) = 77, 7 7 2 -{Pl(tl)®P2(t2)+P2(tl)®Pl(t2)-P 3 (t 1 )®P 4 (t 2 )-P A (t 1 )®P 3 (t 2 )), 

WL + THY 

(3.64) 

where P;(i) = A t [Pi]- 

The probability that a meson decays into a final state f\ at proper time t\ while the 
other meson decays into the final state / 2 at proper time t 2 is given by 

V(f 1 ,h-f 2 ,t 2 ) = Tr[(O d fl ®O d f2 )p d (t 1 ,t2)] (3.65) 

where O d x and O d 2 are the 2x2 matrix representations, in the basis of eigenstates of H, 
of the Hermitian operators Of 1 and Of 2 described in Section 1531 

Since the neutral B lifetime is very short, it is difficult to identify the two proper times 
t\ and t 2 : the probability (|3.65|) is experimentally unattractive. It is thus customary to 
define integrated rates, as in the case of the decays of a single meson. In particular it is 
possible to consider a fixed time delay t = t\ — t 2 between the two decays, obtaining a 
one-time probability rate: 



/■oo 

v(f 1 ,f 2 ;t) = / p(/ 1 ,^ / + ^;/ 2 ,^ , )^ , 

J 



(3.66) 
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when t 0, that is when the first meson decays after the second, while, for t < 0, one 
defines 



J o 

where the 6 of Heavyside is required since the evolution is of semigroup type, with forward 
in time propagation. This definition implies 



In the following we shall always assume t ^ 0. It is also possible to define a time 
independent integrated rate as 



in which we sum over all possible delays between the two decays. This quantity is relevant 
in the experimental analysis at the so-called symmetric B factories where the T(4S') 
decays pratically at rest in the laboratory frame, and then the time information is hardly 
attaineable. On the other side, at the so-called asymmetric B factories the decay of the 
T(4S') resonance takes place in a boosted frame respect to the laboratory; in this case the 
temporal information can be extracted and the one-time probability (|3,66[) is suitable to 
the data analysis. 

The explicit expressions for these joint probabilities are rather involved; some partic- 
ularly interesting examples are collected in Appendix IO 

One can nevertheless make a general remark. Consider V(f,f;0), i.e. the probability 
rate for the simultaneous decay of the two B mesons in the same final state /. Since 
the initial state is antisymmetric, the usual dynamics produces V(f, /; 0) = whereas the 
modified dynamics in general gives V(f, /; 0) / 0. This quantity could then unambiguously 
show the presence of non standard effects. 

3.9 Asymmetries for entangled B meson 

As in the study of the decays of single B mesons, the experimental analysis of our system 
is not usually performed directly with the probability rates introduced before; rather, 
some suitable asymmetries are constructed. We use the assumption 5 = to simplify our 
expressions and, as before, we refer to a rescaled time scale r = tT. 
An interesting quantity involving semileptonic decays in both mesons is 




(3.67) 




(3.68) 




(3.69) 




(3.70) 
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In the given approximation, one finds 



K{t) 



(1 + |rffr£,| 2 )(cosh<5r) 



(\ r L\ 2 + | rjy | 2 ) cosh 5t + 2 Re(rLrJ^) coswr 
2 



+ 



+ 



(1 + coswr) 



(1 + TJ— COSUT + — (SIIIWT + U COSLUT) Re 

r Am \ 1 — iu> 



(3.71) 



and its evaluation for simultaneous decays is particularly meaningful: 



*(o) = - 2 



1 



U) 



9J 



r Am 1 + uS 1 



'wlmO'-Re 9') 



(3.72) 



since in the standard case 7Z(0) = 0. Therefore a fit of TZ(t) on actual experimental 
data that produces a non-zero intercept would clearly signal the presence of non-standard 
effects. 

Another quantity built with the double semileptonic decays and relevant in the study 
of dissipative effects is the asymmetry A& m {T), analogous to the one defined in (|3.52j) : 

[V(h+, h--r)+ V(h-,h+; t)] - [V(h+,h+; r) + VjhrjTyr)] 
AmiTj ~ V(h+,h+;T)+V(h-,h~;r)+V(h+,h-;T)+V(h-,h+;T) ' 1 ' 

Its explicit expression is particularly simple and meaningful, even in the case of vanishing 
8: 



r + n> 

— (sinwr 



e n r coswr + (Re6*) 2 (l - coswt) + 



1 



Am 



u) cos ujt) Re 



0' 



1 — iu 



(3.74) 



Thanks to the different time behavior, an experimental study of this observable should 
allow to extract precise information on £1' and the combination wlm9' — ReO'. The 
asymmetry A Am is intensively studied at B-factories [BI1IBB1IBH] since it allows a pre- 
cise determination of the mass difference Am. Indeed, AAm( T ) = cos ujt in absence of 
dissipation. 

Another asymmetry interesting from a dissipative point of view is built with one- 
time integrated rates and it involves a semileptonic channel and an hadronic channel (in 
particular a state with definite CP-parity): 



A f (T) 



V(h+,f;T)-V(h-,f;r) 
V(h+,f;T) + V(h-,f;T) 



(3.75) 
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If we neglect possible violations of the AB = AQ rule, the semileptonic decay identify 
the flavor of the decaying meson. In this case this asymmetry is equivalent to the A fir) 
defined in Q3.54JI . If the final state is / = J/\l/ K these observables are directly expressed 
in terms of one of the angles that parametrize the CP violation in the Cabibbo Kobayashi 
Maskawa matrix; for this reason they are extensively studied. The expression of Af(r) in 
the case of irreversible dynamics is given by: 



AAr) 



2 Re A/ 
: l + |A-i 2 



2C/ 



Re 6 + C f Re(A h - fx h + 2( f y h ) - -%lm A'+ 



+ COS LOT 



+ sinwr 



1 



I \ |2 



21mA/ 



l^Re.-^Re 
1 + |A/| 2 ujT 



uj + 2i 
2 + iuj' 



A' 



2 Si 

LOT 



Im 



u + 2i 
2 + iuj' 



A' 



+ 



(3.76) 



where Xf is defined in (|3.56|) and in ()3.43|) : for simplicity in this formula only first order 
terms in CPT violating parameters are kept. Equation (|3.76|) holds for r ^ 0; for r < 
the expression is 



A f ir) 



2 Re Re + C Re(Ah _ Mfc + 2C f y h ) - ^ Re f ^ A' 



l + |A'l 2 



+ COS LOT 



+ sinwr 



1 + |A/| 2 1 + |A/| 2 wT 



21mA 



/ 



1 + |A/| 2 ' cjT 



+ — ^ Re A' 



+ 



(3.77) 



the two expressions coincide for r = and, in absence of dissipation, one can be obtained 
from the other by letting r — ► — r. The fact that this is not longer true when A' ^ is a 
signal of irreversibility. 

Assuming the validity of the AB = AQ rule and neglecting CPT- violating effects in 
mixing (6 = 0) as well as direct CP violations (i.e. |A/| = 1, see |57j). a fit of (|3.76|) with 
experimental data allows the determination of both Re A' and ImA'. The study of Af 
could result in one of the best tests on the presence of dissipative effects in B physics, in 
particular for the final state J/\I/ K |87| , 

The quantities involving the fully integrated probability rates, eq. (|3.69|) . in general 
are not suitable to study dissipative effects, since the dissipative parameters are multiplied 
by 5; neverthless it is possible to define two ratios unaffected by this problem: 



XB 



V(h + ,h + )+V(h-,h- 



P(h+, h+) + V(h~, h-) + T(h+, h-) + V(hr,h* 



(3.78) 
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and 



Rb 



V(h + ,h + )+V(h-,h-) 
V{h+ : h~)+V(h-,h+)' 



(3.79) 



Assuming the AB = AQ rule, xb represents the total B° - B mixing probability. In 
our framework, one explicitely finds: 



UJ 



XB 



1 - (Re#) 2 + 



Rb 



2(1+ a; 2 ) 

2 [1- (Re0) 2 ] 



UJ 



+ 



UJ 2 (1+UJ 2 

4 



-X 



2 + w 2 [l + (Refl) 2 ] (2 + w 2 ) 2 ' 
the dissipative contribution is the same for the two observables and it is given by 



T v 'Am 



(wlme'-Ree'). 



(3.80) 
(3.81) 

(3.82) 



Independent measures of these two quantities would provide a way to estimate both 
He 9 and X, and therefore give limits on both dissipative and CPT-violating effects. 



3.10 Discussion and conclusions 

In this chapter we have treated the neutral B mesons system in the framework of quantum 
dynamical semigroups. The formalism is very general and can be applied to describe the 
time evolution of this system when it is subject to a weak interaction with an external 
environment. As described in Section l3~2*l quantum gravity could give a motivation for such 
an approach, since in this case the space-time itself would act as an effective environment 
for any physical system. 

As discussed before, the open system framework allows a rough evaluation of magnitude 
of the dissipative parameters; they scale at most as m^/Mp, where me is the meson mass 
while Mp is a large fundamental mass scale. Tipically, Mp coincides with the Planck 
scale so that mf^/Mp ~ 10~ 18 GeV. This value is very small; however, the sofistication 
of the dedicated B-experiments, both at colliders(CDF-II, HERA-B, BTeV, LHC-6) and 
B-factories (BaBar, Belle, CLEO-III), is so high that the sensitivity required to probe the 
presence of non-vanishing dissipative parameters can be reached in just a few years of data 

taking [THIIHOllHZllHHIIHni- 

Remark 3.5 An evaluation of the dissipative parameters for a gravitational bath have 
been alredy performed in the context of the neutral K mesons physics |55U90| . Upper 
bounds for the dissipative parameters have been obtained; these are in the range of what 
suggested by gravitational or stringy effects. More data are expected from ongoing exper- 
iments at the so-called ^-factories, the analogous of B-factories for the K mesons. 
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To conclude, we observe that, besides the study of irreversible dynamics, these experi- 
ments, involving correlated mesons, could give the possibility to investigate the mathemat- 
ical notion of complete positivity from a phenomenological point of view. Let us consider 
a simplified situation with a trivial Hamiltonian sector, H = 0, and a dissipative part with 
vanishing a, b and c. A dissipative contribution of this kind has been used to describe 
non-standard effects in the K meson system, using a positive but not completely positive 
dynamics |fi? H l77 | l78| . The eigenvalues of the evolved singlet state in (|H.62|) , for t\ = ti , are 



where 



A 1)2 = 2 ± {\E+(i)f + [E-{t)] 2 + 2[F(t)] 2 }, 

A 3 ,4 = ±[E+(t) + E_(t)]{[E + (t) - E_(t)] 2 + 4[F(t)] 2 } 1 / 2 , 



E±(t) = + 2a)e" ±1 - + 2a)e^% 

\a (3-84) 
F(t) = _ v [e v - f - e v +% 



with v± = -(a + 7) ± yj(a - 7) 2 + 4/? 2 . 

The condition of positivity of this dynamics reads 07 ^ j3 2 whereas complete positivity 
requires a = 7 and /? = 0. Consequently, in the simply positivity case v± are both negative 
and thus, for t ^ two of the eigenvalues (jH.HMjl are always negative and the statistical 
interpretation of the density matrix p(t, t) is compromised. In the completely positive case 
we find v + = and all the considered eigenvalues remain non negative for any time. 

These results offer a further motivation for studying generalized dynamics of the 
form (|3.10|) at mesons factories. These set-ups are in fact high-performance quantum 
interferometers: at least in principle, they can verify from the experimental point of 
view the role of the condition of complete positivity in the time evolution of correlated 
mesons |;44 , 45 , 46 : 47 ; 48 : 4^] - As shown by the simple example above, this is not just a mere 
technical question: it is crucial for our physical understanding of the quantum dynamics 
of open systems. 
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Chapter 4 

Neutrinos 



Neutrino physics is presently a very active area of reserch. Recent experiments suggest 
neutrinos to be massive and the three flavors, electronic (z^e), muonic (zAj) and tauonic 
{v T ) neutrinos, to mix among themselves. Consequently the Standard Model should be 
modified by the introduction of a mixing matrix in the leptonic sector (analogous to 
the Cabibbo Kobayashi Maskawa matrix in the hadronic sector). The first hint to the 
possibility of neutrino mixing came from the observation of a reduced v e flux from the 
Sun, incompatible with the usual model of solar evolution. The most favoured explanation 
of this phenomenon is the so-called MSW effect which is deduced from the hypothesis that 
electronic neutrinos interact with solar matter [91II92II93| . In this context, flying neutrinos 
behave as open systems in interaction with an environment (the Sun), but the MSW effect 
is only the Lamb-shift term (a redefinition of the neutrino effective Hamiltonian) and no 
dissipative effects appear. 

In other approaches, dissipation has been taken into account [9*Hl9*o1l9*6ll9ni9*51l9*9lllUUl 
ll01|ll02lll03|IT0l|ll05j ■ though paying scarce attention to internal consistency. In this 
chapter, following the lines introduced in Chapters ^ and [21 we will apply the rigorous 
formalism of open quantum systems to describe neutrino time-evolutions in matter, with 
particular care to the derivation of a Markovian reduced dynamics. In this way, we will 
get results unaffected from physical inconsistencies that otherwise can appear |99j . 

We shall consider an environment with constant mean density and derive the general 
form of the generator of the time-evolution in the Markovian approximation. Since, in 
the standard description of the MSW effect, it is the variation of this mean density that 
produces the most relevant aspects in the treatment of the flavors oscillations }92(I93| . we 
shall describe how to generatize our formalism to account for the case of environments 
with a varying density profile. 

The transitions probabilities are explicitly written out in two different situations: a 
flavor-changing interaction and a flavor-preserving one (the case of solar matter). 
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4.1 The physics of neutrinos in vacuum 

In this section we shortly outline the principal facts of neutrino physics in vacuum. For sake 
of simplicity, we consider the case of two neutrino's flavors, v e and and we describe the 
neutrino system by means of the bidimensional Hilbert space of flavor, 9)2 |106U107lllU8| : 
indeed, in vacuum the neutrino momentum can be kept fixed and the related degrees 
of freedom does not appear in the description. This is the standard phenomenological 
approach to this system and the experimental data are usually analyzed in this context. 
This approach is justified by more rigorous treatements; see 

As in the case of the neutral B mesons system we distinguish between two basis for 9)2- 

• the basis of physical states (or flavor states), {|^ e )> these are the states that 
are actually detected in experiments; 

• the basis of the eigenstates of the Hamiltonian of the system Hg, {l^i), |^2)}> char- 
acterized by 



where the H$ is Hermitian since neutrinos are stable particles and 2 are real 
eigenvalues. 

Hs generates the time-evolution of egenstates in vacuum: 



and the eigenvalues e- are the energies of these states, that depend on the neutrino momen- 
tum p and mass mj. Indeed, in agreement with experimental evidence |1U9) . the states v\ 
and V2 are assumed to be massive and we denote with mi, respectively m2, their masses. 
In the Standard Model they vanish; for later reference we assume they are not degenerate. 

Since mi and m2 are believed to be very small, e[ can be expanded in the ultrarela- 
tivistic limit, i.e. neglecting 1 the terms of higher order than m\ and m\: 



where E ~ p is the neutrino energy. 

The two basis are connected, in vacuum, by an unitary transformation Uq (the so- 
called mixing matrix) that can be chosen to be orthogonal by a suitable redefinition of the 
states and \v2)- 

1 These terms are supposed smaller than the dissipative effects that we shall consider in the cominig 
sections 




(4.1) 




(4.2) 




(4.3) 
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Remark 4.1 The choice of an orthogonal Uq matrix is a consequence of the impossibility 
to violate the CP symmetry in a two flavors context: in fact the CP violation is related 
to the existence of unavoidable complex phases in the mixing matrix and they appear 
necessarily in a three flavors context. However, when dealing with dissipative dynamics, 
in general the complex phase appearing in Uq is relevant also in the two flavors context, 
because it produces physical consequences |117j . In our context, however, this phase can 
be consistently reabsorbed; indeed, the particular dissipative effects we shall consider in 
what follows does not produce any evidence of this phase. Therefore, in the following, we 
do not account for possible violations of CP in the leptonic sector. 

We have: 

v e ) = C/oki), 



Uq\v 2 ); 



cos V sin 6 
Uq = [ . p ; 4.4 

- 8111 9 COS f 



where the parameter 9 is the mixing angle in vacuum between the two families considered. 
The Hamiltonian H$ in the basis above reads: 

rd _ ( £ 1 \ IT _ ( ~ COS 20 Sin 20 



Hs= {o e 2 J' Hs = ^e{ sm29 cos 29 ) (45) 

where 2 H$ = UoHgU^ 1 . Since we prefer to deal with a traceless matrix H&, we have 
redefined the eigenvalues of Hs as e% = —Am 2 /AE and e 2 = Am 2 /AE, where Am 2 = 
m 2 ~ m i an d E is the neutrino energy. 

The distinction between the two basis of the Hilbert space produces the phenomenon 
of flavor oscillations, i.e. periodical transitions between v e and z^. 

Since a physical state is given by a superposition of eigenstates of Hs, an initially pure 
state \v e ) will present, for t > 0, a non-zero component along \ For instance 

A 2 

\ Ve (t)) = ( e iAm2t / 4E cos 2 9 + e -^m 2 i/4£ gin 2 + f is i n 29sm=^t)\^). (4.6) 



The probability of transition (in vacuum) is given by |106[I107|I1()8] : 

Pu^At) = \{u^\u e {t))\ 2 = sin 2 29 sin 2 -^-t (4.7) 



4E 



and the surviving probability is 



P Ve ^ Ve (t) = l-P Ve ^ tl (t). (4.8) 
Vacuum neutrino oscillations were the first explanation of the lack of v e in the solar flux. 



2 We use the notation introduced in Chapter |5]to characterize the operators expressed in this basis 
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Remark 4.2 This phenomenon should lead to lepton numbers violations, forbidden in 
the Standard Model. 

However, the transition probabilities are modified by the interactions of neutrinos in 
matter |91U130| : in a homogeneous medium a matter mixing angle and two redefined 
masses are needed and the probabilities of transition are formally equal to 1)4.7)) and 1)4. 8J) 
but with the redefined parameters 6 m and Am^,. These transition probabilities are used 
to fit the parameters 9 and Am 2 to the experimental data )106U107lll08j . 

In what follows, we shall describe neutrinos in flavor space by means of statistical 
operators since we are interested in decohering evolutions. In particular, p Ve = |z/ e )(z/ e | 
and p Ufl = |z/ /J )(f At |, whence 1)4.7)1 reads 



We shall also adopt the vector representation, developed in Chapter [2J for the density 
matrices. 

4.2 The model of interaction 

In the previous section we have described the behavior of neutrinos in vacuum. Here we 
consider the case of neutrinos interacting with a suitable environment, consisting of a slab 
of matter of volume V. The Hamiltonian describing both neutrinos and environment has 
the form Ht = H$ + Hb + Hi, where H$ is the neutrino Hamiltonian (|4.5|) . Hb is the 
Hamiltonian of the environment, and Hi is the contribution containing the interaction 
between neutrinos and environment. In this section we shall derive its form. 

Remark 4.3 In order to have a general Hi, it is not possible to consider a two-dimensional 
$)S, as i n vacuum; momentum dependence needs also to be taken into account. Neverthe- 
less, we shall see that, under certain circumstances, a two-dimensional effective description 
is available also in the open system context. 

The interaction Hi is determined by a Hamiltonian density fj by 



V 

where x = (x, t) (the dependence on time in Hi disappears if we suppose the layer of 
matter to be uniform). We assume a linear coupling of the form 



P Ve ^(t) = Tr[p Ue (t)p Vfl ]. 



(4.9) 




(4.10) 




(4.11) 



o 



with 




(4.12) 
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where g is a suitable coupling constant (dimensionally it is an inverse energy squared), 
Vg B (x) are the neutrino and bath currents that characterize the interactions and a is a set 
of flavor and Dirac indices. The 7«'s, /i = 0, . . . 3 are the Dirac matrices and 75 = ^70717273 
while vi(x), I = e, fi, are the fermionic neutrino fields. To guarantee the Hermiticity of Hj 
we require 

V%*\x) = V^"{x), (4.13) 

otherwise, V^(x) is completely arbitrary. The indices S and B denote the neutrino system 
and the environment respectively, according to the notation of Chapter^ 

Remark 4.4 The assumption (|4.12|) is quite general; for instance, it includes transitions 
forbidden in the Standard Model, like neutrinos flavor-changing interactions |112II113II1141 
I115j . These latter are induced by suitable environment operators V 1 ^ 1 ^ with I m. 

Remark 4.5 To describe neutrino interactions in the Sun, a low-energy Hamiltonian den- 
sity is actually introduced; it contains neutral currents and charged currents contributions 
in the point-interaction approximation (Fermi theory) to the Standard Model jll)8j : 

\){x) = ^J=Y1 jy^h^ 1 - 7sMM ® [6i e e(x)^(l - 7s)e(x)+ 

1 " , (4-14) 

+ [/ 3 (1 -75) -2Qsin 2 ^] /(*)}, 

f 

where I = e,/i identifies the neutrino fields; / = e,p,n characterizes the fermionic fields 
(electrons, protons and neutrons respectively); I3 and Q are the isospin and charge opera- 
tors; 9w is the weak angle and Gf the Fermi constant. In this case, only flavor-preserving 
transitions are allowed. 

In the representation for which 75 is diagonal, the left component of neutrino field can be 
written as 

^) i = i(l-75Wx) = ^( (4.15) 

correspondingly the neutrino currents become 

V^(x) = 2<pl(x)a^ m (x), (4.16) 

where {a^, fj, = 0, . . . , 3} is the set of Pauli matrices plus the 2x2 identity. The Weyl 
spinors 4>i(x) can be expanded only in the sector of interest over components of negative 
helicity and positive energy |116j so that at t = 

<Mz, 0) = JL £ e ^(3(k)Mk), (4.17) 

k 
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where V is the volume of the layer of matter crossed by the neutrino, and (5{k) is the 
negative elicity eigenstate, 

a -^m = (4.i8) 

4>i(k) are the annihilation operators of neutrino states with flavor I and momentum k and, 
together with the creation operators 4>\{k), they satisfy anti-commutation relations: 

4(k)\0) = \l,k); {Mk)Aj(k')} = S n%n {Mk),Mk')} = {<t>j(k),<f>j(k')} = 0. (4.19) 
Using (j4.17j) the currents (|4.12j) read 

V s a (x,0) = i^ e ^'-^%(fe,fc04(^)^(^), (4-20) 

where a = (l,m,p) and S^(k,k') = ft (k)<j^f3(k'). The behavior of this quantity under 
complex conjugation is 

[S„(M')r = S^(k',k). (4.21) 

These quantities shall be employed in the following sections. 

The generic neutrino state ps in the finite volume V can be written in terms of flavor 
and momentum eigenstates: 

^(*) = EE^( r >"*;*)i^O s i = ^2Y2 PtJ (f,s;t)4(r)\o)mj(^, (4.22) 

where the creation and annihilation operators satisfy Q4.19j) . 

4.3 The Markovian approximation 

Following the standard steps introduced in Chapter ^ in order to deal with the environ- 
ment and to end up with a Markov time-evolution for the neutrino state ps, we have to 
perform a suitable approximation. 

In Section ll.4l we have described two different procedures to obtain a Markovian ap- 
proximation for the evolution of a system in a thermal bath. Introduced the typical times 
Ts and tb, characterizing system and environment respectively, a Markovian approxima- 
tion amounts to t$/tb — ► +00, that is either ts — ► +00 or tb — > 0. In the physical context 
we are actually considering, the typical time of the system is of the order of the period of 
oscillations, T5 ~ E/Am 2 . On the other hand, the characteristic time of the environment, 
tb, is related to the decay time of its correlations. To get a Markovian approximation, 
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we consider the case t b — ► 0; that is, we consider environments whose correlation func- 
tions decay in sufficiently short times. As we have seen in Section 12. 4( this procedure 
corresponds to the singular coupling limit. 

Remark 4.6 We are considering a homogeneous thermal bath in a finite volume V, thus 
the Markovian approximation, that requires an infinite volume, is to be performed in the 
hypothesis of a volume V big enough so that the Poincare recurrencies can be neglected pQ. 

The singular coupling amounts to considering two-point correlation functions behaving 
like Dirac-deltas in time: 

h aal (x,x') = TT B [p B Vg(x,t)V§\x',t')}^c aa/ (x,x , )5(t-t > ). (4.23) 

According to the discussion of the previous chapters, the resulting Markovian evolution 
equation is: 



p S (t) = -i [H e s ff , PS (t)]+lYl jjc aa ,(x,x')([Vg(x,0)p s (t),V s ?\x',0)] + 



V 

+ [Vg(x , 0),p s (t)Vg'(x>, 0)])dxdx' 



(4.24) 



and the statistical operator ps lives in the neutrino Hilbert space #5, containing flavor 
and space degrees of freedom. We shall later try to reduce further the description onto the 
flavor degrees of freedom. In (|4.24|) . fig is the effective neutrino Hamiltonian obtained 
from Hs by taking into account the mean field term produced by the interaction with the 
environment: H s = Hs + Tr B (p B Hi), with Hj defined in (|4.1U|) and p B the thermal 
equilibrium state of the bath; see the discussion in Section fl. 31 It is the effective Hamilto- 
nian usually introduced in the description of the interactions of neutrinos in matter leading 
to the MSW effect pTll92ll93UlU8j . In the two-dimensional flavor space, the representation 
of H e / f over the physical states is 

jeff _ Aro^ ( -cos 26 sin 29 \ ^ ( V ee V e „ 

where 



H&s ~ AE { sm26 cos2# ) + { V* e ' 



V ee = (v e \ Tr B (p B H I )\v e ), V efJL = {v e \ Tr B {pbHi) \ u N , 



= (l/ftl Tr b(pbHi)Wh), = (z^| Ti B (p B H I )\v e 

are contributions of first order in Hi and V ee , S R, V^ e = V*^. 



(4.26) 



Remark 4.7 The off-diagonal entries vanish for flavor-preserving interactions. Indeed, in 
this case in (|4.11|) we have 

V B m »(x) = S lm V^(x), (4.27) 
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where l,m = e, ii. For instance, for neutrinos in solar matter, we find: 

V ee = V2G F (n e -^n^ , V^ = -^=G F n n , V efl = V fie = 0, (4.28) 

where n e and n n are the density of electrons, neutrons respectively, and Gp is the Fermi 
constant. 

The effective Hamiltonian is diagonalized by an orthogonal transformation U that 
connects the eigenstates of Hff to the physical states: 

U=( COS0 ™ S[nd Q m ), (4.29) 
y — sin v m cos 9 m J 

where {|^ m i)> Wm2)} are the eigenstates of H% and 9 rn is the mixing angle in matter. The 
Hamiltonian eigenvalues in the adopted relativistic approximation are e m \ = —Am^/AE 
and e m 2 = Am^/AE where Am^, = rr? m2 — m 2 ml and m m i, m m 2 are the neutrino masses 
in the medium. They are related to the equivalent quantities in vacuum; if V e ^ = V^ e = 0, 
as in the case of solar matter, we have |92| I93| 




Am 2 sin 26 



Am^ = yj A 2 + (Am 2 ) 2 - 2AAm 2 cos 29, (4.30b) 

where 



A = 2E(V ee - V). (4.31) 

Remark 4.8 The resonance in equation (|4.3Ua|) produces an enhancement of the oscillat- 
ing behavior that, combined with a suitable variation of the density of the crossed medium, 
generates the so-called MSW effect, actually the best candidate to the explanation of the 
problem of solar neutrinos |92[ l93j . 

To first order in the interaction (that is neglecting the dissipative contributions of 
equation (|4.24jl ). the transition probability in the uniform medium considered reads 

P Ve ^ (t) = sin 2 29 m sin 2 (4.32) 
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The dissipative effects are of second order in Hj; they result in Hamiltonian, damping 
and noise terms as results from rearranging equation (j4.24j) : 



p s (t) = -i [H e s ff , Ps (t)\ - dxdx% aa ,(x,J)[V s a (x,0)V s a '(x',0),ps(t)} + 

a,a' y 

jj dxdx'a aa ,(x,x"){V s ?(x,0)V s a \x / ,0),p s (t)}+ (4 33) 
a,a' y 

+ Z) // d ^' a ^(x, x*)Vg'(af, 0)p s (t)Vg(x, 0), 



a. a 



V 



where the new coefficients are defined as 



/oo 
h aa /(x,x';u)du; (4.34a) 

/•oo rO 

b aa >(x,x')= / h aa '(x,af;u)du — / h aa '(x,x';u)du. (4.34b) 



o 



According to Section 11.51 in equation (|4.33j) we recognize the typical form of the gen- 
erator of a Markov semigroup. 



Remark 4.9 Accordingly to Remark 11.171 in the first row of Q4.33JI we note that the 
vacuum Hamiltonian H$ is redefined by contributions of first order in the interaction, 
already accounted for in Hg , and of second order, determined by the generalized matrix 
[b aa '(x,x')}. 

Remark 4.10 Since the time-evolution (|4.33|) is obtained from a singular coupling, it is 
completely positive and the generalized matrix [a aa /(x, x')] is positive definite. 



4.4 The generator of the dissipative evolution 

We want to write the form of the generator of the dissipative evolution, eq. (|4.33|) in the 
flavor space 9)2 only; we thus have to get rid of the momentum degrees of freedom. We 
start by considering the expression (j4,22|) for the statistical operator p$- By rearranging 
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the summation indices, Q4.33J) results in the following equation for the entries of ps{t)' 
pij(f,s;t) =~i^2 [( H s ff ) ik Pkj(r,S]t) - p ik (r,s-,t)(H e s ff ) kj 



+ 



k 
-,2 



T7n ^ yZ { \ A imlm(r, q,k,q) + B im i m (f, q, k, q) pij(k, s; t) + 

frn %j J (4.35) 

+ [Aimjm(q, k, s, k) - B lmjm (q, k, s, k) p u (r, q; t)+ 
- 2Aiju(k, s, q, r)pim{q, k;t)\ , 



where 



A lrnl , m ,(k, k',q, ^) = Y,jj dxdx'a aa ,(x, ^e^-^e'^-^ S^%, fyS^tf, q), 



(4.36a) 
(4.36b) 

and a = (I, m, p), a' = (I', m', p'). 

Equation (|4.35f) is rather complex; however, under suitable conditions, it is possible 
to simplify it, reducing the effective Hilbert space to a two-dimensional one. The Hilbert 
space describing our system is S)s = -f)2®£>p> where f)2 is the two-dimensional flavor Hilbert 
space whereas Sj p is the momentum Hilbert space. The reduction to Sj2 is performed by 
adapting the projector technique developed in Section 11.31 We shall see that such a 
reduction will produce a nicer equation than (|4.35|) . under conditions that are not too 
restrictive. 

Remark 4.11 Usually the formalism of projectors is applied to eliminate the environment 
degrees of freedom, leaving only those of the system of interest. In the present case, 
however, the trace over the environment has alredy been performed, and an adapted 
version of the projection technique allows us to discard the momentum degrees of freedom. 

Let us start with the time-evolution (|4.35j) : for later convenience we shall adopt the 
compact notation 

Ps (t) = L[p s (t)] = (L H + L D )[p s (t)], (4.37) 

in which we have separated the Hamiltonian and the dissipative part of the genera- 
tor ()4.33|) . In particular, Lh\ps{P)} represents the first row of ()4.33|) while L£>[ps{t)} 
the latter two. 
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The state ps acts on f)2 &>£)p and we want to extract the time-evolution of the reduced 
statistical operator acting on f)2 only. We will denote this latter by p; it is obtained by 
tracing p$ over the momentum degrees of freedom: p = Tr p (ps). 

In order to apply the projector technique, we need a reference state for the momentum 
degrees of freedom: we choose p p = \p){p\, where p is the momentum of the incoming 
neutrino; accordingly we define the projectors P[-] = Tr p ( ■ )<8>Pp and Q = I— P. Adapting 
relations (|1.29j) and (|1.3()j) to this case, from (j4.37j) we extract the following equation of 
motion: 

p{t) = P{L H + L D )p{t) + / dt'P{L H + L D )e QL ^QL D p(t'), (4.38) 

J o 

where we have used the property QLhP = and supposed the initial state to be factorized: 
Ps {0) = p(0)<g)p p . 

This is not a manageable expression yet because of the convolution term in the right- 
hand side. In order to get a time-evolution having the structure of semigroup, for flavor 
degrees of freedom only, we have to make further assumptions. 

After explicit computation, the dependence of A(k,k' ,q,q') and B(k,k' , <f, <f ) on the 
momenta, in the operators PLd and QLq, is rather different: in PLq these coefficients 
present in any case two pairs of equal momentum indices, while this never happens for 
QLd- To understand the physical meaning of this, observe that these coefficients are 
Fourier transforms, for example 

A(U',m * /J (4.39) 
v 

where a(x, x') is a spatial correlation function of the bath; in general, it can be written as 
the sum of two functions: /(x, x') = f(\x — x'\), traslationally invariant, and a contribution 
g{x,x f ): 

a(x,x') = f{\x-x'\)+g{x^). (4.40) 

Therefore (|4.39|) becomes 

A{k,k',q^)^\^ w _ f) J dzf(\z\)e^-^ + J! dxdx'g(x^)e^'-^e-^-^' . 

V V 

(4.41) 

The contributions to PLd come from both the terms in the right-hand side of this 
equation whereas the contributions to QLd are produced only from the Fourier transform 
of g(x,x'). Translational invariance means momentum conservation, therefore the term 
PLd does not change momentum, QLd instead does. Thus, if the correlation function 
a(x,x') is traslationally invariant (i.e. g(x,x') = 0) automatically QLd = and the 
integral term in eq. (|4.38|) vanishes. In general it could be neglected if, in a(x,x'), the 
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traslationally invariant part f(\x — x'\) is dominant with respect to the non invariant one, 
g{x,x'). 

In the following, we work within the hypothesis of almost translational invariance 
within V, so that (|4.38|) becomes: 

p(t)=P(L H + L D )p(t), (4.42) 
and the evolution equation for the entries of the 2x2 statistical operator p reads 



Pij(t) 



e s n ) ik P kj {t)-p lk (t)[H[ 



eff 



k 
2 

^2 Z){[^(P) + B u (p)}p l3 (t) + [Ajip) - B lj (p)}p il (t) + 2j2A mj u(p)Plm(t)}, 



(4.43) 



where 



Airnji(p) — Aj m ji {jp, k,p,k), Aij(p) — Ai m j m (p); 
k m 

B ij(P) = ^2Y1 B imjm(P, k,P, k). 



(4.44a) 
(4.44b) 



k m 

In vector representation, using the physical basis for S)2, this equation assumes the form 

d 



dt 



\p(t)) = £\p(t)) = (H + V)\p(t)), 



(4.45) 



where 



H 



( \ 

h 3 -h 2 

-h 3 /ii 

V h 2 -h x J 



V 



( \ 

u a b c 

v b a (3 

\ w c (3 7 ) 



(4.46) 



The matrix C = TC + T> is the generator of the reduced time-evolution we are looking 
for. The parameters characterizing the Hamiltonian part are 



Am 2 



sin 26 + Re V e u + oj\ , 



AE 

h 2 = -ImV^ +uj 2 , 

Am 2 1 

h 3 = cos 20 + - (V ee 

4E 2 K 



(4.47) 



Vaa) + UJ 3 , 
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where uj\ , u>2 and lu^ are the Hamiltonian contributions of second order in the interaction 
(see Remark 14.9)1 . whose expressions are reported in Appendix IDl 

Following the discussion of Sections 11.81 and 12.31 the coefficients u, v and w in the 
dissipative part T> are related to entropy behavior and to the existence of asymptotic 
states different from the maximally mixed state \p2) = 1/2(1,0,0,0). The remaining 
parameters, instead, produce decoherence and dissipation. The expressions of all of these 
parameters are listed in Appendix IDl 

These coefficients satisfy the inequalities ()2.21j) that guarantee the complete positivity 
of the evolution because they are obtained via a singular coupling limit |.S9j . 

In the neutrino's case the impact of complete positivity cannot be fully appreciated 
since this property is needed for the physical consistency of the time-evolution of bipartite 
entangled systems, and there are not phenomenological contexts of this kind in neutrino 
physics. 



Remark 4.12 If the interactions are flavor-preserving most of these coefficients vanish 
since we have 

Aimji — fiim&jiAnjj, Aij — 5ijAa, Bij — 8ijBa. (4.48) 

Considering the expressions reported in Appendix[Dj we find u>\ = L02 = 0, u = v = w = 0, 
a = a and b = c = 7 = 0; the dissipative contribution is fully parameterized by only one 
coefficient, a, in agreement with [9"5] . 



4.5 An example of flavor changing dissipative evolution 

As already mentioned (see Section f2.4j) . the exact integration of equation (|4.45|) is always 
possible, but the expression of Qt is not transparent enough to be amenable to physical 
interpretations. It is more convenient to expand Qt to first order in the dissipative pa- 
rameters appearing in T>; this computation has already been motivated and developed in 
Section |2~H and applied to neutrinos in |117U118j . 

In this section, however, we wish to give the exact transition probability in a simplified 
case |117) : /i2 = 0, u = v = w = 0, a = a = 7 and c = 0. As explained in Remark l4.121 this 
latter case amounts to flavor-changing interactions since b, (5 and 7 are not vanishing and 
it is compatible with the inequalities expressing the complete positivity of the evolution 
provided that a 2 ^ b 2 + (3 2 . 

Remark 4.13 In our vector representation the transition probability can be expressed as 
an inner product: 

V v ^(t) = 2{p^\p Ue (t)) = 2(p^\Q t \p Ue ), (4.49) 

as can be seen by putting in ()4.9|) the expansions of both p Ve and p v over the basis of 
Pauli matrices plus identity. 



so 
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After exponentiation of the generator C, the time-evolutor Qt = e a amounts to 

( 1 \ 


Qt 

Vo 

where the 3x3 sub-matrix Qt is given by 



-lot 



' sin 29 m t 

l3 + ^e — • M + 



sin 2 9 m t 2 
402 yW 



with 



(4.50) 



(4.51) 



6 + /i 3 
6-/13 P + hi 
/3 - fci 



(4.52) 



and I3 the 3x3 identity matrix. We have introduced O m = y/h'f + ti£ — b 2 — (3 2 , that can 
be either real or complex, depending on the relative size of the coefficients appearing in 
it. Since in the present context the dissipative parameters are small, we consider a real 
m ; the general case is presented in |117| . 

Remark 4.14 Note that, that, in absence of dissipation (i.e. b = /3 = 0) Q m reduces to 
Am^/AE, and then, in absence of matter, to Am 2 /4E. 

In vectorial notation, the pure v e state (the initial state) and the pure Vn state are 

1 



0,0,1), 
^(1,0,0,-1) 



whence 



/ N 1 



1 



-2at 



P 2 



hi 



-2at 



Q2 



sin m i 



(4.53a) 
(4.53b) 

(4.54) 



Consistently, if there is not dissipation Q4.54JI gives the transition probability in matter 
to first order in the interaction, that is equation (|4.32[) . In absence of matter, we get the 
vacuum transition probability (14. 7|) , 

The damping exponentials, typical of dissipative phenomena, flatten the oscillations 
in the transition probability so that it approaches the value 1/2 for very long times, and 
this happens regardless of the value of the vacuum mixing angle 9. We remind that in 
vacuum the transition probability has the upper limit sin 2 26 and then tends to vanish for 
9 approaching whereas the asymptotic mixing between the two flavors, in the dissipative 
case, is maximal, regardless of the value of 9. 
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4.6 Non homogeneous matter 

Up to now we have considered a simplified environment with a constant mean density. 
However the physical appeal of neutrinos interacting with matter arises when the density 
of the considered medium is varying along the neutrino path. This is indeed the real 
situation: in the solar case the neutrinos are produced near the core of the Sun, the region 
with higher density, then they travel through regions with gradually decreasing densities. 
Actually the more plausible explanation of the missing neutrino solar flux depends strongly 
on the varying profile of the Sun density |92U93| . 

Since any slab of matter with a not-constant density profile can be well approximated 
by a sequence of homogeneous layers with suitable thickness, the time evolutor can be 
obtained by composing elementary evolutors computed using the formalism described 
before. That is, if the inhomogeneous medium consists of n homogeneous layers, one has 

& = 0£o...o0£ 00*, t = h+t 2 + ... + t n (4.55) 

where ti, i = 1, . . . n are the times needed to traverse the successive layers. The composition 
of maps in 1)4.55(1 is completely positive by construction. 

Then in our formalism all we need to describe the neutrino evolution in the real mat- 
ter is to divide the traversed medium into suitable slabs, possibly of different thickness, 
depending on the rapidity of the mean density change. 

Traditionally, in the study of neutrino transitions in matter two approaches are possible 
to account for the change of the medium density: it can vary slowly, and the transitions 
are called adiabatic, or rapidly. In the adiabatic approximation, neutrino state evolves in 
time like eigenstates of the effective Hamiltonian that, in the neutrino frame, is varying in 
time; otherwise neutrinos do not necessarily remain in eigenstates of H e J^{t) but can pass 
from an eigenstate to another. This standard formalism has been adapted to the case of 
a dissipative evolution for statistical operators in |117j . 

4.7 The solar matter case 

In this section we apply our formalism to describe the dissipative time-evolution of a 
neutrino crossing a layer of solar matter with constant mean density and to evaluate the 
impact of this modified evolution on the flavors transition probability. We adopt the 
Hamiltonian density (|4.14jl . 

The autocorrelation functions reflect the decoupling of flavors: 

ai m ^i'm'fi'(x,x) = Si m 5ii m iSiiiai m fj,,i'm' fi'(x,x'), 
bim^,i'm'^'(x,x') = 5i m Si' m 'Su/bi mfl: ii m i^(x,x'); 

5u' means statistical independence of the bath operators of different flavors. As a conse- 
quence of (|4.56|) most of the coefficients in Appendix El vanish and the evolution equation 



(4.56) 
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in the flavor basis, eq. Q4.45|) . takes the simple form 



where, following (|4.4 



/ \ 

h 3 

-h 3 hi 

\ 4i / 

4.28jl and Appendix [D 



/ \ 

a 

a 

\ / 



\p(t)), 



(4.57) 



Am 2 
G F 



V2 



Am 2 
AE 



cos 20 I + W3 



(4.58) 



w 3 = ^ Im [ 5 w(^> - 5 ee(p)] , 
Q= 2 [Aee(p) + A m {0)} . 

In order to evaluate the time-evolution of this state, we use the formalism developed 
in Section 12.41 with r\ = — cot # m and ri = — tan# m . It is convenient to perform the 
computation in the basis {\v m \), |^ m 2)} that diagonalizes H§ (introduced in (|4.29|) ): 
indeed P Ue ^ Uti (t) is basis independent by definition. 

The resulting expression of the transition probability in the dissipative regime is thus 



Xt) = \{l~e 

- sin 2 2(9, 



2asin 2 2e micos 2 



-a(l+cos 2 29 m )t Am 2 ^ 

2E 



2E 
Am 2 



2 Am 2 ] i 

,1 + 3 cos 2t/ m ) sin - — t >, 



2E 



(4.59) 



where the parameter uj 3 has been absorbed in a redefinition of the parameters appearing 
in the Hamiltonian part of the evolution. For a = we obtain the transition probability 
at first order in the interaction for an homogeneous medium, eq. (|4.32|) . It reduces to the 
vacuum oscillation probability in absence of matter. 

In order to obtain the transition probability in the real, not uniform, solar matter, 
the considerations of Section 14.61 should be applied; however, using a magnitude estimate 
based on dimensional analysis, the impact of dissipation for the phenomena that can be 
described in our framework is negligeable, as we shall shortly explain. 



4.8 Discussion and conclusions 



We have treated the topic of neutrino oscillations in matter from the rigorous point of view 
of open quantum systems, using rigorous Markovian limits that produce a fully consistent 
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dynamics for the described system and avoids problems of physical consistency. 

The interaction between neutrinos and environment has been considered to be as gen- 
eral as possible, including some scenarios forbidden in the Standard Model, like flavor 
changing currents. The formalism of the singular coupling has been employed, producing 
a request on the bath correlation time: it must be much shorter than the characteristic 
time of the system. 

The conditions to obtain a convolutionless dynamics in the flavor space have been 
discussed, showing the tight connection between the semigroup structure of the reduced 
dynamics and the isotropy of the bath correlation functions (i.e. their invariance under 
translations). No specific models of fluctuations have been introduced; we have rather 
assumed that the contributions breaking this isotropy be negligeable. If they are not so, 
flavor and momentum degrees of freedom need to be treated together and an effective 
two-dimensional prescription is unavailable. 

In this framework, several dissipative phenomena can be described. We have evalu- 
ated the transition probability between two neutrino flavors in two different contexts: a 
flavor-changing and a flavor-preserving interaction and explicit results for homogeneous 
matter have been obtained. To account for matter inhomogeneities, in general, numerical 
computations should be performed. 

The typical terms introduced by dissipation are damping exponentials whose lifetimes 
are the reciprocals of the dissipative coefficients; this means that the effect of the irre- 
versible dynamics should be relevant on a time scale of order tdi ss ~ ^/(9 n ) 2 TB, the inverse 
of the typical magnitude of dissipative parameters in non-relativistic matter, where n is 
the density of scatterers and g the coupling constant in the interaction term. 

For example, the thermal fluctuations in the Sun satisfy the limits imposed by our 
formalism; a satisfactory evaluation of the correlation time is given by the inverse of 
the thermal energy at temperature T: tb = K/kT, where h and k are the Planck and 
Boltzmann constants respectively. The explicit evaluation for the Sun core (T ~ 10 K) 
gives tb ~ 10~ 21 s, much shorter than ts ~ E/Am 2 ~ 10~ 4 s (obtained considering the 
values E ~ 1 MeV and Am 2 ~ 10 -5 eV 2 ). However, neutrinos cross the Sun in a time so 
short that the thermal contributions are completely negligeable: tdi ss ~ 10 14 s (using the 
electron density n e ~ 10 26 cm -3 and Gf ~ 10~ 5 GeV~ 2 ), versus a crossing time of about 
3 seconds. 

Does it exist a physical instance in which the dissipative effects give a strong effect de- 
spite the weak interaction with the surroundings? It is difficult to imagine a real situation 
of this kind. However the gravitational thermal bath described in Section 13.21 could give 
such an environment: it should be everywhere present, and plausibly homogeneous and 
isotropic. In this case the bath energy is parametrized by the Plank mass, Mp ~ 10 19 GeV, 
and for a neutrino energy of order Eg ~ lMeV, we can roughly evaluate tdiss ~ 10 7 s, that 
is of order of the time needed for a neutrino to cross intergalactic distances. The asymp- 
totical behavior for the transition probability, with and without dissipation, should be 
very different. In the first case it would approach the 1/2 value; in the second one it would 
oscillate with a maximal amplitude of sin 2 28 (observe that in this scenario it is impossible 
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to distinguish between and 6 m ). The averaged transition probability would be certainly 
increased and its time dependence destroyed. In general, decoherence in neutrino physics 
has been studied in |TT7IITTgilT5nilT2T| . 

The last remark is about the discrimination of the Dirac or Majorana nature of neutrino 
fields |122U123| : although dissipative dynamics in general allow to distinguish between 
these two possibilities |117| . this is not the case in this chapter, since the microscopic 
interactions just presented are diagonal in the flavor neutrino fields. 



Chapter 5 

Neutron interferometry 



In this chapter we apply the formalism of open quantum systems to study neutron inter- 
ferometry. There are two main motivations for this kind of approach. First of all, physical 
setups that are subjected to irreversible evolutions are actually available: they give the 
opportunity to study several models of Markovian dynamics, and to test their physical 
consistency. In the second place, since experiments based on quantum interferometry 
strongly rely upon the property of quantum correlations between different degrees of free- 
dom, they can be suitably employed to investigate the notion of complete positivity (CP) 
of the reduced time-evolution A^. As we have seen in Chapter^ this property, tightly 
connected to quantum correlations expressed by entanglement, has a somewhat abstract 
and technical definition; it could seem at first sight a mathematical artifact rather than 
a physically motivated request. However in this context the relevance of CP becomes 
apparent since we are dealing here with reduced dynamics of the form I ® At, precisely 
the one used in its definition (see Section IT3|) . Indeed, as explicitely shown in the models 
discussed below, lack of CP necessarily leads to the appearance of unacceptable negative 
eigenvalues in the time evolving physical state describing entangled degrees of freedom. 

After a brief introduction to the topic of neutron interferometry, and to the mathemat- 
ical framework usually employed to describe it, we shall study non standard dissipative 
evolutions for neutrons inside the interferometer. These reduced dynamics are the result 
of the interaction with an external environment in which the interferometric device is 
immersed. 

We shall be very general, allowing for evolutions that are not completely positive. 
Three different prototypes will be taken into account, all of them experimentally achievable 
by choosing appropriately the properties of the external environment (typically a thermal 
bath): they cover instances of completely positive, simply positive (but not CP) and even 
not-positive time-evolutions. Some of these dynamics are commonly used in some branches 
of physics (for example, in chemical physics |125U126ll2"7] ^. 

Neutron interferometry has been proposed to test the non-contextuality hypothesis 
in quantum mechanics |129j . that is the possibility that all the members of any set of 
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commuting observables posses a value independent on the choice of the commuting ob- 
servables measured at the same time (that is, independent on the context in which the 
measure of that observable is performed). This hypothesisis can be tested by studying 
violations of the CHSH inequality |128j (a particular Bell-like inequality); in the final part 
of the chapter, using the previously developed results, we will study how this relation is 
modified when the interferometer is immersed in a stochastic background. 

5.1 The physical apparatus 

For our discussion, we shall refer to an idealized, but standard model of neutron interfer- 
ometer 1 |129| . The apparatus, sketched in Figure I5TT1 consists of three basic components: 
an initial beam splitter (bs%), dividing an incoming, polarized, neutron beam in two com- 
ponents; a spin flipper (sf), acting on one of the two components emerging from the 
beam splitter; a final beam splitter (bs2) that recombines the two components, equipped 
with a phase shifter (ps). The neutrons of the incoming beam are polarized along a fixed 
direction, say the z-axis; the beam splitter bsi is characterized by a transmission p and 
reflection q coefficients satisfying the condition \p\ 2 + \q\ 2 = 1 (in the ideal case where no 
neutrons are absorbed). The two components follow spatially separated (by distances of 
the order of centimeters) trajectories, denoted by u (for up) and d (for down). Then, the 
translational degree of freedom may be described by a two-dimensional Hilbert space 9)\ 
generated by the orthogonal states \tp u ) and \ipd), that are associated to the two macro- 
scopically distinguished paths. 

Remark 5.1 Since the neutrons are 1/2-spin particles, their spin degree of freedom is 
also described by a two-dimensional Hilbert space Sj\, spanned by the states | \ z ), spin- 
up, and | iz), spin-down along the z-axis. On the whole, neutrons in the apparatus may 
be appropriately described by the 4-dimensional tensor Hilbert space Sj?, (8) Sfy — #4- 

Inside the interferometer, the component following the u trajectory is subjected to a 
spin- flip, | | z ) — ► | [ z ); as a result, the state describing a neutron travelling in the inner 
region of the interferometer is entangled (where the entanglement relies upon a quantum 
correlation between internal and spatial degrees of freedom of the same system): 



In the density matrix formalism, this state is represented by the 4x4 statistical operator 



M =p\^u) ® I iz) + <#<*} ® | T*). 



(5.1) 



p ab = IVXV'I = \p\ 2 Pl ®Q2 + \q\ 2 P2 ® Qi +pq*P 3 ® Qa+P*qPa O Q 3 , 



(5.2) 



similar set-ups can be constructed using polarized hotons sec 130 
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Figure 5.1: Typical model of neutron interferometer for non-contextuality tests, bsi and bs2 are the two 
beam splitters, sf is the spin nipper, ps the phase shifter, ci and C2 the position plus spin counters. 

where 

Pi = iV'uX^ui, p 2 = m^ d \, P3 = \^u)(M, p* = \iPd)(M, (53) 

Qi = \h)(h\, Q2 = \iz)(Ul Q 3 =IT*>U*|, Q4 = \iz)(1z\- 

Remark 5.2 We stress that the entanglement in the state (|5.1|) is different from the 
usual one, for instance the one described in the case of neutral B mesons, where quantum 
correlations connect the degrees of freedom of two spatially separated physical systems, a 
neutral B meson and its antiparticle, both produced in the decay of a resonance. In the 
present case, on the contrary, the entanglement is between different degrees of freedom of 
the same physical system: translational and spin degrees of freedom. 

By denoting with A^ the dynamical map describing the time-evolution undergone inside 
the apparatus, at the second beam splitter, after a time t, the initial state p ab in Q5.2|) is 
changed according to a law p ab — > p ab (t) = A t [p ab ]. 

At the second beam splitter, the u-component is further shifted by an angle <p, with 
respect to the d-component, by the phase shifter ps. The effects of the whole operation 
bs2 + ps are described by a unitary transformation 

m -> = U(6,<p)\1> u ), , , ( e-^ S m9 e-^cosfl 

where cos 2 9 and sin 2 8 are the reflection and transmission probabilities for bs2 (the ana- 
logue of q and p); their values can be modified by changing the angle 6, that is by changing 
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the physical properties of the beam splitter bs2 (the parameters p and q, characterizing 
the first beam splitter, are instead kept fixed). The initial statistical operator thus evolves 
as follows: 

p ab (U(0, tp) h)p ab (t)(U\6, if) ® I 2 ), (5.5) 

where 1 2 is the 2x2 identity in the translational Hilbert space 

The exiting neutron beam in the state 1)5 .5|) then impinges on two counters (ci and 
C2), equipped with spin-analyzers that can be set along any direction n = (n x ,n y ,n z ) in 
space. These counters are used to measure the relative frequencies with which neutrons are 
detected in one of the two exiting beams, with spin up or down along n. These frequencies 
of counts, denoted by Ol' n (9, (p) are given by mean values of the projectors Pj(0, (p) (8) Qh 
in $j% <g>$j\ over p ab (t): 

Of \6, p) = Tr[p ab (1)^(6, p) ® Qft], (5.6) 

where 

p 3 {e,p) = u\e,p)p 3 u{e,p), (5.7) 

with j = 1,2, are the spatial projectors in (|5.3j) . whereas Qn = | TnKtn I is the projector 
over the eigenspace spanned by the state | with spin-up along the n-direction. 

Remark 5.3 The last beam splitter together with the counters plus spin analyzers con- 
stitutes the measuring device, since it sets the values of the angles 9 and p, that can be 
varied by the experimenter. 

Neutron interferometry has proved to be an extremely powerful tool to investigate 
gravitational, inertial and phase shifting effects occurring inside the interferometer (1311 
I137j : they are described by a time-evolution for the spatial density matrix. Instead, in 
what follows, we shall assume a dissipative time-evolution to affect the spin degrees of 
freedom, because of the presence of a stochastic magnetic field inside the apparatus [39} 
I125j . The neutron beam thus behaves like an open system with time-evolution given by 
At = 1 2 (8> A b , A b being the dynamical map describing the non-trivial evolution in the spin 
sector. Explicitely, while traversing the interferometer, the neutron beam state changes 
according to p ab — ► p ab (t) = I2 (S> A b [p ab ]. Neutron interferometry can then be used to 
investigate the property of complete positivity. In particular, the form of At is the one 
appearing in the definition of complete positivity in Section ll.2l to guarantee the positivity 
of p ab {t) = At[p ab (0)] for any time t we need necessarily a completely positive map A b . 

In Section 15.41 we shall show explicitely that lack of complete positivity generates 
inconsistencies in the physical interpretation of the statistical operator p ab (t), that is the 
appeareance of negative eigenvalues. Noticeably, these eigenvalues can be directly accessed 
in experiments by setting appropriately the angles 9 and p and the direction n. Indeed, the 
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entries of p ab (t) can be expressed by means of the frequencies measured at the counters. 
Writing p ab {t) = £? 7=1 Pij{t)Pi <8> Qj, with P and Qj as in (j5,3[) . it turns out that 



P n(t) = Ol' z (0,0), p 12 (t) = <#-*((), 0), 
p 2 i(t) = 0^(0,0), p 22 (t) = O, 2 '- 2 (0,0); 

the remaining 12 entries pjj cannot be evaluated directly since in these cases the P, £§> Qj 
operators are not Hermitian; nevertheless they can be written by using suitable combina- 
tions of observables involving the following projectors, see (|5.7jl : 



Q ±x = i(| t*) ± | i*»«T* I ± (I* I), Q±, = 1(| T*> ± i| U)«t* I T 

Then, from 



(5.9) 



P 3 = -[(P+-P-)+i(P+i-P-i)] = PI Qs = -l(Q x -Q- x )+i(Q y -Q-y)} = Ql (5.10) 

the expressions of the remaining entries can be explicitely evaluated (see Appendix [EJ); for 
example 
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(5.11) 



Also these entries are thus, in principle, accessible to experimental measures. 



5.2 The stochastic thermal bath 



We assume that, inside the interferometer, the neutron beam is affected by a weak stochas- 
tic magnetic field; as a result, the dynamics is that of a dissipative open system. 
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We use the formalism introduced in Section II. 71 namely, the spin degree of freedom 
is described by a 2 x 2 stochastic density matrix satisfying the time-dependent 

Liouville-Von Neumann equation 

±m = v*+Lmmt)], (5 , 2) 

L„H = -i[H ,-]; L,{t)\] = -i[H,{t),-} 

where Hq = uj$o z /2; Hj(t) = Ylk=i Vk(t)&k is the fluctuating part of the Hamiltonian and 
Vk(t), k = 1, ... 3, are proportional to Gaussian stochastic magnetic field components. We 
assume (Vfc(t)) = V/c and the two-point correlation functions to be real, stationary and 
symmetric, that is 

hij(t -s) = (Vi(t)Vj(s)) = h* 3 {t -s) = hjiis - t). (5.13) 

The covariance matrix, defined as W(t) = [hij(t)], is thus real, stationary and symmetric. 

The effective 2x2 spin density matrix of the neutron beam when it impinges upon the 
second beam splitter is given by p{t) = (£(t)), where an average over the noise has been 
performed; its equation of motion has the form 

3 

p{t) = -i[H ,p{t)} + d kj (t)[a k , [aj,p(t)]] (5.14) 
k,j=i 

with 

3 nt 

dkj{t) = h k i(s)uij(-s)ds (5.15) 

l=i Jo 

where the Uij are the entries of the orthogonal matrix 



U{t) = [uij(t)\ = ( smujQt cosLOot | (5.16) 

such that <ji(—s) = Yl^=i u ij(~ s ) (J j- Equation 1)5.1 4|) is the form assumed by the general 
equation ()1.70|) in the two-dimensional case and we refer to Section ll.7l for details. The 
only assumptions required to extract the reduced dynamics 1)5.14)1 are the Gaussian nature 
of the stochastic magnetic field and the initial decoupling between spin and noise, that is 
p(0) = <E(0)) = S(0). 

Our assumptions on the correlation functions make the 3x3 matrix D(t) = [<%(£)] 
real; thus it can be written as the sum of a symmetric and a skew-symmetric part. Defining 
the two 3x3 matrices 

D SA {t) = 1 [D(t) ± D T {t)] = [d S,A {t)l (517) 
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where T means transposition, the dissipative part in equation (|5.14|) can be written as the 
sum of a Hamiltonian contribution plus a purely dissipative one. Then 



Ho ~ Yl dfj(t)tijk<Tk,p(t) +2 niWjiPi*)} - °jp(t)tri • (5-18) 



ij,k=l 



In the limit of weak stochastic magnetic field, it is possible to perform a Markovian 
approximation for the reduced dynamics Q5.18JI as discussed in Chapter ^ in practice it is 
usually performed by extending to +oo of the upper limit of the integral in (|5.15jl . The 
matrices D s and D A become then time-independent and we obtain dynamical equations 
of motion in Lindblad form: 



where C 



3 i 1 3 / 1 

H 0~ Y d ij e ijk°k, P(t) +~ Y Ci 3 ( a iPi t ) Cf j ~ TjWj^iiPit)} ) = ( 5 - 19 ) 



i,j,k=l 



Remark 5.4 Following the result of Theorem ll.5| the dynamics (|5.19j) is completely 
positive if and only if the matrix C = [cij] is positive definite. In general, the brute 
force Markovian approximation adopted above may provide a Markovian equation (|5.19|) 
that generates time-evolutions completely positive, simply positive or even not-positive, 
depending on the behavior of the correlation functions, in (|5.15j) . 



Equation (|5.19j) can be rewritten using the vector representation \p(t)) for density 
matrices p(t) presented in Section 12.11 The Schrodinger-like equation that follows is: 



±\p(t)) = (H + V)\p(t)) 



(5.20) 



where H and T> are 4x4 real matrices representing the Hamiltonian, respectively dissi- 
pative contribution in (|5.19j) : 



n 














\ 










-h 2 







-h 3 





hi 




\o 


h 2 


-hi 








where 



V 



( \ 

a b c 

b a j3 

V c p 7 / 



(5.21) 



hi = 2d 



32- 



1d A 



21 



(5.22) 
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and 



2(4^ + 40, 


6 = 2df 2 , 


■2(d?i + 4»), 


/3 = 2df 3 , 


2(df 1+ df 2 ), 


c = 2df 3 . 



(5.23) 

7 

Notice that, in this parametrization 

(a + 7 -a -26 -2c 
-26 a + 7-a -2/3 (5.24) 
-2c -2/3 a + a-7 

Prom Theorems 12.41 and 12.51 it immediately follows that 

1. The reduced dynamics generated by (|5.19j) is positive if and only if the parameters 
entering in T> satisfy the inequalities: 

a > [ a))!) 2 

a ^ ; < 07 > c 2 ; 007 + 26c/3 - ac 2 - 7& 2 - a/3 2 > 0; (5.25) 
7^0 ( aa^ P 2 

2. The reduced dynamics generated by (|5.19j) is completely positive if and only if the 
parameters entering in T> satisfy the inequalities: 

2R = a + j-a^0 ( RS^b 2 

2S = a + -/-a^0 ;< RT > c 2 ; RST - 2bc(3- Sc 2 - R(3 2 - Tb 2 ^ (5.26) 
2T = a + a- 7 ^0 [ ST ^ (3 2 

that is, if C is positive. 



5.3 Three different stochastic backgrounds 

In order to investigate the notion of complete positivity in more detail, we consider different 
decaying properties of the covariance matrix W(i). We focus our attention on three 
particular choices, reproducing equations of motion usually employed in the study of open 
systems |125|ll26| l2T|: the aim of this work is to give again direct evidence to the fact 
that the absence of complete positivity in the reduced dynamics generates necessarily 
physical inconsistencies in the interpretation of the formalism. In particular, these reduced 
dynamics may be experimentally realized by suitable stochastic magnetic fields and thus 
describe concrete physical contexts. 
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5.3.1 White noise 

The first case is that with white noise correlations for the stochastic magnetic field: 



(5.27) 



h 



j i 



hij(t) = hijS(t), 

where hij (i, j = 1, . . . 3) are time-independent, such that h. 

The covariance matrix W(t) is thus time-independent, positive and symmetric. Con- 
sequently D A = and D s = — W ^ 0. Then C ^ and the evolution is automatically 
completely positive; the dissipative coefficients are related to the correlation functions by 



/)*. 



a = -2(h 2 2 + h 



33 I 



2h 



12^ 



a = -2(h u + h 33 ), f3 = 2h 23 , 
j = -2(h n + h 2 2), c = 2h 13 . 



(5.28) 



5.3.2 Diagonal covariance matrix 

In this case, we suppose W(t) diagonal with 

h ll (t)=h 22 (t)=g 2 B 2 1 e- x W 



h 33 {t)=g 2 Ble-^\ 



(5.29) 



where the B^s are the time-independent components of the magnetic field, g is a constant 
proportional to the neutron magnetic moment and A, \i are decay coefficients. 
We find 




C 




where 



2g 2 B( 2g 2 Bf\ 
a = — — i + 1 



A 2 + 



2 ' 



7 = 



4g 2 BfX 



2 ' 



The 4x4 matrices 7i and T>, in Q5.2UJI take the form 



n 



( o o o o \ 

-u 

to 

V o o o o / 



V 



( \ 

a 

a 

\ 7 / 



where u = uj + -j^^ 



(5.30) 



(5.31) 



(5.32) 
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Therefore, the entries of the spin density matrix p(t) satisfy the Redfield-Bloch [I] 
equations of motion: 



Pl(t) = ~7 Pi (*) + 7P2(t), P2(t) = -f>l(t), 

fa{t) = -iup3(t) ~ 2ap 3 (t), p 4 (t) = [p 3 (t)]* , 
that can be explicitely solved yielding 



(5.33) 



Pl (t) = \[{l + e- 2 *) pi(0) + (1 - e" 3 *) pa(0)] , 

p 2 (t) = l - [(1 - e- 2 ^) Pl (0) + (1 + e- 2 ^) p 2 (0)] , (5>34) 

P3 (t) = e-( 2a -^Vs(0), 

P4 (t) = e -( 2 «+^v 4 (o). 

Since a and 7 are real and positive, the positivity of the dynamics is assured; complete 
positivity requires C ^ and then 2a — 7 ^ 0. From (|5.34|) we observe that diagonal and 
off-diagonal elements of p decay with different relaxing times, T\ and T2 respectively; they 
are the reciprocals of the coefficients 27 and 2a. Then the condition of complete positivity 
reads I/T2 ^ l/2Ti (HUES]- In |125[I126| it is shown that this typical order relation 
can be reversed by setting B% = and keeping fourth-order terms in the general master 
equation (|1.67|) leading to (|5.14j) . In such a case we have a positive, but not completely 
positive time-evolution. 

5.3.3 Single component field correlation 

In this case the magnetic field is supposed to have only the x-component: V\(t) 7^ 0, 
V 2 {t) = V 3 (t) = with 



consequently 



with 



h u (t) = 


9 2 B 2 e 


-x\t\. 






( 1 






2a 


-2b 


-1 


:)■ 


H 


-2b 





\ 












2g 2 B 2 X 


b = - 








A 2 + u; 2 ' 




X 2 +col 







(5.35) 



DA = - 2 ( X T+l l) I "I I , C=\ -2b I (5.36) 



a=-r- — 6 = --5- — 2- ( 5 - 37 ) 
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In the vector representation (15.20(1 . the dissipative matrix has the form 



V 



( \ 

6 

6 a 

\ a ) 



(5.38) 



whereas the Hamiltonian contribution 7i, expressed as in equation 1)5.32 J) , is characterized 
by u = u>o + 2 ^ 2 ^J . The corresponding Redfield-Bloch equations are 



pi(t) = -api(t) + ap 2 (t), p 2 (t) = -pi(t), 

p 3 (t) = -iupz{t) - ap 3 (t) + ap±(t) + 2ibp 4 (t), p±(t) = [p 3 (t)\* , 
with solutions 



(5.39) 



Pl(t) = 


i[(i + e -*-) Pl (0) + 


(l-e- 2 -*)p 2 (0)], 




P2(t) = 




-e- 2at ) Pl (0) + 


{l + e- 2at )p 2 (0)], 




Ps(t) = 


e~ at 


OL 

P4(0) cosh St H 


— 2ib 

— — P4(0) sinh St + 




^-p 3 (0) sinh St 




Pi(t) = 


e~ at 


at 

p 3 (0) cosh 5t -\ 


p3 (0) sinh St — 




^-p4(0) sinh St 




= \J a 2 


+ 4b 2 


-UJ 2 . 







(5.40) 



In this case, positivity and complete positivity asks for 6 = (and thus uj = ujq), while 
for b the time-evolution is not even positive. 



5.4 Positivity and complete positivity 

The covariance matrices in Subsections 15.3.21 and 15.3.31 do not guarantee the complete 
positivity of the time evolution. In this section we address the problems that arise in 
these cases. 

In the diagonal covariance matrix approach, let us consider the maximally entangled 
state, p = —q = l/v2 in ()5.1)) : it evolves into 



p a \t) = g [ p i ® Q*® + p z® " p 3 ® Q 4 (t) - ^4 ® Q 3 (t)} (5.41) 
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where Qi(t) are computed from equations (|5.34|) . The 4x4 statistical operator has the 
form 



P a \t) 



( £?_(*) \ 

E+(t) F(t) 

F*(t) E + (t) 

\ E-{t) J 



where 



E ± (t) = -(l±e- 2 *), F(t) = 



1 -(2o-iw)t. 



the eigenvalues of p (t) at the exit of the interferometer are thus 



(5.42) 



(5.43) 



Ai,2(i) = ^(l-e- 2 ^), 



A 3 ,4(*) 



1 + e 



-2-yt 



± 2e" 2a *) . 



(5.44) 



Observe that A4(0) = and dX±(Qi)/dt = (2a — 7)/2; then, by continuity, if a < 7/2, 
there is whole range of t where Xn(t) < and p ab (t) loses its physical meaning. On the 
other hand, if the complete positivity condition is fulfilled, X^t) for any time t and 
no problem arises. 

Similar conclusions hold for the single component field correlation case. 



Remark 5.5 It is necessary an entangled initial state to see the problems generated by 
positive dynamics that are not completely positive. Indeed, suppose to have a factorized 
initial state, i.e. in the form 

y 

with Cjj 0, and let it evolve following the factorized dynamics I2®A^; then the statistical 
operator at time t is positive for any t, because of the positivity of A^: 

p ab (t) = (I 2 ® Ai)[p ab ] = Ci jP ? ® Ai[p b ] > 

Remark 5.6 Besides the appearance of negative eigenvalues, we observe that, in the sin- 
gle component field correlation case, for 5 > a, ps(t) and p±(t) diverge exponentially, 
see (|5.4U|) : the positivity (and complete positivity) condition, 6 = 0, avoids this unaccept- 
able behavior. 

Remark 5.7 Although apparently formal, these resuls are far from being academic. In- 
deed, as alredy stressed at the end of Section f5.1( the entries of p are in principle directly 
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Figure 5.2: Three different behaviors for Xi(t): for a > j/2 (green line), for a — 7/2 (blue line) and for 
a < 7/2 (red line). We see that in the latter case there is a whole range of t where A<i(i) < and then the 
density matrix loses its physical meaning; indeed, the time-evolution is not completely positive. 

accessible to the experiment. For example, in the diagonal covariance matrix case, it fol- 
lows from eq. £3J that = 0^(0,0) and E+(t) = O^ z {0,0) while F(t) coincides 
with the expression in equation 1)5.11(1 . 

By modulating a background magnetic field close to the stochastic properties investi- 
gated in the previous three cases, one might reproduce experimentally the conditions for 
three different reduced dynamics and check their consequences. 

Then, one may conclude that reduced, Markovian time evolutions A^ must be not only 
positive, but also completely positive, since the lack of any of these constraints results in 
experimentally detectable inconsistencies. 

The use of one reduced dynamics instead of another depends on the Markovian ap- 
proximation used to derive it and whether, given the properties of the stochastic field, it 
was justified or not. It thus seems appropriate to conclude that, whenever a semigroup 
composition law is expected, the physically appropriate Markovian approximations are 
those which lead to reduced dynamics consisting of completely positive maps |35) . 

5.5 Noncontextuality: the Clauser-Horne-Shimony-Holt in- 
equality 

Quantum Mechanics has been shown to be a non-local and contextual theory, as it violates 
the so-called Bell inequalities that have to be necessarily satisfyed by local, realistic and 
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noncontextual theories |127j . 

Neutron interferometry is an optimal tool for investigating noncontextuality |129j : 
indeed in this case the entanglement does not involve spatially separated parts, as in the 
usual tests exhibiting violations of Bell-like inequalities which focus on locality. 

It is customary to use a particular form of Bell inequality, due to Clauser, Home, 
Shimony and Holt |128j . Using the previously defined projectors in configuration and spin 
space, one constucts the observables 

A(e, ip ) = p 1 (e,ip)-p 2 (e,i P ), B(n) = Qh- Q-n- (5.45) 

These observables commute as translational and spin observables do; they are called 
dichotomic since their eigenvalues are ±1. Choosing two different values for both 9 and 
(p and two distinct directions for the versor n, we obtain the set of 4 observables A(6i, <pi) 
and B(fij) with i, j = 1,2. 
From the mean values 



C t (6,<p;n) = Tr[p ab (t)A(9,<p)®B(n)] = 

= 1 /{9, <p) + 0^- H (6, <p) - Ol>- H (9, <p) - 2 ' H (9, <p), 



with four possible configurations of the control parameters 9, (p and n, the Clauser-Horne- 
Shimony-Holt inequality is obtained |128j : 



|C t (0i, <pt; ni) + C t (9x, tpr, n 2 ) + C t (9 2 , cp 2 ; ni) - C t (9 2 , ^ 2 ;n 2 )| ^ 2. (5.47) 

This inequality holds if we suppose the two possible outcomes ±1 of a measurement of 
the observables A(9i, tpi) and B(fij) are predetermined by the state p ab (t) independently of 
which other commuting observable is simultaneously measured (for example, the value of 
^4(#i, ipi) is unchanged if we measure simultaneously either B(n{) or B{fi2)). A violation of 
this relation implies necessarily the contextuality of the underlying theory. Equation (|5.46|) 
shows that this analysis can be directly performed by means of the frequencies of counts 
in the interferometric apparatus. 

Remark 5.8 In the configuration described in |129|I130] the state p ab is not subject to 
any dynamics inside the interferometer, but in our context it does: this is the reason why 
we have considered an explicit time dependence in the Ct(9,(p,n) mean values. 

For the maximally entangled initial state (p = —q = l/\/2) these quantities take the form 



C t (9,<p;n) = sin 2 9 Tt b [Q 2 (t)B(n ) +cos 2 0Tr 6 Qi (t)B(n)} + 

(5.48) 

- sin 29 Re {e~ lLp Ti b [Q 4 (t)B(n)]) 
where Tr& is the trace over the spin degree of freedom. 
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In vectorial representation, the time-evolution in the spin sector reads \p(t)) = Qt\p(fi)}\ 
in the following we shall consider the first order approximation of Qt in the dissipative pa- 
rameters, supposed to be small with respect to the parameter characterizing the standard 
Hamiltonian evolution, loq (for the expression of the entries of Qt, see Section l2.4j) . We 
introduce two vectors, G(t) and F(t), whose components are defined by 



Gi{t) = {g t ) a , Fl {t) = {Q t ) tl -t{Q t ) l2 , i = l,2,3. 
Then, the traces in (|5.48j) can be written as: 

Tr b [Q h2 (t)B(n)\ = ±G(t) ■ n, Tv b [Q 4 (t)B(n)] = F(t) • n; 
whence the mean values (5.2) read 

C t (6,ip;ii) = n- \G(t) cos26 - Re (e~ iip F(t)) sin 26 



(5.49) 



(5.50) 



(5.51) 

We will now discuss the explicit behaviour of Ct(0,ip;n) in the three cases introduced 
in Section 1531 to leading order in the dissipative parameters. 



5.5.1 White noise 

Though analytic expressions of Qt are obtainable in the general case of a stochastic mag- 
netic field with white noise correlations, these are rather involved and scarcely illuminating. 
More conveniently, we consider the expression of the time evolutor Qt to first order in the 
parameters, appearing in the dissipative part V of (|5.21|) . assumed to be small. One 
calculates 



Gl (t) = -— sin — cos [— + 4>c), 

4\C\ . uj t . /wot \ (5.52) 

G 2 {t) = sin — sin — + <p c , 

too 2 V 2 I 



G 3 (t) = e~ 2 ^, 
where |C| 2 = c 2 + j3 2 and tan cj)c = (3/c\ 



\B\ 

Re le-^FUt)] = e -^ a+a)t cos(w t - <p) + — sin w * cos Up + 4>b), 

I B I 

Re [e~ iip F 2 {t)\ = e - {a+a)t sin(w t - ip) - — sinw t sin (ip - <p B ), (5.53) 



Re [e-*F 3 (t)] = sin ^ cos ( ^ - <p - ^) : 

LOQ Z \ Z J 
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where \B\ 2 = (a — a) 2 + Ab 2 and tan</>£ = 2b/ {a — a). It thus follows that, in leading order 
with respect to a, b, c, a, f3 and 7 



G t {v, ip;n) = n z | — — cos 29 sin — cos I — + (p c ) 



sin 26* 



-{a+a)t 



cos (wo* — p) H sinwot cos Up + <Pb) r + 

u 



c4\C\ 



+ n. 

- sin 29 



u> t . ( u t 

cos 20 sin sin I h 0c ) + 

wo 2 y 2 



-(a+a)t , 



1^1 1 1 

sin (u)Qt — (p) sinwot sin (tp — \ + 

+ nje- 2 ^ cos 26 + M sin 2 sin ^ cos f ^ - ^ - ^ ). 
I 2 \ 2 J J 



(5.54) 



5.5.2 Diagonal covariance matrix 

For the stochastic magnetic fields with correlation functions as in ()5.29j) we have 
G(t) = (0, 0, e~ 27 *) , F(t) = e -( 2a -Ht( l5 _j 9 o), 

whence 



(5.55) 



Ct(9,p;n) = n z e 27< cos2# — e 2ai sin26'[n :r cos (wi — p) + n y sin (u;r — tp)]. (5.56) 
5.5.3 Single component field correlation 

In the case of magnetic field with correlation functions as in (|5.35|) the computation is 
similar. One finds: 



G(t) = (0,0, e~ 2at ), 

F(t) = e~ at (cosh. St + ij sinh 5tj (1, -i, 0) + - + ^ %b sinh 6t(l,i,0) 
whence, assuming 5 > 0, one obtains: 



(5.57) 



C t (9,(p;n) 



n z e~ 2at cos 29+ 



-at 



sin 29 



if a ■ i r S \ w + 26 

cosh os + — sinh 01 cos 93 H — smh or sm 

0/0 



+ 



, . a . \ . u! — 2b 

cosh or + — smh or sm 93 H — smh or cos 99 

0/0 



}■ 



(5.58) 



5.6. DISCUSSION AND CONCLUSIONS 



101 



The lack of positivity preservation which characterizes the time-evolution leading to 
equation (|5.58[) manifests itself in that the quantities Ct(6, (p; n) diverge in t when 5 > a. 

Expressions ()5.54jl . (|5.56j) and (|5.58j) agree with those used when there is no dissipation, 
namely putting a = b = c = a = f3 = 'j = 0: 



Ct(9, ip; n) = —n x sin 29 cos (uot — ip) — n y sin 29 sin (uot — <p) + n z cos 29. 



(5.59) 



Notice that the unitary time-evolution generated by the Hamiltonian Hq contributes 
to a time- varying redefinition of the angle (p. 

Concerning the issue of complete positivity versus simple positivity, in expressions (|5.54j) 
and Q5.56JI the two possibilities result in different relaxation properties due to whether in- 
equalities ()5.26j) or (|5.25jl are fulfilled. No physical inconsistencies may affect the mean 
values Ct (9, ip; n); indeed, negative probabilities may result in negative mean values of pos- 
itive observables only if the latter are entangled. In the case of the quantities involved in 
inequality (|5.46|) . the observables are factorized, Pi t 2(9, l p)®Qn an d the positivity of their 
mean values is preserved even when the time evolution is only positive and not completely 
positive. 

This can be seen as follows. By definition of dual map A|, to the Schrodinger time- 
evolution p ab {t) = (I2 ® At) [p ab (0)], there corresponds the Heisenberg time-evolution of 
observables A{t) = (I 2 ® A$)[A], 



Tr 



(I 2 ®A t )\p ab (0)]A 



Tr 



a;)[A] 



(5.60) 



The maps A|, dual to At, form a semigroup of dynamical maps that transform positive 
observables into positive observables, if the A^'s preserve the positivity of states. Conse- 
quently, even when the initial state p ab is entangled and the A^'s positivity preserving, but 
not completely positive, it turns out that 



Tr 



p ab (t)P 1 , 2 (9,i P )^Q H 



Tr 



p ab P j (9,ip)^)A* t [Q H ] 



> 0. 



(5.61) 



5.6 Discussion and conclusions 



Complete positivity is a property of quantum time-evolutions which is enjoyed by the 
standard dynamics of closed quantum systems generated by Hamiltonian operators, but 
not automatically by the more general reduced dynamics describing time-evolution of 
open quantum systems in interaction with suitable environments. Complete positivity is 
intimately related to the phenomenon of quantum entanglement, between two different 
systems, but also between two different degrees of freedom of a same physical system. 

In this chapter we have considered the two entangled degrees of freedom, translational 
and rotational, of a beam of neutrons travelling through an interferometric apparatus. 
We have studied the consequences of placing the interferometer in a stochastic, gaussian 
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magnetic field weakly coupled to the spin degree of freedom, providing an experimentally 
controllable environment. As explained in the three cases of Section f5.3( the same Markov 
approximation naively yields a semigroup of dynamical maps 1 2 (8> Af , where only the spin 
degree of freedom evolves in time; by varying the decay properties of the external field 
correlations, these maps turn out to be alternatively completely positive, simply positive, 
not even positive. 

The noncontextuality tests proposed in [12911130] are based on the Clauser-Horne- 
Shimony-Holt inequality (|5.47|) without time-dependence. The presence of a fluctuating 
magnetic field induces relaxation on the spin degree of freedom with strength and prop- 
erties depending on those of the field. Typically, the mean values in the inequality are 
damped and make it more difficult to be violated. However, in presence of stochastic 
fields yielding reduced dynamics that do not preserve positivity, the inequality might be 
dramatically violated because of possible mean values diverging in time. 

This latter possibility is a manifestation of the fact that any physically consistent time- 
evolution A^ must preserve the positivity of spin states in order that the eigenvalues of 
the corresponding spin density matrices might at any time be used as probabilities, in 
agreement with the statistical interpretation of quantum mechanics. If A^ preserves the 
trace of spin density matrices, but not their positivity, spin states may evolve in time 
in such a way that some of their eigenvalues become negative, while others greater than 
1, without upper bounds. It is this physically unacceptable phenomenon that leads to 
diverging mean values. 

The request of positivity preservation by the maps A^ with respect to spin states is 
thus unexcapable, but it is not enough to avoid physical inconsistencies when the time- 
evolution maps I2 ® Aj act on states p ab with correlations between spin and translational 
degrees of freedom. 

Inequality (|5.47|) does reveal the difference between completely positive and simply 
positivity preserving A£, but only as long as the relaxation characteristic are concerned, 
without any further effect (as the divergence of some contributions to the inequality). In 
fact, the positive observables in Q5.46JI are factorized, that is they incorporate no entan- 
glement between the translational and spin degree of freedom. Even if the initial state 
does incorporate entanglement, it nevertheless follows that the mean values of factorized 
observables remain positive and bounded. 

However, the interferometric apparatus proposed in |129U13Uj might also be used to 
measure the entries of the states of the neutron beam at the exit of the interferometer. In 
this way, one might have access to the spectrum of an initially entangled state after being 
subjected to the effects of the stochastic magnetic field. 

In the case of fluctuating fields yielding reduced dynamics that preserve positivity, 
but are not completely positive, the theoretical predictions indicate the appearance of 
negative eigenvalues, that is of negative probabilities, in the spectrum of the entangled 
exiting state. The fact that they are, in line of principle, detectable experimentally, does 
not allow to dismiss such an occurrence as practically negligible. Rather, it forces to 
reconsider the Markov approximation used to derive the time-evolution and to select as 
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physically consistent only those providing completely positive reduced dynamics. 
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Appendix A 

Dissipative parameters in the 
eigenstates basis 

For an unstable system, the dissipative parameters entering into 1)2.46(1 are denned by: 



Ai = , ■ ^\-w(l + \r 1 \ 2 )(l-\r 2 \ 2 )-2(l + \r 1 \ 2 )Re(Ur 2 )\, 

Fi + Tiy V ) 

A 2 = ] l -—{- w {l + \r 2 \ 2 ){\- |r 1 | 2 ) + 2(l + |r 2 | 2 )Re(C/r 1 )), 

\r\ + r<i\ y J 

a = L^{-u;(l + |r 2 | 2 )(l-|r 2 | 2 )-2(l + |r 2 | 2 )Re(C/r 2 )), 

Fl + r 2\ > 

i = | j |-^(l + |r 1 | 2 )(l-|r 1 | 2 ) + 2(l + |r 1 | 2 )Re(^r 1 )|, 
Fi + r 2 \ z I J 

Si = 1 ^—^{w(l-r 1 r* 2 )(l-\r 2 \ 2 ) + 2(l-r 1 r* 2 )Re(Ur 2 )\, 

Fl + r 2\ ) 

$2 = ] l -^{-w(l + \r 2 \ 2 ){l + r 2 rl) + {l + \r 2 \ 2 ){U\l-Ur 2 )), 

Fl + r 2 \ z L J 

01 = | j |2 { W (l-r 1 r 2 *)(l-|r 1 | 2 )-2(l-r 1 r 2 *)Re(^r 1 )|, 

02 = , ^- 1I {- U ;(l + |ri| 2 )(l + r 2 ^) + (l + |r 1 | 2 )(i7M-C/r 2 )j, 

Fi + r 2\ J 

w = | ^^{w(l-rir 2 f )(l+r 2 r 1 *)-(l-r 1 ^)(C/*^-C/r 2 )), 

Fi + r 2 \ z I J 

= 1 l -—{ w (l-r 2 rl){l + r 2 rl)-(l-r 2 rl)(Wrl-Ur 2 )}, 
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A. DISSIPATIVE PARAMETERS IN THE EIGENSTATES BASIS 



A = i 1 —^\j(l-\r 1 \ 2 )(l-\r 2 \ 2 )-2ARe(r 2 r* 1 ) + 2Re(Br* 1 r* 2 ) + 

In + r 2 r 

-2(1 - |r 2 | 2 ) Re(CrJ) + 2(1 - |n| 2 ) Re(Cr 2 *)}, 
S = , ri T~|2 {t(1 ~ ln| 2 ) 2 + 2^|n| 2 - 2Re(Brf ) - 4(1 - |n| 2 ) Re(Cn*)}, 

5 = 1 ^{ 7 (l-|r 2 | 2 ) 2 + 2^|r 2 | 2 -2Re(i?r* 2 ) + 4(l-|r 2 | 2 )Re(Cr*)), 

In + t 2 \ m i J 

A = -. ^—^{- 1 (l-\r 2 \ 2 )(l+r 1 r* 2 )-A(\r 2 \ 2 -r 1 r* 2 ) + Br* 2 2 -B*r 1 r 2 + 

In + r 2 \ z i. 

+(C*n - Cr 2 *)(l - |r 2 | 2 ) - 2(1 + n r* 2 )Re(Cr* 2 )} , 

$ = | j |2 {- 7 (i-|n| 2 )(i + nr 2 *)-A(|n| 2 -n^) + i?V 2 -i?r 1 V 2 * + 
In + r 2 r I 

+(C*n - Cr 2 *)(l - |n| 2 ) + 2(1 + nr£)Re(Cri)}, 

= | j |2 {- 7 |l+nr 2 *| 2 -A(|n| 2 + |r 2 | 2 )-2Re(i?n*r 2 t ) + 
In + r 2 \ z i 

+2Re [C (n*(l-H 2 )-r 2 *(l-|n| 2 ))]}, 

6 = : l -^{- 1 {l + r 2 rl) 2 + 2Ar 2 rl + Br{ 2 + B*rl-2{l+r 2 rl){C*r 2 -Crl)\. 

In + ^2r J 

For a stable system, they are: 

Ai = [w(l-\r l \ 2 )-2Re{Ur l )}, 

A 2 = -{ w (l-|n| 2 )-2Re(C/n)}, 

o = ^2{ w (l-|n| 2 )-2Re(C/n)}, 

£ = -|n| 2 {^(i-|n| 2 )-2Re(C/n)}, 

5 2 = -{2 W -U*r* 1 + ^}, 

<fo = -\ ri \ 2 {2w-U*rl + ^y 
5\ = 4>\ = uj = 6 = 0, 
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A = (i + | ri |2)2 ~ l^il 2 ) 2 ~ 2A \^\ 2 + 2Re(Brf ) + 4(1 - |n| 2 ) Re(Crt)}, 

E = ~ (1 + IriP) 2 {" 7(1 ~ ln|2)2 " 2Alnl2 +2R < Br ^ + 4 ^ " \n\ 2 )MCrl)}, 

S = ~ Uj2(i + | ri |2)2 {-7(1 ~ ln| 2 ) 2 - 2^|rx| 2 + 2 Re(i?rf) + 4(1 - |n| 2 ) Re(Cr 1 *)}, 

A = , □ ,^ (2 7 (l-|r 1 | 2 ) 2 -A(l-|r 1 | 2 ) + 4-^V 2 + 



(l + lnl 2 ) 2 1 1 11 ' y 1 11 ' r\ 
(cVi-^)(l-|rx| 2 )-4Re(C^)}, 

7 T ^ FF {27(1 - Inl 2 ) 2 - A(l - N 2 ) + | - B*r\ + 



-(cV 1 -^)(l-|r 1 | 2 )-4Re(Cr 1 *)}, 
= (i + | ri |2)2 {~47|ri| 2 - ^(N 4 + 1) - 2Re(Brf) - 4(1 - | n | 2 ) Re^)}, 

Note that, in this case, the following relations hold: 

A 2 = -Ai, £ = -|n| 2 Ai, fo = \ri\ 2 5 2 , 
H = -h| 2 A, S = -^A, $ = -|n| 2 A. 
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Appendix B 

Expressions of the time evolutors 



The entries of the time evolutor Qf in (|2.30[) for an unstable system are 1 : 



{Q d t)n 






{G d t)i2 


= J_(S + a)(e-^-e- 


-71*), 




- A 2 r+ (A + ^)(e- r -- 




{Qi)u 


= (A* + 5*)(e -r+ * 


- e-^) 




= J_(H + 0(e-^-e- 


-7i* )) 








(Gt)23 


= J^($ + 1 )(e- 72 *- 







2 

AT, 


-($* + #)( e -^ e 


-r +t)) 


(Qihi = 


2 
AT, 


-(^ 2 -$*)(e- r -*-e 




{Qt)32 = 


2 

AT. 


-(5 2 -A*)(e-^ t -e 




(G?)33 = 








(#)s* = 


1 

Am 


(6 + 9)e~ rt sin Amt, 




(G%i = 


2 

Ar_ 


-(^-$)(e" r + 4 -e- 


71*), 



1 By definition, A = (Ai + A 2 )/2 
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B. EXPRESSIONS OF THE TIME EVOLUTORS 



{QtU = ^(S 2 *-A)(e-^- e - r +*), 
( ^ )43 = ^( e * + r ) e_r<sinA ^ 



For a stable system they take the form: 



(#)n = \{e 2M + l) + ^(e 2M ~l) 
(Gt)i2 = ^(e 2At "l) 
(GfU = -^A(e^-l) 

= ^(e 2At + l)-^(e 2A '-l) 

{GtU = -^$*(e- iAmt - 1) 
Am 

(^)si = -^(h-^)(e lAmt -l) 
(G?)32 = -^(<5 2 -A*)(e iAmt -l) 

(9t)u = -r— 6 sin Ami 
Am 

(Gthi = ^m-<s>)(e- iAmt -i) 

{G d t)A2 = ^(5* 2 - A)(e~ iAmt - I) 

(Qf) 43 = J—e*smAmt 
Am 

d\ (Q—iAm)t 



(G?)u = e 



Appendix C 

Neutral B mesons decays 



C.l Semileptonic decays 



7V(B°;t) = \M h \ 2 e- T { S mhSr^- 

+ cosh 5t 



\ r L\ 



Re Aft 
- Re y h + 



+ sinwr 



+ COS OJT 



\th\ 2 + \tl\ 2 
\th + r L \ 2 
lm(r H rl) , + 
+ r L \ z 

2 MTL^l e -(n'+A') 
- kif + r L \ 2 



AT 

2 



+ 



r 5 2 + 
1 



Re A' 



+ 



2Am 



ReG' 



+ 



Mi 

r 



Re y h - - 



rs 2 + iv 2 



Re A' 



7V(B°;t) = \M h \ 2 e~ T ^smh5T^-Refj, h 

+ cosh 5t 
+ sin ojt 

+ COS LOT 



2 W T A/ A 

— n n Im A' - — — 

r 5 2 + Ar 



+ 



kjj + r L | 2 



+ Re Vh 



+ 



2 6 1 
Im/ift + — o ImA' + — — Re6' 



r s 2 + 
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-(fi'+A')f 



2Am 
Reyh 



+ 
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7\+(B°;r 



+ cosh 5t 



r L \ 2 ~ \r H \ 2 „ A' 

■ 12 Refi h - — 

L \r H + r L \ z AT 

\r H \ 2 + \r L \ 2 , D 2 



+ 



+ sihwt 



+ COS LOT 



lm{r L r* H ) 
2- ■ — % - Imnh 

\th + tjX 

2 Re(r g r£) 



r <5 2 + w 2 

ReG' 



2Am 



Re A' 



+ 



+ 



e-(^ + AOf +Re2/7i+ 2_|_ ReA ,]| ; 



P h -(B ;r) = |M fe | 2 e- r sinh^r 



2 w . A' 

ReA/l + r^T^ ImA "AT 



+ 



+ cosh <5r 
+ sin wr 
+ cos wr 



+ 



Im Xt 



2 <5 



r 5 2 + w : 



ImA' + 



2Am 



ReG' 



+ 



-2- pre 

\rH + tl\ 



(n'+A')£ 



+ Re y/t 



C.2 Decay in a CP eigenstate / 



^\M h \ 2 e~ T {sinh 5t 



|l + A{| 2 |r L | 2 



r^ J 5 2 +lo- 



. i 4A' 4 
ReA' - — - - 



-+ 



+ cosh <5r 



Ar r <5 2 + a; 2 

f 

L 



ImA' 



+ 



|1- A^| 2 |r^| 2 + |l + A{| 2 |r L | 2 



LO 



+ sm LOT 
4 



lm(r L r* H (l + \{)(l-\ f H *)) . 



a; 



5 2 +w 



ReA' + 



4 <5 



cos Loir 
4 



Y5 2 +lo 2 



ImA' 



+ 
+ 

(n'+A')f + 
+ 



1 Re(7- J ,r^(l + Aj)(l - Ag» e _ ((y+A , )g | 



r s/ j 2 +w 2 



l^ij + r L | 2 
4 5 

ImA'--- .ReA' 

T5 2 + io 2 
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V f (B°;r) 



^\M h \ 2 e 



sinh 5t 



2\r L r H \ 2 
6 



1-A 



/ 12 



H + a{| 2 



+ cosh 5t 



2\r L r H \' 



T, A/ 4A ' 4 W 

'1-A£| 2 + |1 + A{| 2 



ImA' 



Re A' 



■ImA' 

+ 
+ 



+ 



+ sin ujt 

+ COS LOT 



-4\r L r H \ 



2 Im((l-A^)(l + A{*)) 



5 2 +uj'- 



A\r L r H \ 



Re A' 



Vh + rL\ 2 
4 5 



-(sr+A')f 



ImA' 



2 Re((l-Ai)(l + A{*)) 



+ 

e -(Q'+A')f + 



5 2 +lo> 2 



ImA' + 



r^ 2 + 



Re A' 
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Appendix D 



Dissipative parameters in the 
neutrino case 

g 2 
g 2 

^2 = -^72 Im \ 2B ^if + A eefie (f - A efllifl (p) - A eeei Jf) + A^^ip)} , 

g 2 

^ = - Im i B ee if ~ B W if ~ 2A ee^ if] , 

g 2 

u = - Re \A e>iee if) + A mile if - A eiM if] , 

g 2 

V = y2 Im [Ae^eeif + A^^if - A efl (f)] , 

g 2 

C — 2,V^ [^ee/xe(jO Ae^^^if -\- A eee ^ 1 

g 2 
g 2 

P = - Im [^4 eeeA1 if - A^^ if - A eefMe if + A e/W if] , 

g 2 
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Appendix E 

Frequencies of counts 



P2i{t) 
P3l(t) 
P32(t) 
Pil(t) 
042 (*) 



\{[0]' x (0,0) - O^ x '(0,0)] +i[ol'~ y (0,0) - o\' v (0,0)]} 
i{[Oi' x (0,0) - O^^O)] -i[ol~ y (0,0) - o\' v (0,0)]} 
O 2 ' x (0, 0) - 0]:~ x (0, 0) + i 2 '~ y (0, 0) - 0} v (0, 0) } 
2 t ' x (0, 0) - o}~ x (0, 0)] - i \o 2 '~ v (0, 0) - 2 t ' y (0, 0)] } 

-4^(;.-f)-^(H)]} 



<#-'(£,o)-tf-(£,o 



;.-f!L,_?Uo, 2 --«^,|)]} 



o;"(j.o)-o t 2 '«(j,o 



+ i 



4' 2 



4' 2 



4' 2, 



P33(*) 



O 



1,3 



+ 



o 



2..r 



1,2/ 







2.;/ 



7T 



7T 

4'' 

7T 

4' 

7T 

4' 



7T 



1 4 



1,2/ 



+ i 

7T 

2 

7T 



7T 



2,2/ 



+ i 







O 



7T 

4 

1,3 
i 



2.x 



o 



1-2/ 



7T 



-,o 



+ 



7T 

4' 

7T 

4' 



7T 



-,o 

1, - 

2, - 



+ 

7T 

4' 

7T 



7P 

2, 

7T 



+ 
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4' 



O 



l..r 



"IT 



2,.r 



+ 



1..'/ 



o 



2,y 



O 



M-oi 



o 2 t 



TT 



7T 

4'° 

y (1 

U' 

y 



TT 

4' 



+ i 



dp 



TT 



2.y 



o 



+ i 



7T 
I' 



7T 

4'' 

7T 

4'' 



7T 

4'' 



+ 



+ 



O 



TT 

4'' 

7T 

4'' 



7T N 
7T > 



+ 



]} 



4' 



o 



1,1 



'7T 



+ 



2,.r 







7T 

2 

7T 

2 



7T 



.0 



.0 



o 2 t 



4 

7T 

4 

» (- 

U' 

4' 



+ i 



TT 



2.y 



+ i 



.0 



7T 

I' 1 

,l,a; 



f 

2..r 

t 



TT 
7T 

4'' 



O 



+ 



+ 

TT 

4' 

TT 
1' 



TT S 
2, 



+ 



!)]} 







l..r 



' TT 



o 



+ 



2,.r 



0; 



1,1/ 







2..V 







,0 )-oi 



o 



TT 



1° 
y [1 
U' 

y en 

4' 



dp 



+ i 



TT 



TT 



2.y 
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+ i 



TT 

4'° 
)2 „ 
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4' 
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4' 
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